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Abstract 

Anisoplanatism  in  adaptive  optics  (AO)  systems  is  a  performance-degrading  effect 
that  arises  whenever  light  from  the  wave-front  sensor  beacon  and  light  from  the  object 
of  interest  sample  different  volumes  of  optical  turbulence.  This  effect  occurs  if  there  is 
either  a  spatial  separation  between  the  object  and  the  beacon,  or  a  spatial  separation 
between  the  wave-front  sensor  and  phase-compensation  aperture,  or  if  both  types  of  sepa¬ 
ration  are  present  in  the  AO  system.  Anisoplanatism  results  in  an  increased  value  of  the 
aperture-averaged  residual  phase  variance  after  AO  compensation,  which  causes  an  expo¬ 
nential  decrease  in  system  performance.  This  dissertation  offers  a  theoretical  framework 
in  which  the  most-general  form  of  anisoplanatism  can  be  analyzed.  In  this  framework, 
anisoplanatism  is  characterized  by  the  correlation  of  Zernike  coefficients  of  object  and 
beacon  turbulence-induced  phase.  It  is  shown  through  theoretical  calculations  validated 
with  experimental  data,  that  while  most  Zernike  coefficients  are  uncorrelated  in  the  ab¬ 
sence  of  anisoplanatism,  moderate  amounts  of  anisoplanatism  result  in  an  enhanced  cross¬ 
correlation  of  the  Zernike  coefficients.  This  fact,  along  with  a  theoretical  description  of  the 
performance-degradation  process  motivates  the  application  of  optimal  estimation  theory 
to  AO  compensation.  Using  the  general  anisoplanatism  analysis  framework,  the  optimal 
compensation  theory  is  applied  to  imaging,  beam-projection,  and  interferometric  systems. 
In  imaging  systems,  optimal  compensation  reduces  the  eflfect  of  anisoplanatism,  increasing 
beacon  separation  limits  by  up  to  three  times,  and  enabling  estimation  of  wave-front  tilt 
from  a  laser  beacon,  which  cannot  be  measured  conventionally.  For  beam-projection  sys¬ 
tems  such  as  Airborne  Laser,  optimal  compensation  can  decrease  the  AO  update  rate  by 
over  70%.  In  interferometric  systems,  a  fringe  visibility  increase  of  25%  is  realized  by  ap¬ 
plying  optimal  compensation,  resulting  in  enhanced  image  reconstruction  at  large  baseline 
separations. 
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OPTIMAL  ATMOSPHERIC  COMPENSATION  FOR 


ANISOPLANATISM  IN  ADAPTIVE-OPTICAL  SYSTEMS 

I.  Introduction 

1.1  Overview 

A  fundamental  result  in  the  study  of  optics  is  that  all  optical  systems  are  limited  in 
resolution  by  the  effects  of  diffraction.  For  a  system  unaffected  by  fixed  aberrations,  the 
limiting  angular  resolution  due  to  diffraction,  designated  9^  is  often  expressed  as  [35]: 

oc  (1) 

where  A  is  the  propagation  wavelength  of  the  electromagnetic  radiation,  and  D  is  the 
aperture  diameter  of  the  optical  system.  Thus,  the  limiting  angular  resolution  of  an  optical 
system  is  decreased  by  imaging  at  shorter  wavelengths,  or  by  increasing  the  aperture  size. 
For  imaging  at  visible  wavelengths  (A  =  500  nm),  the  diffraction  limit  is  approximately 
5  //rad  or  1  arcsec  for  a  10  cm  aperture,  0.1  arcsec  for  a  1  m  aperture,  0.01  arcsec  for  a 
10  m  aperture,  and  so  on. 

The  diffraction  limit  in  Eq.  1  applies  strictly  to  an  optical  system  operating  in  free- 
space  or  vacuum,  which  is  a  reasonable  description  of  most  laboratory  systems.  When 
imaging  through  extended  atmospheric  paths,  however,  the  assumption  of  vacuum  prop¬ 
agation  is  no  longer  valid.  In  this  case,  atmospheric  effects  must  be  considered.  Even 
when  imaging  through  clear  air  so  that  atmospheric  scattering  and  absorption  effects  are 
negligible,  phase  aberration  effects  caused  by  atmospheric  turbulence  degrade  the  resolu¬ 
tion  of  an  imaging  system.  It  has  been  found,  both  theoretically  and  experimentally,  that 
the  resolution  of  an  optical  system  is  limited  to  approximately  1  arcsec  no  matter  how 
large  the  aperture  diameter  [20].  Increasing  the  size  of  a  telescope’s  aperture  will  increase 
the  irradiance  incident  on  the  system’s  focal  plane,  but  will  not  provide  higher-resolution 
imagery  for  reasonable  light  levels.  Therefore,  any  optical  system  operating  at  visible  wave- 
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UNCOMPENSATED 


ADAPTIVE  OPTICS  COMPENSATION 


Figure  1.  Example  of  adaptive  optics  compensation  in  astronomical  imaging.  These 
images  show  the  binary  star  k-Peg  as  seen  through  the  3.5  m  telescope  at 
the  Starfire  Optical  Range  without  atmospheric  compensation  (left)  and  with 
adaptive-optical  correction  applied  (right).  Actual  binary  star  separation  is 
0.3  arcsec.  (Photo  courtesy  Air  Force  Research  Laboratory) 

lengths  with  an  aperture  larger  than  approximately  10  cm  is  turbulence-limited  when  used 
for  imaging  over  extended  atmospheric  paths,  since  the  limiting  resolution  angle  imposed 
by  turbulence  is  larger  than  that  imposed  by  diffraction. 

To  improve  the  resolution  of  turbulence- limited  optical  systems,  adaptive  optics  (AO) 
compensation  technology  has  been  developed  [69].  In  an  ideal  AO  system,  the  turbulence- 
induced  phase  incident  from  the  object  being  observed  is  measured  using  a  wave-front 
sensing  device,  and  the  conjugate  phase  is  formed  on  the  surface  of  a  compensation  device. 
Reflection  from  this  surface  removes  phase  aberrations  from  the  incident  wave-front,  and 
system  performance  is  restored  to  its  diffraction  limit.  An  example  of  AO  compensation 
applied  to  astronomical  imaging  is  shown  in  Figure  1.  These  long-exposure  images  show 
a  close  binary  star  pair  as  seen  through  the  3.5  m  telescope  operated  at  the  Air  Force 
Research  Laboratory’s  Starfire  Optical  Range.  Without  AO  compensation,  the  binary 
pair  cannot  be  resolved,  since  the  separation  is  0.3  arcsec.  The  star  pair  is  imaged  as  a 
bright  spot  approximately  1  arcsec  in  diameter.  When  AO  compensation  is  applied,  each 
of  the  stars  is  clearly  resolved,  and  the  actual  separation  is  apparent.  The  image  of  each 


star  represents  a  system  point  spread  function  (PSF)  which  is  nearly  diffraction-limited 
(approximately  0.03  arcsec  for  this  system). 

The  performance  of  a  practical  adaptive-optical  compensation  system  can  be  de¬ 
graded  by  a  number  of  factors  related  to  system  design.  Inaccuracies  in  wave-front  mea¬ 
surement  due  to  finite  light  levels,  inadequate  wave-front  spatial  sampling,  or  poor  phase 
reconstruction  techniques  lead  to  reduced  AO  system  performance  [65].  Phase  conjuga¬ 
tion  errors  caused  by  deformable  mirror  physical  constraints  and  control  defects  also  limit 
the  benefits  of  using  adaptive-optics  [64].  Even  when  wave- front  sensor  and  phase  conju¬ 
gation  effects  are  negligible,  adaptive  optics  compensation  may  still  be  degraded  due  to 
the  nature  of  atmospheric  turbulence,  and  the  manner  in  which  turbulence  measurements 
are  obtained  in  an  AO  system.  This  effect  is  classified  broadly  as  turbulence-induced 
anisoplanatism  [21,22,68] 

The  term  “anisoplanatism”  as  it  applies  to  adaptive  optics  is  borrowed  from  termi¬ 
nology  associated  with  fixed  optics.  An  isoplanatic  region  A  is  any  region  in  which  the 
imaging  PSF,  designated  here  as  K,  is  invariant  with  respect  to  a  spatial  translation  [9]: 

(2) 

where  Ka  represents  the  diffraction-limited  PSF  centered  at  the  position  f.  The  term 
“anisoplanatic”  is  therefore  used  to  imply  that  the  isoplanatic  condition  of  Eq.  2  does  not 
hold. 

In  the  context  of  adaptive  optics,  anisoplanatism  typically  refers  to  a  performance 
degradation  that  results  when  the  optical  source  used  as  a  beacon  for  wave- front  sensing  is 
separated  from  the  object  being  observed.  To  understand  the  physical  origin  of  turbulence- 
induced  anisoplanatism,  consider  the  diagrams  shown  in  Figure  2(a)-(b).  For  the  ideal  AO 
system  of  Figure  2(a),  the  object  of  interest  is  bright  enough  to  be  used  as  the  AO  beacon 
for  wave-front  sensing.  Thus,  neglecting  wave-front  sensing  errors,  the  measured  phase 
aberration  in  the  system  aperture  is  the  same  as  the  object  phase  aberration,  and  ideal 
phase  conjugation  removes  all  wave-front  errors  induced  by  propagation  through  optical 
turbulence.  When  the  object  of  interest  is  too  dim  to  provide  for  wave-front  sensing,  a 
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Figure  2.  Physical  origin  of  turbulence-induced  anisoplanatism  in  adaptive-optical  sys¬ 
tems.  (a)  An  ideal  system,  where  the  object  of  interest  is  also  used  as 
the  AO  beacon,  (b)  Anisoplanatism  introduced  by  separation  between  ob¬ 
ject  and  beacon,  (c)  Anisoplanatism  introduced  by  separation  between  wave¬ 
front-sensing  aperture  and  wave-front-compensation  aperture,  (d)  Generalized 
anisoplanatism  in  an  AO  system,  with  both  object/beacon  source  and  sens¬ 
ing/compensation  aperture  separation. 
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separate  bright  source  or  beacon  must  be  used  to  obtain  the  phase  aberration  information, 
as  depicted  in  Figure  2(b).  Light  propagating  from  the  object  and  the  beacon  sample 
different  volumes  of  turbulence.  As  a  result,  the  turbulence-induced  phase  incident  in 
the  aperture  from  the  object  and  the  source  will  not  be  equal,  but  instead,  only  partially 
correlated. 

The  degree  of  correlation  between  the  object  and  beacon  phase  depends  upon  the 
difference  in  the  turbulence  volumes  sampled  by  the  light  as  it  propagates  from  each  source 
to  the  system  aperture.  According  to  conventional  theory  [44,78],  turbulence  is  stochas¬ 
tic  in  nature,  and  must  be  characterized  by  its  statistical  properties.  Stochastic  models 
for  turbulence  specify  second-order  statistics  such  as  the  spatial  correlation  of  the  index  of 
refraction  fluctuation,  which  depends  upon  the  spatial  separation  between  turbulence  sam¬ 
pling  points.  It  therefore  follows  that  the  statistical  nature  of  turbulence-induced  phase 
can  also  be  characterized  by  spatial  correlation  properties,  which  are  determined  by  the 
spatial  separation  of  fleld  points  along  the  propagation  path  from  each  source.  The  degra¬ 
dation  in  AO  compensation  performance  associated  with  the  conventional  anisoplanatism 
geometry  shown  in  Figure  2(b)  results  directly  from  the  partial  correlation  of  object  and 
beacon  phase. 

This  understanding  of  the  physical  origin  of  anisoplanatism  leads  naturally  to  an  un¬ 
derstanding  of  this  phenomenon  encompassing  a  more-diverse  class  of  AO  system  configu¬ 
rations.  Figure  2(c)  illustrates  the  turbulence-sample  volume  differences  resulting  from  a 
separation  of  the  wave-front-sensing  aperture  and  wave-front-compensation  aperture  when 
the  object  is  bright  enough  to  be  used  as  the  AO  beacon.  In  this  configuration,  the  sensed 
wave-front  will  again  be  only  partially  correlated  with  the  conjugate  phase  required  for 
proper  compensation.  The  most  general  form  of  this  decorrelation  process  occurs  when 
both  the  object  and  beacon  as  well  as  the  sensing  aperture  and  compensation  aperture  are 
separated,  as  depicted  in  Figure  2(d).  The  system  configurations  shown  in  Figure  2(b)-(d) 
are  affected  by  anisoplanatism  in  a  general  sense;  that  is,  the  object  and  beacon  phase  are 
only  partially  correlated.  This  generalized  view  of  anisoplanatism  naturally  gives  rise  to 
a  unified  approach  to  the  analysis  of  diverse  AO  system  configurations.  In  this  unified 
approach,  the  aperture  and  source  geometry  is  generalized,  and  the  phase  correlation  is 
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quantified  as  a  function  of  the  spatial  position  of  the  object,  beacon,  wave-front-sensing 
aperture  and  wave-front-compensation  aperture.  The  object  and  beacon  phase  correlation 
quantifies  the  manifestation  of  anisoplanatism  in  the  given  system  configuration. 

The  generalized  view  of  anisoplanatism  espoused  in  this  research  is  enlightening,  as 
it  exposes  the  mechanism  of  anisoplanatic  degradation  in  AO  compensation  performance. 
The  ability  to  quantify  the  spatial  correlation  of  the  turbulence-induced  phase  leads  directly 
to  the  application  of  optimal  estimation  techniques  [43]  in  AO  compensation.  When  ap¬ 
plied  to  AO  compensation  for  anisoplanatic  systems,  the  wealth  of  knowledge  represented 
by  optimal  estimation  theory  produces  immediate  dividends.  The  result  is  an  optimal 
atmospheric  compensation  analysis  technique  which  indicates  how  the  performance  of  sys¬ 
tems  affected  by  anisoplanatism  may  be  improved.  Thus,  the  limits  on  optical  performance 
imposed  by  the  nature  of  turbulence  are  expanded,  and  diffraction-limited  performance 
over  extended  turbulence  paths  is  made  more  accessible. 

1,2  Historical  background  and  research  motivation 

The  effects  of  the  atmosphere  on  optical  propagation  are  rooted  in  everyday  experi¬ 
ence  with  natural  phenomena.  Negative  temperature  gradients  resulting  from  solar  loading 
on  the  earth’s  surface  give  rise  to  “mirages”  in  the  distance,  commonly  observed  by  desert 
travelers  or  highway  drivers  looking  over  an  asphalt  surface.  These  mirages  are  actually 
images  of  the  sky  and  surrounding  landscape  produced  by  the  refractive  qualities  of  heated 
air.  At  night,  positive  temperature  gradients  produce  a  “looming”  effect,  which  enables  us 
to  see  objects  below  the  horizon.  Indeed,  we  are  all  familiar  with  the  “twinkling”  of  stars 
or  distant  city  lights  which  is  caused  by  refractive  effects  of  optical  turbulence. 

Atmospheric  effects  were  considered  in  the  earliest  of  modern  treatises  on  optics. 
Christiaan  Huygens  (1629-1695)  noted  the  steady  variations  in  the  optical  properties  of 
the  atmosphere  [41]: 

If  then  you  look  through  [a  telescope]  at  different  hours  of  the  day,  leaving  it 
always  fixed  in  the  same  way,  you  will  see  that  the  same  spots  will  not  always 
appear  at  the  middle  of  the  aperture  of  the  telescope,  but  that  generally  in  the 
morning  and  in  the  evening,  when  there  are  more  vapours  near  the  earth,  these 
objects  seem  to  rise  higher,  so  that  the  half  or  more  of  them  will  no  longer  be 
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visible;  and  so  that  they  seem  lower  toward  mid-day  when  these  vapours  are 
dissipated. 

Although  Isaac  Newton  (1642-1727)  attributed  a  different  fundamental  nature  to  light  than 
Huygens,  he  seemed  to  understand  the  limits  that  atmospheric  effects  place  on  optical 
systems  [53]: 

If  the  Theory  of  making  Telescopes  could  at  length  be  fully  brought  into  Practice, 
yet  there  would  be  certain  Bounds  beyond  which  Telescopes  could  not  perform. 

For  the  Air  through  which  we  look  upon  the  Stars,  is  in  a  perpetual  Tremor;  as 
may  be  seen  by  the  tremulous  Motion  of  Shadows  cast  from  high  Towers,  and 
by  the  twinkling  of  the  fix^d  Stars. 

After  the  time  of  Huygens  and  Newton,  research  in  optics,  electricity,  and  magnetism  de¬ 
veloped  almost  independently,  until  James  Clerk  Maxwell  (1831-1879)  was  able  to  describe 
all  known  electromagnetic  phenomena  in  terms  of  a  system  of  field  equations,  leading  to 
the  conjecture  that  light  waves  are  electromagnetic  waves.  This  insight,  and  the  experi¬ 
ments  that  followed  to  test  the  theory  has  lead  directly  to  our  current  understanding  of 
optical  phenomena,  including  optical  propagation  through  turbulent  media. 

Maxwell’s  field  theory  also  advanced  the  development  of  electronic  components,  en¬ 
abling  a  practical  design  theory  for  electronic  circuits  to  be  developed.  By  the  early  1950s, 
Babcock  [3, 4]  understood  that  an  electronically-controlled  system  could  be  constructed 
to  correct  for  phase  aberrations  introduced  by  propagation  through  turbulence,  resulting 
in  improved  optical  performance  when  the  corrected  light  was  focused  by  a  lens.  Bab¬ 
cock  also  pointed  out  that  the  reference  star  must  be  of  adequate  brightness  to  provide 
for  wave-front  sensing,  and  that  good  compensation  performance  is  limited  to  a  small 
angular  extent  around  the  reference  star.  Thus,  anisoplanatism  was  understood  to  be  a 
performance-limiting  factor  from  the  earliest  days  of  adaptive  optics  technology. 

The  science  of  atmospheric  optics  was  advanced  in  the  early  1960s  by  Tatarski  [78], 
who  drew  upon  the  turbulence-characterization  work  of  Kolmogorov  [44]  to  develop  a 
stochastic  model  for  the  effects  of  turbulence  in  optical  systems.  Tatarski’s  work  was  used 
by  Hufnagel  and  Stanley  [40]  and  Pried  [20]  as  a  basis  for  understanding  the  degradation 
of  optical  resolution  caused  by  atmospheric  turbulence.  Specifically,  Fried  introduced  the 
concept  of  an  atmospheric  coherence  diameter  Tq,  demonstrating  that  the  resolution  of 
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optical  systems  with  aperture  diameter  less  than  Tq  is  limited  by  diffraction,  whereas 
turbulence  effects  limit  the  resolution  of  systems  with  aperture  diameter  larger  than  tq* 
Fried’s  work  illustrates  theoretically  that  optical  systems  reach  a  limiting  performance 
determined  by  the  statistical  properties  of  the  turbulence,  not  by  the  physical  extent  of 
the  aperture.  Pried  also  introduced  modal  analysis  techniques  in  the  analysis  of  turbulence- 
induced  phase  [19]. 

With  continued  growth  in  a  theoretical  understanding  of  atmospheric  turbulence  [37, 
47,54,92]  and  technology  advances  in  the  late  1960s  and  1970s  [34],  adaptive  optics  became 
a  viable  technology  for  imaging  applications,  especially  in  astronomy.  With  practical  AO 
systems  being  tested,  effects  of  anisoplanatism  began  to  be  investigated  theoretically  [21, 
45].  Following  Fried’s  convention  [22],  anisoplanatism  effects  were  characterized  by  an 
isoplanatic  angle,  designated  Oq.  Pried  argued  that  optical  performance  is  degraded  by 
a  factor  exp[— where  0  represents  the  angular  separation  between  the  object 
of  interest  and  a  natural  guide  star  (NGS)  used  as  the  AO  beacon.  Thus,  the  extent 
of  “angular  anisoplanatism”  scales  according  to  the  isoplanatic  angle.  This  fact  has  also 
been  exploited  in  theoretical  works  regarding  the  temporal  analysis  of  AO  systems  [84]. 
Others  [74]  have  questioned  the  appropriateness  of  this  measure  for  angular  anisoplanatism. 
Alternative  methods  of  establishing  performance  limits  imposed  by  angular  anisoplanatism 
have  been  addressed  by  Troxel  et  al.  [82,83]. 

The  period  of  the  mid-1980s  through  the  mid-1990s  was  marked  by  significant  ad¬ 
vances  in  adaptive  optics  technology.  In  1985,  Foy  and  Labeyrie  first  introduced  the 
concept  of  laser  guide  star  (LGS)  adaptive  optics  in  the  open  literature  [18]  to  provide 
increased  sky  coverage  for  atmospheric  compensation.  The  first  experimental  results  were 
reported  by  Thompson  et  al.  in  1987  [80].  Atmospheric  composition  using  the  laser  guide 
star  was  first  reported  by  Primmerman  et  al.  [56]  and  Fugate  et  al.  [25].  The  finite  altitude 
of  the  LGS  for  adaptive  optics  introduces  a  form  of  performance  degradation  referred  to 
as  “focal  anisoplanatism.”  This  effect  has  been  investigated  in  detail  in  works  by  Gardner 
et  al  [27],  Welsh  and  Gardner  [93],  Tyler  [87],  Fried  and  Belsher  [23],  Esposito  et  al.  [17], 
and  Molodij  and  Rousset  [51]. 
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Concurrent  with  many  of  the  recent  advances  in  adaptive  optics  technology,  theoret¬ 
ical  modeling  of  compensation  performance  has  lead  to  the  application  of  optimal  control 
in  AO  systems.  The  use  of  optimal  estimation  techniques  in  atmospheric  compensation 
dates  back  to  the  work  of  Valley  and  Wandzura  [88],  who  pointed  out  that  an  optimal  com¬ 
pensation  could  be  achieved  in  a  system  whose  performance  is  degraded  by  anisoplanatism 
by  estimating  the  object  phase  given  beacon  phase  measurements  and  knowledge  of  turbu¬ 
lence  statistics.  Optimal  phase  reconstruction  from  wave-front  sensor  slope  measurements 
has  also  been  explored  in  the  spatial  domain  [16,86,91]  and  in  the  time  domain  [15,90,97]. 
Optimal  modal  compensation  for  NGS  and  LGS  adaptive  optics  has  recently  been  devel¬ 
oped  by  Whiteley  et  al.  [96]  to  partially-recover  performance  degradation  in  astronomical 
imaging  systems  degraded  by  both  angular  and  focal  anisoplanatism. 

Throughout  the  history  of  adaptive  optics  research,  the  concept  of  anisoplanatism  has 
been  expanded  due  to  changes  in  adaptive  optics  technology.  For  instance,  consideration 
of  the  focal  anisoplanatic  effect  with  a  laser  guide  star  represented  an  expansion  in  the 
original  concept  of  anisoplanatism  as  introduced  by  Babcock  and  Pried.  Performance 
degradation  was  no  longer  related  to  an  angular  field  separation,  but  rather  to  an  altitude 
difference  between  the  object  and  LGS  beacon.  Prior  to  the  advent  of  laser-guide  star 
adaptive  optics,  focal  anisoplanatism  had  not  been  addressed  in  system  analysis  efforts. 
Focal  anisoplanatism  became  relevant  only  after  LGS  beacon  technology  was  developed. 

Today,  adaptive  optics  is  being  applied  to  new  problems,  including  dynamic  beam 
projection  systems  like  the  Airborne  Laser  [26],  multiple  beacon  compensation  [42,72], 
and  optical  interferometric  systems  [75].  Ground-based  and  spaced-based  variants  of  the 
Airborne  Laser  are  being  considered  for  weapons,  propulsion,  and  laser  communication  ap¬ 
plications.  Interferometric  imaging  applications  with  adaptive  optics  are  promising  tech¬ 
nologies  for  advanced  surveillance  of  space  assets  and  astronomy.  While  conventional 
analysis  of  anisoplanatism  has  a  rich  history  dating  back  to  the  earliest  days  of  adaptive 
optics  technology,  current  analysis  methods  are  inadequate  to  treat  the  diverse  nature  of 
emerging  AO  applications.  The  concept  of  anisoplanatism  and  the  associated  analysis 
methods  must  be  expanded  to  treat  these  new  systems  so  their  performance  capabilities 
and  limitations  may  be  properly  assessed. 
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1.3  Problem  statement  and  research  scope 

Although  adaptive  optics  has  its  origin  in  astronomical  applications,  and  much  of 
the  current  research  applies  directly  to  those  applications,  AO  technology  is  now  being  in¬ 
corporated  into  advanced  weapon  and  surveillance  systems.  The  current  analysis  methods 
for  anisoplanatism  are  inadequate  to  treat  the  broadened  scope  of  adaptive  optics  tech¬ 
nologies.  The  desire  to  establish  the  performance  capabilities  and  limits  of  this  expanded 
class  of  AO  systems  is  the  key  motivation  for  the  problem  addressed  in  this  dissertation: 

Establish  a  unified  framework  for  the  analysis  of  anisoplanatism  in  a  broad  class 
of  adaptive-optical  ground-based  imaging  systems,  dynamic  beam-projection 
systems,  and  interferometric  imaging  systems;  and  quantify  performance  im¬ 
provement  obtained  through  optimal  atmospheric  compensation  techniques  ap¬ 
plied  to  each  system  type. 

To  this  end,  it  is  assumed  that  the  most-general  form  of  anisoplanatism  is  mani¬ 
fested  in  an  interferometric  system,  where  both  the  object  source  and  beacon  source  as 
well  as  the  wave-front-sensing  aperture  and  wave-front-compensation  aperture  may  be  sep¬ 
arated  spatially.  This  assumption  serves  as  the  primary  method  of  unifying  the  analysis 
of  imaging,  beam-projection,  and  interferometric  systems.  Performance  degradation  due 
to  turbulence-induced  anisoplanatism  and  the  associated  optimal  compensation  theory  are 
developed  for  the  case  of  an  interferometric  system,  with  imaging  and  beam-projection 
systems  being  a  special  case  of  this  more-general  analysis.  The  results  obtained  in  this 
research  are  for  multi-layer,  wide-sense  stationary,  isotropic  turbulence  with  finite  outer 
scale.  For  ground-based  AO  imaging,  NGS  and  LGS  compensation  methods  are  explored, 
demonstrating  the  performance  improvement  obtained  with  optimal  compensation.  Tilt 
compensation  performance  is  investigated  in  the  case  of  the  Airborne  Laser  system  both 
with  and  without  optimal  compensation.  For  an  adaptive-optical  interferometer,  the  effect 
of  angular  anisoplanatism  on  system  performance  is  investigated,  and  optimal  compensa¬ 
tion  theory  is  applied  to  provide  enhanced  visibility  measurements  with  AO  compensation. 

1.4  Key  results 

The  theory  and  analysis  conducted  in  this  research  and  presented  in  this  document 
have  lead  to  the  following  key  results: 
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•  The  adaptive  optics  compensation  performance  of  imaging,  beam-projection,  and 
interferometric  systems  may  be  analyzed  in  a  unified  context,  where  imaging  and 
beam-projection  applications  are  viewed  as  a  limiting  case  of  an  adaptive-optical  in¬ 
terferometer.  An  optimal  modal  compensation  theory  for  anisoplanatism  in  adaptive 
optics  is  developed,  where  the  generalized  effect  of  anisoplanatism  is  characterized 
by  the  cross-correlation  of  Zernike  phase  coefficients. 

•  A  theoretical  and  experimental  analysis  of  the  cross-correlation  of  Zernike  phase 
coefficients  with  finite  turbulence  outer  scale  demonstrates  that  while  most  Zernike 
coefficients  are  uncorrelated  when  no  anisoplanatism  is  present  in  an  AO  system, 
the  cross-correlation  of  many  Zernike  modes  is  enhanced  when  moderate  levels  of 
anisoplanatism  exist. 

•  In  ground-based  AO  imaging  systems,  application  of  the  optimal  compensation  the¬ 
ory  results  iti  increased  Strehl  ratio  performance  when  angular  and  focal  anisopla¬ 
natism  are  present  in  the  system.  For  NGS  compensation,  equivalent  performance 
can  be  obtained  at  larger  angular  beacon  offsets.  For  LGS  compensation,  equivalent 
performance  can  be  obtained  at  lower  LGS  altitudes.  The  optimal  compensation 
theory  leads  to  a  new  class  of  LGS  tilt-sensing  methods,  where  tilt  components  are 
estimated  from  higher-order,  off-axis  LGS  measurements. 

•  When  time-delays  are  present  in  a  dynamic  beam-projection  system  like  the  Air¬ 
borne  Laser,  application  of  optimal  modal  compensation  results  in  increased  Strehl 
ratio  performance  at  longer  time  delays.  Using  the  optimal  estimator  for  tilt  com¬ 
ponents,  higher-order  modal  measurements  can  be  used  to  enhance  beam-steering 
performance.  Consequently,  the  AO  bandwidth  for  the  Airborne  Laser  system  can 
be  reduced  by  over  70%,  while  maintaining  a  constant  performance  level, 

•  The  theoretical  generalized  anisoplanatism  framework  permits  the  analysis  of  angu¬ 
lar  anisoplanatism  in  adaptive-optics  for  interferometry.  Angular  anisoplanatism  is 
shown  to  impose  more  severe  limitations  in  an  adaptive-optical  interferometer  than  in 
an  imaging  AO  system.  Application  of  the  optimal  compensation  theory  in  adaptive- 
optical  interferometry  leads  to  a  simplified  system  design,  as  the  optimal  conjugate 
phase  can  be  formed  on  a  single  deformable  mirror,  based  on  wave-front  measure- 
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ments  in  each  baseline  aperture.  Optimal  compensation  also  results  in  increased 
fringe  visibility  for  larger  AO  beacon  angular  offsets,  leading  to  a  greater  capability 
to  reconstruct  object  irradiance  information  from  interferometric  measurements. 

1,5  Dissertation  organization 

The  remainder  of  this  dissertation  is  organized  into  six  chapters.  Chapter  II  presents 
the  background  and  analysis  required  to  understand  the  effects  of  turbulence  and  aniso- 
planatism  on  the  performance  of  optical  systems.  Turbulence  statistics  and  modeling  are 
discussed,  and  wave  propagation  in  turbulent  media  is  explored.  On  the  basis  of  this 
standard  theory,  optical  system  performance  is  analyzed  in  imaging,  beam-projection,  and 
interferometric  applications.  Adaptive  optics  technology  and  modal  analysis  of  compensa¬ 
tion  performance  are  introduced,  leading  to  a  quantification  of  the  performance  degradation 
that  results  from  both  angular  and  focal  anisoplanatism. 

Optimal  compensation  for  anisoplanatic  systems  is  discussed  in  Chapter  III.  In  this 
analysis,  the  compensation  performance  of  an  adaptive-optical  interferometer  is  consid¬ 
ered  using  modal  analysis  techniques.  Compensation  is  shown  to  be  suboptimal  when 
beacon  measurements  are  used  directly.  Optimal  compensation  requires  that  the  con¬ 
jugate  phase  be  estimated  from  modal  beacon  measurements,  given  knowledge  of  the 
spatial  cross-correlation  of  the  Zernike  modes.  The  form  of  the  optimal  estimator  is 
specified  and  its  performance  is  quantified,  indicating  the  need  to  evaluate  the  Zernike 
coefficient  cross-correlation  for  a  generalized  analysis  geometry.  The  cross-correlation  cal¬ 
culation  is  detailed,  and  the  resulting  expression  is  validated  by  comparison  with  theoretical 
and  experimental  results.  The  results  of  Chapter  III  are  also  applicable  to  imaging  and 
beam-projection  systems,  by  considering  limiting  cases  of  the  more-general  results  for  the 
adaptive-optical  interferometer. 

Chapter  IV,  Chapter  V,  and  Chapter  VI  deal  with  the  application  of  the  optimal  com¬ 
pensation  theory  to  adaptive  optics  for  ground-based  imaging  systems,  a  beam-projection 
system  (Airborne  Laser),  and  an  optical  interferometer,  respectively.  For  each  application, 
the  general  theory  is  specialized,  refining  notation,  and  considering  the  peculiar  features  of 
each  system  type.  Optimal  compensation  is  shown  to  yield  enhanced  performance  and  to 
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expand  the  limitations  placed  on  system  operation  by  turbulence-induced  anisoplanatism. 
Conclusions  from  this  research  and  recommendations  for  future  research  are  summarized 
in  Chapter  VIL 

1,6  Summary 

In  this  chapter,  atmospheric  compensation  through  adaptive  optics  has  been  briefly 
introduced,  and  the  concept  of  anisoplanatism  was  discussed  as  a  factor  which  degrades 
compensation  performance.  Since  AO  technology  has  its  roots  in  instrumentation  for 
astronomical  imaging,  the  analysis  of  anisoplanatism  has  conventionally  focused  on  the 
system  geometries  experienced  in  those  applications.  Directed-energy  weapon  systems 
like  the  Airborne  Laser  and  interferometric  systems  for  high-resolution  surveillance  also 
experience  performance  degradation  with  the  same  physical  origin  as  conventional  aniso¬ 
planatism,  but  these  effects  are  not  easily  explained  by  conventional  analysis  methods. 
In  order  to  analyze  these  systems,  the  concept  of  anisoplanatism  must  be  generalized  to 
include  any  physical  conflguration  which  results  in  a  partial  correlation  of  the  sensed  wave- 
front  aberration  and  the  required  conjugate  phase.  To  achieve  this  goal,  a  unified  approach 
to  understanding  turbulence  degradation  and  anisoplanatic  effects  will  be  taken  for  a  broad 
class  of  optical  systems  in  Chapter  II.  This  unified  analysis  leads  to  the  development  of  an 
optimal  compensation  theory  for  adaptive  optics  in  Chapter  III  which  requires  the  gener¬ 
alized  form  of  anisoplanatism  be  quantified  for  a  modal  phase  decomposition.  Application 
of  the  optical  compensation  theory  for  anisoplanatism  in  adaptive  optics  will  be  discussed 
in  Chapters  IV,  V,  and  VI  for  ground-based  imaging,  beam-projection  (Airborne  Laser), 
and  interferometric  imaging  systems. 
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IL  Atmospheric  and  Adaptive  Optics 

2.1  Introduction 

The  purpose  of  this  chapter  is  to  review  the  fundamental  concepts  associated  with 
atmospheric  and  adaptive  optics,  and  to  provide  a  theoretical  framework  in  which  the 
effects  of  turbulence  on  conventional  and  interferometric  optical  systems  can  be  under¬ 
stood.  Modal  compensation  of  atmospheric  aberrations  is  discussed,  and  the  concept 
of  turbulence-induced  anisoplanatism  in  adaptive  optics  is  introduced.  Notation  used 
throughout  the  remaining  chapters  is  established  here,  and  quantities  important  to  later 
analyses  are  derived. 

Section  2.2  of  this  chapter  explores  the  effects  of  turbulence-induced  phase  aberra¬ 
tions  on  optical  systems.  The  section  begins  with  a  discussion  of  the  statistics  of  index 
of  refraction  fluctuations  in  the  atmosphere,  leading  to  the  subject  of  wave  propagation 
in  turbulent  media  and  turbulence  modeling.  A  theoretical  development  is  then  presented 
which  can  be  used  to  quantify  turbulence  effects  on  conventional  imaging,  beam-projection, 
and  interferometric  systems. 

Adaptive  optics  technology  is  overviewed  in  Section  2.3.  A  typical  adaptive  op¬ 
tics  system  is  considered,  and  the  function  of  each  system  component  is  described.  The 
unique  characteristics  of  natural  and  laser-guide-star  adaptive  optics  are  discussed,  and 
the  concept  of  modal  atmospheric  compensation  is  introduced  as  a  means  of  quantifying 
the  performance  of  an  adaptive  optics  system. 

Section  2.4  discusses  turbulence-induced  anisoplanatism  in  adaptive  optics.  Aniso¬ 
planatism  is  studied  using  modal  analysis  techniques,  and  its  effects  on  adaptive-optical 
system  performance  is  explored.  This  development  sets  the  stage  for  the  analysis  presented 
in  Chapter  III,  when  optimal  methods  for  modal  compensation  of  atmospheric  aberrations 
are  considered.  A  brief  summary  of  the  chapter  is  presented  in  Section  2.5. 

2.2  Turbulence  effects  on  optical  systems 

The  effects  of  turbulence  in  optical  systems  arise  from  variations  in  the  index  of 
refraction  along  the  propagation  path.  Large-scale  variations  of  n  at  the  spatial  posi- 
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tion  ^  can  be  related  to  slowly-varying  atmospheric  characteristics  such  as  pressure  and 
temperature  by  the  following  steady-state  deterministic  relationship  [71]: 


n 


(f)  -  1  =  77.6  X  10-'  ) 


pji) 

t{() 


(3) 


where  A  is  the  propagation  wavelength  measured  in  /xm,  -P  is  the  atmospheric  pressure 
in  millibars,  and  ^  is  the  atmospheric  temperature  in  degrees  Kelvin.  Smaller-scale 
variations  in  the  index  of  refraction  must  be  regarded  as  stochastic,  depending  upon  random 
changes  in  local  atmospheric  parameters,  most  notably  temperature  and  velocity.  Thus, 
the  analysis  of  turbulence  effects  on  optical  systems  is  inherently  stochastic  in  nature, 
requiring  a  statistical  characterization  of  the  index  of  refraction  fluctuations. 

2. 2 A  Statistics  of  atmospheric  index  of  refraction  fluctuations.  Our  understand¬ 
ing  of  the  phenomenology  of  atmospheric  optical  turbulence  has  its  origin  in  the  work 
of  Kolmogorov  [44]  and  Obukov.  Kolmogorov  argued  that  when  eddies  are  formed  in  a 
turbulent  flow,  those  of  larger  scale  size  break  up  into  smaller  eddies,  and  that  the  eddy 
scale  distribution  follows  a  universal  law  based  upon  energy  conservation.  In  the  standard 
theory,  statistics  of  the  air  speed  field  are  characterized  by  the  structure  function 

Dy  sj,  which  is  defined  as: 


Dv  ({,  ^  (^ )  -  ■«  (C  +  s)]  ^  > 


(4) 


where  f  and  ^  +  s  are  spatial  locations  within  the  field,  and  ( • )  represents  averaging 
with  respect  to  an  ensemble  of  turbulence  realizations.  In  their  development,  Kolmogorov 
and  Obukov  made  two  key  assumptions.  First,  they  assumed  that  the  air  speed  field  in 
the  turbulent  flow  is  wide-sense  stationary  [32]  or  homogeneous,  such  that  the  structure 
function  depends  only  upon  the  separation  of  points  at  which  it  is  sampled,  i.e.,  Dy  5^  = 
Dy(^.  Furthermore,  they  assumed  that  the  field  was  also  isotropic;  that  is  Dy  depends 
only  upon  s  =  |^.  For  the  range  of  separations  Li  <  s  <  Lq  known  as  the  inertial 
subrange,  where  Li  is  the  turbulence  inner  scale  and  Lq  is  the  turbulence  outer  scale,  the 
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Kolmogorov-Obukov  law  says  that: 


Dy{s)  oc 


Index  of  refraction  fluctuations  n'  associated  with  a  turbulent  flow  are  analyzed 
in  a  manner  similar  to  the  air  speed  field.  The  structure  function  of  the  index  of  refraction 
fluctuations  is  designated  and  is  defined  by: 


Dn'  ^  [n'  +  ^]  Y 

By  employing  the  concept  of  conservative  passive  additives,  Tatarski  [78, 79]  showed  that 
for  homogeneous  and  isotropic  turbulence,  varies  with  separation  in  the  same  manner 
as  in  the  inertial  subrange.  Specifically,  Dn*{s^z)  is  given  by: 


D„.{s,z)  =  Cl{z)s^l^, 


(7) 


where  C^{z)  is  the  index  of  refraction  structure  constant,  which  specifies  the  steady  vari¬ 
ation  of  the  turbulence  strength  with  altitude  z.  The  units  of  C^(z)  are  meters“^/^. 


For  the  purpose  of  analysis,  it  is  often  useful  to  express  turbulence  statistics  in  the 
spatial- frequency  domain  instead  of  the  spatial  domain.  The  spectral  equivalent  of  Dn'  (s,  z) 
is  the  power  spectral  density  (PSD)  of  the  index  of  refraction  fluctuations,  designated 
^)-  If  ^  represents  a  spatial  frequency  (in  units  of  inverse  meters),  then  VFn/(A:,  z)  is 
related  to  D„/(5,z)  as  follows  [71]: 


Wn'{k,z)  =  {V^k)-^ 


rOO 

/  ds 

'sin(27rfcs)‘ 

d 
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ds 

ds 

Using  the  form  of  Dn'is^z)  given  in  Eq.  (7),  the  integral  over  $  in  Eq.  (8)  may  evaluated 
to  yield  the  following  expression  for  Wn'{k^z): 


Wn^{k,z) 


\/3r(8/3) 

2(27r)8/3 


Cl{z)k- 


-11/3 


(9) 
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The  form  of  Wni{k,z)  specified  in  Eq.  (9)  is  commonly  referred  to  as  the  Kolmogorov 
spectrum,  since  the  corresponding  structure  function  scales  with  distance  according  to 
Kolomogorov’s  law. 

y 

While  the  mathematical  form  of  Eq.  (9)  is  simple  (a  power  law),  and  has  been  used 
widely  in  the  analysis  of  atmospheric  optics,  its  application  is  limited.  When  used  in 
propagation  calculations,  Eq.  (9)  results  in  non-physical  theoretical  predictions,  such  as 
an  infinite  variance  for  the  turbulence-induced  phase,  and  a  non-zero  phase  correlation 
at  large  spatial  separations.  Note  that  as  it  is  defined  in  Eq.  (9),  Wn'{k^z)  — ^  oo  as 
fc  — ^  0.  To  provide  a  finite  limit  on  Wnf{k^z)  as  ^  0,  the  effect  of  the  turbulence  outer 
scale  is  incorporated  into  W’„/(A;,  z),  A  common  way  to  include  the  outer  scale  effect  is  as 
follows  [89] : 

=  +  «-“/•,  (10) 

where  ko  =  Lq  This  form  of  Wn'{k,z)  is  referred  to  as  the  von  Karman  spectrum.  The 
functional  dependence  of  Wn>  on  k  is  shown  in  Figure  3.  Throughout  this  research,  the  von 
Karman  spectrum  will  be  employed  as  an  appropriate  statistical  model  for  atmospheric 
optical  turbulence.  Use  of  the  von  Karman  spectrum  places  this  research  in  concert  with 
the  majority  of  other  research  efforts  employing  finite  outer  scale  models  [68,71,76,98]. 

2,2,2  Wave  propagation  in  turbulent  media.  Under  reasonable  physical  assump¬ 
tions,  Maxwell’s  equations  in  the  presence  of  a  dielectric  media  (such  as  the  atmosphere) 
imply  a  scalar  wave  equation  for  each  component  of  the  vector  electromagnetic  field  given 
by  [71]: 

V2u(|')  +  (y)'n2(|')  u(e)  =  0,  (11) 

where  the  notation  u  is  used  to  indicate  that  the  scalar  field  is  a  complex  quantity.  To  solve 
this  wave  equation,  the  Rytov  method  [78]  may  be  employed,  where  U  is  expressed  as: 

u(f)  =exp  ['J' (f)]  ,  (12) 
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Figure  3.  Effect  of  outer  scale  Lq  on  the  von  Karman  spectrum.  For  Lq  oo,  the 
Kolmogorov  spectrum  is  recovered  (dashed  line). 


and  Eq.  (11)  is  transformed  into  the  nonlinear  Riccati  equation  [71]: 

v^®(f)  +  [v»(e)]’  =  -(|)^*(«).  (13) 

Now,  if  the  index  of  refraction  fluctuations  are  assumed  to  be  small,  then  we  may  say  that 
”  (0 "  ^  en'  where  e  is  some  small  number.  A  solution  to  Eq.  (13)  is  obtained  by 
a  perturbation  method  where  it  is  assumed  that  ^  can  be  written  as  a  power  series  in 
e: 

^  (e  )  =  ^^0  (e  )  +  (i)  +  (^)  +  •  •  •  •  (14) 

Retaining  only  the  first  two  terms  results  in  the  second  order  solution  for  the  scalar  field: 


u  (f )  =  u  )  exp  [’J'l  )]  , 


(15) 
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where  represents  the  solution  to  the  wave  equation  in  free  space  [Eq.  (11)  with 

^  =1];  that  is,  the  wave  that  is  undisturbed  by  turbulence. 

The  perturbed  field  is  obtained  with  the  Green’s  function  for  free-space  prop¬ 

agation  as  [71]: 


27ru 


»(£) 


exp  ( 

j  X 

£-£'1) 

(16) 


where  the  integration  is  over  the  source  volume  denoted  by  the  primed  coordinate.  The  real 
and  imaginary  components  of  are  identified  as  x  ?  the  log-amplitude  fiuctuation 

and  the  turbulence-induced  phase,  respectively,  and  may  be  expressed  as: 


Thus,  X  (^)  computed  as  follows: 


(17) 


Hi)  =  («)-*:(«) 


(18) 

(19) 


In  this  research,  weak  turbulence  is  assumed,  and  therefore  log-amplitude  fiuctuations  are 
assumed  to  be  negligible  (x  (^)  =0).  Thus,  according  to  Eq.  (15),  the  turbulence-induced 
phase-aberrated  field  is  written  as: 


U  (f)  =  u  (f)  exp  )  . 


(20) 


If  diffraction  effects  are  also  assumed  to  be  negligible,  then  the  turbulence-induced  phase 
at  position  f  transverse  to  the  propagation  direction  represents  the  integrated  effect  of 
index  of  refraction  fluctuations  along  the  path: 


<^(f)  =  ^  dz  (y )  n'(f,  z), 


(21) 
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where  n\r^z)  is  the  index  of  refraction  fluctuation  at  a  point  specified  by  the  transverse 
vector  r  at  position  z  along  the  propagation  path  of  length  L. 

Prom  Eq.  (21),  the  statistical  properties  of  the  turbulence-induced  phase  may  be 
derived.  First,  since  can  be  viewed  as  the  sum  of  a  large  number  of  random  variables, 
all  of  which  are  zero-mean,  then  by  the  central  limit  theorem  [36],  will  be  zero- mean 
Gaussian.  Furthermore,  by  employing  the  Fourier-transform  relation  between  the  PSD 
and  the  autocorrelation,  the  PSD  of  the  turbulence-induced  phase  W^{k)  associated  with 
plane-wave  propagation  may  be  computed  as: 


w^{k)  =  TTmmr+m 

L  (t)  Jo  (t)  +  ^)} 

=  £dz(^^'^"wAKz),  (22) 


where  ^T{*}  denotes  the  two-dimensional  Fourier  transform  operation,  and  the  index 
of  refraction  fluctuations  have  been  assumed  to  be  uncorrelated  in  altitude.  Using  the 
Kolmogorov  PSD  given  in  Eq.  (9),  W^{k)  is  computed  as: 


w^k)  = 


V3r(8/3)  [/27ry 

2(27r)8/3  [xj  Jo 


dz  Cliz) 


(23) 


For  isotropic  phase,  the  phase  structure  function  is  computed  directly  from  W^{k) 
as  [71]: 

D^{s)  =  47r  r  dk  k  W^ik)[l  -  Jo{2Trks)],  (24) 

Jo 

where  Jq  is  the  Bessel  function  of  the  first  kind  of  order  zero.  Substituting  Eq.  (23)  into 
Eq.  (24),  and  performing  the  integration  with  respect  to  k  yields  the  following: 


D^{s)  =  2.91 


dzCliz) 


.5/3 


(25) 


Thus,  if  the  atmospheric  coherence  diameter  Tq  is  defined  as: 


^0  = 


(26) 
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following  the  convention  of  Pried  [20],  then  D^{s)  is  expressed  as: 

D^{s)  =  6.88  (^-j  .  (27) 

This  expression  for  the  phase  structure  function  has  been  widely  used  to  describe  the  effects 
of  turbulence- induced  phase  on  the  performance  of  optical  systems  [68].  While  Eq.  (27) 
strictly  applies  only  to  turbulence  with  infinite  outer  scale,  the  atmospheric  coherence 
diameter  ro  defined  in  Eq.  (26)  is  useful  in  characterizing  the  optical  properties  of  the 
atmosphere  even  when  finite  turbulence  outer  scale  effects  are  considered.  Since  ro  depends 
explicitly  upon  the  value  of  integrated  over  the  propagation  path,  it  can  be  used  to 
parameterize  the  turbulence  profile  in  propagation  calculations. 

2,2.3  Turbulence  modeling.  To  simplify  calculations  associated  with  propagation 
through  turbulent  media,  a  layered  atmospheric  model  is  often  assumed  [61,68,81].  In 
a  layered  atmospheric  model,  is  regarded  as  being  approximately  constant  within  the 
layer,  and  index  of  refraction  fluctuations  are  assumed  to  be  uncorrelated  from  layer  to 
layer.  For  a  layered  atmospheric  model,  the  expression  for  the  turbulence-induced  phase 
equivalent  to  Eq.  (21)  is  given  by: 

27r 

=  (28) 
^  i 

where  each  layer  is  indexed  by  the  subscript  /,  is  located  at  an  altitude  Z/,  and  is  of  thickness 

J 

Azi,  With  a  layered  model,  W^{k)  is  computed  similarly  to  Eq.  (22): 

w^{k)  =  :rT{(<^(f)^(f+s))} 

=  ^  Az,  (y)  ^  Az,-  (y)  TT  {(n'(f,z,)n'[(f,z,0  +  s])} 

=  (29) 
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The  quantity  within  the  summation  is  regarded  as  the  PSD  of  the  phase  in  the  Zth  layer, 
designated  zi).  Thus,  for  a  layered  turbulence  model,  it  may  be  said  that: 


W^{k)  =  J2W4k,z,).  (30) 

I 

For  finite  turbulence  outer  scale,  Eq.  (10)  is  substituted  into  Eq.  (29)  to  yield  the  von 
Karman  spectrum  for  phase  in  the  Ith  layer: 


Troxel  et  al.  [82]  have  shown  that  a  four-layer  turbulence  model  is  adequate  for 
many  imaging  and  propagation  calculations.  For  the  four-layer  model,  the  layer  altitude 
and  thickness  are  chosen  such  that  the  moments  of  for  continuous  turbulence  match 
the  moments  of  for  a  layered  model  through  seventh  order  (0  <  m  <  7): 


dz  z'^CKz)  =  ^  Azi  z]^  Cl{zi). 


l=l 


(32) 


For  m  =  0,  Eq.  (32)  implies  that: 


U!^dzCl{z) 


(33) 


The  term  inside  the  summation  of  Eq.  (33)  is  identified  as  a  turbulence-dependent  layer 
weight  Wi\ 

Az,  C‘(z,) 

‘  It  CJ(z) 

Using  Eq.  (34)  along  with  the  definition  of  Tq  in  Eq.  (26),  it  may  be  said  that: 


C'^(^/)  = 


6.88 

2.91 


(35) 


Thus,  substituting  Eq.  (35)  into  Eq.  (31)  yields  the  following  parametric  form  for  W^{k,  Zi): 


W4k,zi) 


v/3r(8/3)  6.88 

2(27r)8/3  ^ 


(36) 
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model 

Wi 

W2 

Ws 

W4 

ro  (m) 

HV-21 

0.8902 

0.0443 

0.0591 

0.0064 

0.050 

HV-54 

0.6877 

0.0204 

0.2554 

0.0365 

0.042 

SLC-N 

0.4965 

0.4623 

0.0299 

0.0113 

0.103 

SLC-D 

0.7397 

0.2513 

0.0048 

0.0042 

0.051 

GRNWD 

0.8615 

0.0980 

0.0394 

0.0011 

0.077 

Table  1.  Layer  weights  and  associated  atmospheric  coherence  diameter  ro  for  4-layer  at¬ 
mospheric  turbulence  models.  For  each  model,  layers  are  located  at  altitude 
Zi  —  200  m,  ^2  =  2  km,  =  10  km,  and  z^  =  18  km. 

Troxel  et  al.  [82]  calculated  the  weights  corresponding  to  a  four-layer  model  with  Zi  = 
200  m,  2^2  =  2  km,  Zs  =  10  km,  and  2:4  =  18  km.  The  layer  weights  and  corresponding 
ro  for  various  models  have  been  calculated  [82],  and  are  presented  in  Table  1.  The 
models  designated  HV-21  and  HV-54  are  the  Hufnagel- Valley  model  [39]  with  21  and  54 
mph  upper  atmospheric  winds,  respectively.  The  models  designated  SLC-N  and  SLC-D 
are  based  on  data  collected  at  the  Air  Force  Maui  Optical  Site  (AMOS),  and  modified  to 
simulate  nighttime  (N)  and  daytime  (D)  conditions  [68].  The  model  designated  GRNWD 
is  Greenwood’s  “good-seeing”  model  [33].  As  is  indicated  in  Table  1,  rp  ranges  from  4  cm 
to  over  10  cm  depending  upon  which  model  is  used  for 

2. 2,4  Temporal  turbulence  statistics.  The  discussion  of  the  statistics  of  turbulence- 
induced  phase  up  to  this  point  has  dealt  specifically  with  the  spatial  evolution  of  turbulence, 
assuming  that  any  temporal  evolution  occurs  on  time  scales  that  are  large  compared  to  the 
period  of  the  phase  measurement.  This  analysis  may  be  easily  extended  to  account  for  the 
first-order  temporal  effects  introduced  by  turbulence  motion.  To  accomplish  this  result, 
Taylor’s  frozen-flow  hypothesis  is  assumed  [32,68].  This  hypothesis  states  that  over  short 
intervals  of  time  r,  the  temporal  evolution  of  the  index  of  refraction  fluctuation  n'[(f ,  Zi),t] 
is  accurately  described  as  a  uniform  bulk  motion  due  to  atmospheric  winds  with  velocity 
vi  for  the  Zth  atmospheric  layer.  Mathematically,  this  hypothesis  is  expressed  as: 

n'[{r,z,),t  +  T]  =  n'[(f,z,)-ViT,t].  (37) 
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Thus,  temporal  delays  correspond  to  an  equivalent  spatial  separation  vit.  The  temporal 
autocorrelation  of  the  turbulence-induced  phase,  designated  is  therefore  calculated 

as  follows: 

< + 't)) 

=  (y)  +  r]) 

=  ^  Az,  (y) 

=  ^Az,  ^y^  Bn'{-ViT),  (38) 

where  Bn'  is  the  spatial  autocorrelation  of  the  index  of  refraction  fluctuations,  and  the 
turbulence  has  been  assumed  to  be  uncorrelated  layer  to  layer.  Thus,  using  the  frozen- 
flow  hypothesis,  the  temporal  autocorrelation  of  the  turbulence-induced  phase  is  related 
directly  to  the  spatial  autocorrelation  of  the  index  of  refraction  fluctuations.  Therefore, 
temporal  turbulence  statistics  may  be  understood  in  the  same  context  as  spatial  turbulence 
statistics  with  proper  accounting  for  the  time-equivalent  spatial  separation  ViT, 

2,2,5  Optical  system  performance.  The  turbulence-induced  phase  acquired  upon 
propagation  from  an  optical  source  to  the  aperture  of  an  optical  system  affects  the  perfor¬ 
mance  of  that  system.  For  a  given  realization  of  the  turbulence,  the  atmospheric  phase  will 
act  as  a  fixed  aberration,  thereby  degrading  the  performance  of  the  optical  system  from  its 
diffraction  limit  [9].  Often,  the  irradiance  from  a  low- light  level  source  on  the  focal  plane 
of  an  optical  system  must  be  integrated  temporally  to  the  enhance  signal-to-noise  ratio.  If 
this  is  the  case,  then  an  ensemble  of  turbulence  realizations  are  averaged  by  the  system, 
resulting  in  a  further  degradation  of  system  performance  [32]. 

In  order  to  quantify  the  performance  degradation  associated  with  averaging  over  an 
ensemble  of  atmospheric  phase  realizations,  consider  the  representative  multi-aperture  op¬ 
tical  system  shown  in  Figure  4.  As  shown  here,  this  diagram  is  a  schematic  representation 
of  a  spatial  optical  interferometer,  which  would  be  used  for  interferometric  imaging  ap¬ 
plications.  This  system  collects  light  with  two  baseline  apertures  of  radius  R  and  vector 
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Figure  4.  Representative  multi-aperture  system  used  for  analyzing  the  effect  of 
turbulence-induced  phase  aberrations  on  optical  system  performance.  As 
shown  here,  the  entire  system  is  an  optical  spatial  interferometer,  as  would 
be  used  for  interferometric  imaging  applications.  If  aperture  2  is  removed,  the 
remaining  system  (within  the  dashed  line)  is  a  conventional  optical  system,  as 
would  be  used  for  imaging  or  beam-projection  applications. 

separation  s,  collimates  the  light  through  lenses  of  radius  iij?,  and  combines  the  collimated 
light  using  suitable  optics  with  a  variable  path  difference  h.  The  combined  irradiance  is 
then  focused  onto  a  detector  located  on  the  focal  plane  of  a  lens  also  of  radius  Rp,  The 
analysis  which  follows  is  conducted  for  this  general  system.  It  is  recognized,  however,  that 
the  results  of  the  general  analysis  will  also  be  applicable  to  a  conventional  optical  system, 
if  only  those  quantities  associated  with  aperture  1  are  considered.  As  indicated  in  Figure  4, 
the  system  contained  within  the  dashed  line  represents  a  conventional  optical  system,  such 
as  would  be  used  for  imaging  or  beam-projection  applications. 
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The  quantity  of  interest  to  assess  optical  system  performance  is  the  long-exposure 
irradiance  on  the  focal  plane,  designated  and  calculated  as: 


{ipix))  =  (Uf  (f)u;.(f)) , 


(39) 


where  u^(x)  is  the  instantaneous  field  at  a  point  on  the  focal  plane  located  by  x.  From 
the  Fresnel  diffraction  formula,  np{x)  is  proportional  to  the  Fourier  transform  of  the  field 
incident  on  the  lens  [31]: 

exp  fjfrlxpl  r 

Upix)  = - - i  J  df}\J{fj)  P{rf)  exp 

where  /  is  the  focal  length  of  the  lens,  U {if)  is  the  field  incident  on  the  lens  at  a  point 
located  by  ?/,  P(^)  is  the  lens  pupil  function  defined  by: 


.27r  ^  ^ 


(40) 


Piv)  =  < 


1 

0 


if  1^1  <  Rf 
otherwise 


(41) 


and  the  integration  is  taken  to  be  over  all  space.  Thus,  (tir(x))  is  calculated  as: 

(iip(^))  =  ^  /  drf  I  drf'{l]m*{ff'))P{ff)P*{rf')eKp  -  rf')  ■  .  (42) 

Due  to  the  construction  of  the  interferometer,  U(^)  is  the  superposition  of  fields  from 
aperture  1  and  aperture  2,  Ui(7f)  and  U2(i^,  respectively: 


U(^  =  Ui(^)  +  U2(i^  exp  » 


(43) 


where  the  exponential  term  accounts  for  the  optical  phase  resulting  from  the  path  difference 
h.  Suppose  there  is  a  turbulence-induced  phase  associated  with  each  aperture,  designated 
(j)oi  for  propagation  from  a  source  object  o  to  aperture  1,  and  </>o25  for  propagation  from 
source  object  o  to  aperture  2.  In  this  case,  according  to  Eq.  (20),  the  expression  for  V{if) 
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in  Eq.  (43)  may  be  expanded  as; 


u(^)  =  Uoi  in)  exp[j>oi  {if)]  +  exp[j^<,2(^]  exp  y  >  (44) 

where  Voiiv)  and  U^2(^)  are  the  fields  associated  with  free-space  propagation  from  the 
source  object  to  aperture  1  and  aperture  2,  respectively.  Substituting  Eq.  (44)  into  Eq.  (42) 
and  assuming  that  the  unaberrated  field  is  independent  of  the  turbulence-induced  phase 
results  in  the  following  expression  for  • 

(if(x))  =  j  dif  j  dif'P{rf)P*{if')exp 

(U<,i(^)Uoi(^'))exp{-|  {[(l>oi{rf)  -  4>oi{rf')?)} 

+  (Uo2(^Ul2(^'))  exp  {-|  {[<f>o2{i])  -  <f>o2iii')?)} 

^  / .  . .  (  (45) 

+  (Uoi(^Uo2(^'))  exp{-|  {[(poiiv)  -  <f>o2{v')?)}  exp  {-j^h) 

,  +  (Uo2(^U*i(^'))  exp  {-4  {[<f>o2iif)  -  <f>oi{v')V)]  exp  {j^h) , 

where  the  zero-mean  Gaussian  nature  of  the  turbulence-induced  phase  has  been  used  to 
simplify  the  expected  value  of  the  exponential  terms  [20];  that  is,  for  example: 

(exp{i[^!><,i(77)  -  <t>oi{'n')]})  =  exp  |-i  {[(poiiv)  -  •  (46) 

The  expression  in  Eq.  (45)  can  be  simplified  if  the  aperture-specific  field  and  phase  quanti¬ 
ties  are  expressed  as  field  and  phase  quantities  in  the  plane  of  the  baseline  apertures,  with 
the  vector  arguments  scaled  by  the  ratio  R/Rp  to  account  for  the  effect  of  the  collimating 
lens,  and  referencing  all  vector  arguments  to  the  center  of  aperture  1: 

^o[-§^v)=Voi{v)  and  4>o{j^rf^=4>oi{rf),  (47) 

Uo  s)  =  Uo2(^)  and  (^^+s)  =  <^<,2(1?)-  (48) 

Employing  the  definition  of  the  mutual  intensity  [32]  of  the  unaberrated  fields,  designated 
J^,  and  definition  of  the  structure  function  of  turbulence- induced  phase  for  wide-sense 
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stationary  statistics: 


(49) 

(50) 

then  Eq.  (45)  can  be  rewritten  as: 


J<,(n-f2)  =  {v,{n)VM)), 
D^Xri-r2)  =  {[Mn)  -  , 


{ipix))  =  J  j  d^'^(^-P*(’?')exp  |-jy(j?-^')  •*] 

[  2J,[;^(^-7f')]exp{-|i)^,  [:^(^-^')]} 

-^iv-  v')-s]  exp  {-j^h) 
-^iv-'n')  +  s}exp{j^h). 


+^0 


+Jn 


^iv-v')-s 
-i;(v-v')  +  S 


exp{-|l>^„ 

exp{-§I>^„ 


(51) 


At  this  point,  using  the  variable  transformation  Arf  =  rf—rf'  and  =  (i7  +  ^')/2,  Eq.  (51) 
is  rewritten  as: 

(ipix))  =  j dArfPiArf)  exp  -j^Aff- X 


'25.[ 

exp 

{-kD*.  [^A?]} 

X  ^ 

+Jo 

exp{-|D,^„ 

}exp  i-j^h) 

(52) 

exp{-|I>0„ 

}exp(j^h)  , 

where  7^(A7f)  is  the  autocorrelation  of  the  pupil  function,  defined  as: 

V{Afi)  =  j  dEff  PiEff+  Aff/2)P*{Erf-  A^/2).  (53) 

By  the  definition  of  for  wide-sense  stationary  fields  [Eq.  (49)],  the  following  conjugate 
symmetry  relationship  holds: 


Jo 


R  A-  J 


j. 


^  A  -  - 

--^Atj  +  s 


(54) 


28 


Assuming  isotropic  phase  statistics,  the  phase  structure  function  is  independent  of  the 
direction  of  field-point  separation: 


si  =  +  •  (55) 

Thus,  recognizing  that  the  limits  on  Af)  are  symmetric  with  respect  to  Aff  =  0,  the 
expression  for  {ipix))  in  Eq.  (52)  can  be  rewritten  as: 


(ipix)) 


j  dAff  ViAif) 


.27r 

-j—At]  •  X 


I  Jo  exp{-|I>^„ 

\  +  I  Jo  [^A^  +  s]  I  cos  {^h  -I-  V>  0  exp  [^A^  +  s] }  , 


(56) 


where  '  is  defined  as: 

'  =  arg  I +  s  I  •  (57) 

In  the  considerations  which  follow,  Eq.  (56)  will  serve  as  a  basis  for  quantifying  turbulence 
degradation  in  optical  systems. 


2.2.5. 1  Imaging  systems.  The  performance  of  a  conventional  imaging  sys¬ 
tem  is  obtained  from  Eq.  (56),  given  three  important  considerations.  First,  an  imaging 
system  is  unaffected  by  any  phase  aberration  which  is  uniform  over  the  extent  of  the  aper¬ 
ture,  since  the  only  effect  of  this  phase  is  to  introduce  a  time  delay  in  the  image  formation 
process.  Thus,  the  turbulence-induced  phase  relevant  to  imaging  analysis  is  obtained  by 
subtracting  a  uniform  phase  shift  from  the  actual  phase.  This  phase  is  referred  to  as  the 
piston-removed  phase  [68].  Secondly,  the  quantity  jj^  +  s  must  be  set  to  zero  in 

Eq.  (56),  since  the  field  is  not  sampled  beyond  the  physical  extent  of  the  aperture  in  a 
single-aperture  system.  Finally,  due  to  conservation  of  energy,  the  irradiance  in  a  single¬ 
aperture  system  is  half  of  {iFi^))  for  two  apertures.  Hence,  the  long  exposure  irradiance 
on  the  focal  plane  of  an  imaging  system,  I{x)  is  given  by: 


I{x)  =  j dArjV{Aff)  exp  J 

xexp{-iD,.[^A?]}, 


I  ^  A  J 

- [r;^” 


(58) 


where  is  implicitly  the  structure  function  of  the  piston-removed  phase  (for  a  single¬ 
aperture  system  this  is  identical  to  the  structure  function  of  the  piston-included  phase). 

Imaging  performance  is  quantified  by  the  Strehl  ratio  (SR),  defined  as  the  ratio  of 
the  irradiance  with  aberrations  present  divided  by  the  diffraction-limited  irradiance  on  the 
optical  axis  [9]: 

SR  =  y  (59) 

l{x  =  0)  ^  ' 

where  I{x  =  0)  is  obtained  from  Eq.  (58)  by  setting  =  0.  For  a  distant  point  source, 
Jp  is  a  constant  [32],  and  SR  is  given  by: 


SR  = 


JdAff  ViArf)  exp{-|D^„  [:^^^]} 


/dA^  ViArf) 


(60) 


Since  the  phase  variance  is  a  stationary  quantity,  the  phase  structure  function  is  given 
explicitly  as: 


D. 


3. 

iRp 


Arf 


R 


=  2  -  2  (^<i>oi‘f^4>o  (r  +  ] ) , 


(61) 


and  SR  can  be  written  as: 


SR  =  exp [- <.2,e-))]  „„ 

A  lower  bound  on  SR  is  obtained  from  Eq.  (62)  by  setting  (^4>oi^^o 
zero  [70, 71]  so  that: 

SR  >  exp  [-  {(t>lir^)]  ■  (63) 
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It  is  convenient  to  express  Eq.  (63)  in  terms  of  the  aperture-averaged  phase  variance, 
defined  as: 

£2  =  JdpWi^{4>l{R^),  (64) 

where  p  represents  a  normalized  aperture  vector  whose  magnitude  is  1  along  the  circum¬ 
ference  of  aperture  1,  W (^  is  the  uniform  aperture  weighting  function  given  by: 

0  otherwise 

and  the  integration  is  over  all  space.  Since  constant  over  the  extent  of  the 

aperture,  and  /  dp  W’(p)  =  1,  then  without  loss  of  generality,  it  may  be  said  that  = 
and  the  lower  bound  on  the  Strehl  ratio  is  expressed  as: 

SR  >  exp(— (66) 

Thus,  as  0,  SR  1,  and  the  long-exposure  irradiance  on  the  optical  axis  approaches 
the  diffraction-limited  irradiance. 

2. 2. 5, 2  Beam-projection  systems.  As  discussed  previously,  beam-projection 
applications  employ  the  same  optical  system  as  imaging  applications.  In  fact,  a  beam- 
projection  system  may  be  thought  of  as  an  imaging  system  run  backwards,  as  indicated 
in  Figure  5.  In  the  case  of  an  imaging  system,  light  propagates  from  a  distant  source  and 
is  focused  onto  a  detector.  In  a  beam-projection  system,  light  propagates  from  a  laser 
source,  and  is  focused  onto  a  distant  target.  The  key  difference  between  these  applications 
lies  in  where  the  turbulence-induced  phase  is  acquired  along  the  optical  path.  In  the 
imaging  application,  the  phase  is  induced  before  entering  the  optical  system.  In  a  beam- 
projection  application,  the  phase  is  induced  after  exiting  the  optical  system.  Thus,  the 
analysis  presented  previously  for  imaging  system  performance  may  not  be  applied  directly 
to  a  beam-projection  system. 

Quantification  of  the  effects  of  turbulence  on  system  performance,  however,  is  related 
for  both  imaging  and  beam  projection  systems.  Lukin  has  rigorously  shown  [49]  that  at 
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Figure  5.  Conventional  optical  systems  for  imaging  (top)  and  beam-projection  applica¬ 
tions  (bottom).  Arrows  indicate  the  direction  of  light  propagation.  Shading 
notionally  shows  the  relative  amount  of  phase  induced  by  propagation  through 
turbulence. 


a  given  time,  the  phase  induced  by  propagation  from  an  aperture  to  a  target  object  is 
conjugate  to  the  phase  induced  by  propagation  from  a  source  to  the  aperture.  In  the 
notation  used  here,  this  implies  that  <l>io  =  If  it  is  assumed  that  most  of  the  phase  is 

induced  near  the  aperture,  then  the  long-exposure  irradiance  at  the  target  object  will  be 
given  by  Eq.  (58),  with  /  replaced  by  the  range  of  the  target  object.  Even  if  the  turbulence 
is  concentrated  a  considerable  distance  along  the  propagation  path,  a  Fourier-transform 
relationship  exits  between  the  field  at  that  point  and  the  field  at  the  focal  point  which  is 
similar  to  Eq.  (40),  but  with  scaled  arguments  [31].  Thus,  the  analysis  of  the  long-exposure 
irradiance  at  the  target  object  can  be  carried  out  as  for  the  imaging  system,  and  will  yield 
the  same  lower  bound  for  the  Strehl  ratio.  Hence,  Eq.  66  can  also  be  used  to  quantify  the 
effect  of  turbulence-induced  phase  in  a  beam-projection  system. 
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2. 2, 5, 3  Interferometric  systems.  The  performance  of  an  interferometric 
system  can  be  derived  directly  from  Eq.  (56),  considering  that  for  interferometric  applica¬ 
tions,  it  is  not  necessary  to  spatially  resolve  the  irradiance  on  the  focal  plane.  Instead,  the 
long-exposure  irradiance  is  integrated  over  the  extent  of  a  single-element  detector,  yielding 
the  integrated  irradiance  I  given  by: 


1  = 


(67) 


Substituting  Eq.  (56)  into  Eq.  (67)  and  exchanging  the  order  of  integration  yields  the 
following  expression  for  /: 


/  = 


J dAi]  'P{Aff)  ®  I 


X  < 


h  exp  } 

+  I  Jo  [^Aif  +  s]  I  cos  {^h  +  -0  ')  exp  +  *] }  > 


(68) 


Performing  the  integration  over  x  yields  a  Dirac  delta  function,  which  makes  the  integration 
over  Afj  trivial,  resulting  in  the  following  expression  for  I: 


=  7o  |l  +  |Ac.(^|  cos  +  Ipj  exp 


(69) 


where  Iq  =  2(A/)“^'P(0)Jo(0),  the  magnitude  of  the  normalized  mutual  intensity  |/io(^|  = 
|Jg(^|/Jo(0),  and  •0  =  2irg[J^(^].  Hence,  at  a  fixed  baseline  s,  I  is  modulated  sinusoidally 
as  the  path  difference  h  is  varied.  This  intensity  modulation  is  commonly  referred  to  as 
the  fringe.  Using  measurements  of  the  irradiance  modulation,  the  fringe  visibility  V  is 
calculated  as  [32]: 


^  _  max(/)  —  min(7) 
~  max(7)  +  min(7) 

Using  Eq.  (69),  Eq.  (70)  implies  that  V  is  given  by: 


(70) 


V  =  |Ao(^|exp 


(71) 
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which  says  that  the  fringe  visibility  for  long-exposure  irradiance  measurements  in  an  in¬ 
terferometer  is  determined  by  the  magnitude  of  the  normalized  mutual  intensity  and  by 
the  value  of  the  phase  structure  function,  both  evaluated  at  the  baseline  separation  s. 
If  turbulence  effects  are  negligible,  then  0,  and  the  fringe  visibility  for  free-space 

propagation  is  given  by  F  =  |/io(^|.  In  analogy  with  the  Strehl  ratio  definition  of  Eq.  (59), 
the  fringe  Strehl  ratio  (FSR)  may  be  defined  as: 


V 

FSR  =  —  =  exp 

V 


(72) 


As  indicated  in  Eq.  (72),  the  effect  of  turbulence-induced  phase  aberrations  in  an  interfer¬ 
ometric  system  is  to  reduce  fringe  visibility  beyond  what  would  be  expected  if  aberrations 
were  not  present  in  the  system. 

Since  the  apertures  are  separated  by  s,  then  D^^{s)  can  be  calculated  by  measuring 
the  phase  in  each  aperture  at  corresponding  (but  otherwise  arbitrary)  locations: 


D,t,A^  =  ([<^oi(f)  -  ■ 


(73) 


Consider  the  aperture-averaged  differential-phase  variance  6^  defined  as: 

6^  =  1  dpWi^  .  (74) 

Since  ([(/>oi(RpO  “  is  a  spatially-stationary  quantity,  then: 

=  -  MR^]^) 

=  (75) 


Hence,  the  fringe  Strehl  ratio  may  be  written  as: 

FSR  =  exp  (-i  6^^  . 


(76) 
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Thus,  as  ^  0,  then  FSR  ^  1,  and  the  performance  of  an  interferometric  system 
approaches  that  which  is  expected  for  free-space  propagation. 

2.3  Adaptive- optical  systems 

In  Section  2.2,  it  was  shown  that  the  effect  of  turbulence-induced  phase  is  to  degrade 
optical  system  performance.  The  performance  degradation  is  dependent  upon  the  aperture- 
averaged  phase  variance  (differential  phase  variance  for  interferometric  systems).  Reducing 
the  phase  variance  will  enhance  the  performance  of  the  system.  To  reduce  the  phase 
variance,  the  turbulence-induced  phase  must  be  compensated.  An  adaptive  optics  system 
removes  phase  aberrations  through  reflection  from  a  surface  which  is  phase-conjugate  to 
the  incident  wave-front.  The  resulting  wave-front  is  (ideally)  planar,  and  the  performance 
of  the  optical  system  is  restored  to  its  diffraction  limit  [69]. 

2.3.1  System  components  and  description.  Adaptive-optical  compensation  re¬ 
quires  that  the  turbulence-induced  phase  can  be  sensed,  and  that  the  sensed  phase  can  be 
removed  by  a  suitable  device.  The  schematic  diagram  in  Figure  6  depicts  a  typical  AO 
system.  In  this  configuration,  the  compensating  device  is  a  deformable  mirror  (DM).  Upon 
reflection  from  the  mirror,  light  passes  through  a  beam  splitter.  Part  of  the  light  is  focused 
onto  the  focal  plane  of  the  system,  while  the  other  part  is  passed  onto  the  wave-front  sen¬ 
sor  (WFS),  which  measures  the  turbulence-induced  phase.  Reconstructed  phase  from  the 
WFS  is  used  to  drive  the  mechanical  actuators  of  the  DM.  As  depicted  in  Figure  6,  the 
system  is  closed-loop,  with  the  WFS  measuring  the  residual  phase  after  partial  conjugation 
by  the  DM. 

There  are  two  main  types  of  wave-front  sensors  used  in  AO  systems.  These  two  types 
are  referred  to  as  Hartmann  wave- front  sensors  and  shearing  interferometers  [68].  While 
the  instrumentation  used  for  these  two  WFS  types  is  different,  both  are  sensitive  to  the 
local  gradient  of  the  turbulence-induced  phase.  Thus,  measurements  from  either  type  of 
instrument  can  be  employed  in  a  similar  manner. 

There  are  also  two  main  types  of  deformable  mirrors.  Deformable  mirrors  are  gen¬ 
erally  classified  as  being  either  segmented  or  continuous.  A  segmented  mirror  is  made  up 
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source 

object 


Figure  6.  Typical  adaptive  optics  system.  The  turbulence-induced  phase  is  measured  by 
the  wave-front  sensor,  and  the  conjugate  phase  is  placed  on  the  deformable 
mirror.  When  the  deformable  mirror  is  placed  before  the  wave-front  sensor  as 
shown  here,  the  system  is  closed-loop,  and  the  wave-front  sensor  detects  the 
residual  phase  after  partial  conjugation  by  the  deformable  mirror. 
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of  an  array  of  closely-spaced,  individually  controlled  mirror  segments  that  are  driven  by 
mechanical  actuators.  Continuous  deformable  mirrors  have  a  continuous  reflective  surface 
called  a  face  plate  that  can  be  pushed  or  pulled  from  behind  by  mechanical  actuators  to 
form  a  surface  which  is  phase-conjugate  to  the  incident  wave-front. 

2.3.2  Natural- guide-star  adaptive  optics.  Ideally,  the  object  of  interest  will  be 
bright  enough  so  that  an  accurate  phase  measurement  can  be  obtained  by  the  AO  system. 
If  this  is  not  the  case,  then  another  source  must  be  employed.  Often,  a  naturally-occurring 
nearby  star  is  used,  and  is  referred  to  as  a  natural  guide  star  or  natural  beacon.  Because 
the  maximum  permissible  angular  separation  between  the  object  of  interest  and  the  NGS 
at  visible  wavelengths  is  a  few  arc-seconds,  the  probability  of  acquiring  a  NGS  of  sufficient 
brightness  at  a  particular  region  in  the  sky  is  small  [27].  Therefore,  NGS  adaptive  optics 
may  only  be  used  when  imaging  bright  objects  or  when  the  object  of  interest  is  near  a  suf¬ 
ficiently  bright  source.  Nonetheless,  a  NGS  is  often  used  for  wave-front  tilt  compensation, 
since  the  maximum  angular  separation  is  larger  for  wave-front  tilt  than  for  higher-order 
phase  components  [71]. 

2.3.3  Laser-guide-star  adaptive  optics.  To  overcome  the  sky  coverage  limitations 
imposed  on  NGS  adaptive  optics,  an  artificial  beacon  or  laser  guide  star  is  used  for  adaptive- 
optical  imaging.  The  use  of  a  LGS  in  adaptive  optics  was  first  suggested  in  the  open 
literature  by  Foy  and  Labeyrie  in  1985  [18].  The  LGS  is  formed  by  projecting  a  focused 
laser  beam  into  the  atmosphere  [25,56,80],  where  a  beacon  is  formed  by  scattering  or 
resonance  processes  in  certain  atmospheric  layers.  The  optical  system  used  to  transmit 
the  beam  into  the  atmosphere  is  usually  the  same  as  the  system  used  to  receive  the  beacon 
light.  Since  the  transmitting  and  receiving  apertures  are  the  same,  the  LGS  application 
is  referred  to  as  monostatic.  If  different  apertures  are  used  for  transmitting  and  receiving, 
then  the  LGS  application  is  referred  to  as  bistatic.  The  design  of  AO  systems  employing 
laser  guide  stars  has  been  studied  in  detail  [27]. 

Laser  guide  star  applications  are  classified  by  the  processes  involved  in  beacon  forma¬ 
tion.  Rayleigh  beacons  are  formed  by  Rayleigh  scattering  from  air  molecules  at  altitudes 
from  10  km  to  30  km.  Sodium  (Na)  beacons  are  formed  by  the  resonance  fluorescence  of 
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sodium  atoms  in  the  mesosphere  [27]  at  an  altitude  of  approximately  90  km.  The  per¬ 
formance  of  AO  systems  employing  either  Rayleigh  or  Na  beacons  depends  upon  many 
factors,  most  notably  laser  power  requirements  and  anisoplanatism  [93]. 

As  was  originally  pointed  out  by  Pilkington  [55],  the  LGS  cannot  directly  provide  for 
tilt  compensation.  When  the  LGS  beacon  is  formed,  the  angular  deflection  of  the  beam 
due  to  turbulence  results  in  a  spatial  displacement  of  the  beacon  which  preserves  the  tilt 
induced  on  the  upward  path.  All  higher-order  aberrations  are  not  preserved  in  this  process. 
Since  the  beacon  light  propagates  back  to  the  system  aperture  along  the  same  path  as  the 
upward  propagation,  the  preserved  wave-front  tilt  is  naturally  conjugated.  The  result  is  a 
null  tilt  measurement  with  only  higher-order  aberrations  present  in  the  LGS  signal.  This 
effect  has  been  referred  to  in  the  literature  as  beam  reciprocity  [5-7].  Due  to  the  lack  of 
tilt  information,  the  LGS  signal  is  conventionally  used  only  for  higher-order  wave-front 
compensation  [23]  with  a  nearby  NGS  providing  the  tilt  correction  [24]. 

2,3,4  Zernike  modal  compensation  of  phase  aberrations.  The  turbulence-induced 
phase  that  has  been  discussed  in  the  preceding  sections  varies  randomly  in  space  and  time. 
In  order  to  characterize  the  spatial  variation  of  the  phase  at  a  given  time,  it  is  useful  to 
decompose  the  phase  with  respect  to  a  set  of  basis  functions.  Since  the  relevant  quantity  is 
the  phase  within  a  circular  aperture,  Zernike  polynomials  are  an  appropriate  basis  set  [9]. 
Zernike  polynomials,  designated  Zi{^^  are  defined  on  a  circle  of  unit  radius,  with  the 
normalized  polar  vector  p  specifying  a  location  in  the  plane  of  the  circle.  The  subscript  i 
designates  a  particular  Zernike  mode.  In  this  research,  we  follow  the  ordering  convention  of 
Noll  [54],  where  each  subscript  i  designates  an  ordered  pair  (ni,mi),  specifying  the  radial 
and  azimuthal  order  of  the  polynomial.  With  this  ordering  scheme,  the  orthonormality  of 
the  basis  set  can  be  written  as: 


I  dpW{^Z,{^Zj{^  =  6ij,  (77) 

where  Sij  is  the  Kronecker  delta  and  W{^  is  the  aperture  weighting  function  given  in 
Eq.  (65). 
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Using  Noll’s  ordering  method,  the  Zernike  polynomials  are  defined  as: 


Zi{p,e) 

=  y/rii  +  1  (p)  cos  mi 0,  i  even,  rrii  ^  0 

(78) 

Zi{p,0) 

=  y/ni  + 1  R^.'ip)V2smmi6,  i  odd,  rrii  ^  0 

(79) 

Zi{p,9) 

=  ■v/ni  +  li?°;(p),  mi  =  0, 

(80) 

where  the  radial  function  is  given  by: 


(ni-mO/2 

~  §  g![(ni  +  mi)/2-g]![(ni-mi)/2-g]!' 


(-l)«(ni-g)! 


Jii-2q 


(81) 


The  values  of  n*  and  rui  are  always  integers  and  satisfy  <  rij,  |m^|  =  even.  The 
Zernike  polynomials  through  Noll’s  order  i  =  21  along  with  the  corresponding  values  of  rii 
and  rrii  are  listed  in  Table  2. 

With  the  Zernike  polynomials  as  a  basis  set,  the  turbulence-induced  phase  from  the 
object  of  interest  in  an  aperture  of  radius  R  may  be  expanded  as  follows: 


(82) 

% 

where  a^i  is  the  coefficient  obtained  by  projection  of  the  phase  onto  the  zth  Zernike  poly¬ 
nomial: 

=  I  (83) 

In  addition  to  being  a  complete  set  on  the  unit  circle,  the  low-order  Zernike  polynomials 
correspond  to  classic  fixed  aberrations  in  optical  systems.  For  instance,  Zi{^  represents 
a  piston  aberration.  The  orthogonal  components  of  tilt  (x-tilt  and  y-tilt  respectively),  are 
represented  by  Z2{^  and  Z^i^-  Of  the  higher-order  polynomials,  Z4{^  represents  defo¬ 
cus,  Zs{^  and  Z^{^  represent  orthogonal  components  of  astigmatism,  Zr{^  and  Zs{^ 
represent  orthogonal  components  of  coma,  and  Zii{p)  represents  spherical  aberration.  The 
corresponding  coefficients  aoi  indicate  the  amount  of  each  type  of  aberration  present  in  the 
turbulence-induced  phase.  Graphical  representations  of  several  low-order  Zernike  polyno¬ 
mials  are  shown  in  Figure  7.  From  the  definition  of  aoi  in  Eq.  (83),  it  is  clear  that  since 
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Zernike2:  x-tllt 


Zernike  3: 


Figure  7.  Graphical  representations  of  low-order  Zernike  polynomials.  Light  areas  corre¬ 
spond  to  the  highest  values  of  Zi{p)^  dark  areas  correspond  to  the  lowest. 
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Noll’s 

order 

i 

radial 

order 

rii 

azimuthal 

order 

rrii 

Zernike  polynomial  Zi{p^  0) 

1 

0 

0 

1 

2 

1 

1 

2pcos6 

3 

1 

1 

2psm0 

4 

2 

0 

3.464/92  -  1.732 

5 

2 

2 

2.449p2sin20 

6 

2 

2 

2.449p2cos20 

7 

3 

1 

(8.485p2  -5.657p)sin0 

8 

3 

1 

(8.485p2  _  5.657p)  cos  0 

9 

3 

3 

2.828p2sin30 

10 

3 

3 

2.828p^  cos  30 

11 

4 

0 

3.416p^  -  13.416p2  +  2.236 

12 

4 

2 

(12.649p^  -9.487p2)cos20 

13 

4 

2 

(12.649p^  -9.487p2)sm26> 

14 

4 

4 

3.162/9^  cos  40 

15 

4 

4 

3.162p'‘ sin40 

16 

5 

1 

(34.641p®  -  41.569p2  +  10.392p)  cos  0 

17 

5 

1 

(34.641p®  -  41.569p2  +  lO.392p)sin0 

18 

5 

3 

(17.321p®  -  13.856p2)cos30 

19 

5 

3 

(17.321p®  -  13.856p2)sin30 

20 

5 

5 

3.464p^  cos  50 

21 

5 

5 

3.464p®sin50 

Table  2.  Zernike  polynomials  through  Noll’s  order  i  =  21.  The  radial  and  azimuthal 
orders  of  the  ith  polynomial  are  specified  by  the  ordered  pair  {rti^rrii). 


(po  is  zero-mean  Gaussian,  aoi  must  also  be  zero-mean  Gaussian. 

The  Zernike  phase  expansion  technique  has  been  used  extensively  in  the  analysis 
of  atmospheric  turbulence  [38,54,63],  and  is  often  referred  to  as  modal  analysis.  Modal 
analysis  has  been  applied  to  wave-front  correction  in  adaptive-optical  systems,  since  the 
effects  of  modal  compensation  are  related  directly  to  the  statistical  properties  of  the  Zernike 
expansion  coefiicients  [14,92].  To  see  this,  consider  the  effect  of  applying  a  conjugate  phase 
^0  to  the  deformable  mirror  in  an  AO  system,  for  represented  by  a  finite  series  of  modal 
components  through  order  N: 

N 

~4>o{R^  =  Y,aoiZi{p),  (84) 

i=2 
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where  is  the  ith  Zernike  coefficient  of  the  conjugate  phase.  Thus,  the  iVth  order 
aperture-averaged  residual  phase  variance  is  calculated  according  to  the  definition  of 
given  in  Eq.  (64),  with  replaced  by  “  <^0* 


=  jdpW{^ 


n  2  V 


Mm 


(85) 


By  employing  the  orthonormality  property  of  the  Zernike  polynomials  on  the  unit  circle 
[Eq.  (77)],  6%  may  be  expressed  as: 


f  (^  {4>1{R^)  -  ^  {ali)  +  ^  {{a^i  -  Kif)  .  (86) 

^  i=2  i=2 

In  an  ideal  AO  system,  Uoi  =  Ooi,  and  the  third  term  of  Eq.  (86)  is  zero,  such  that: 


i=2 

Thus,  the  aperture-averaged  phase  variance  is  reduced  by  the  sum  of  the  modal  phase 
coefficient  variances  through  A^th  order.  Since  is  a  positive  quantity  for  all 
then  always  decreases  with  increasing  N, 

For  Kolmogorov  turbulence,  the  aperture-averaged  piston-removed  phase  variance 
for  an  astronomical  source  object  is  given  by  [54]: 

I  dp  W{^  {MRp))  =  1-0299  (0'^'  .  (88) 


The  Kolmogorov  modal  covariance  of  and  ttoj,  designated  (ao^aoj),  is  zero  for  {i—j)  odd, 
and  for  {i  —  j)  even  is  given  by  [54,92]: 

/  n  \  5/3 

{aoiaoi)  =  0.0072  (-l)("‘+«>-2-0/2[(^.  + 

^ _ r(14/3)r[(n,+n,-5/3)/2] _ 

r[(ni  -  rij  +  17/3)/2]r[(n,-  -  +  17/3)/2]r[(ni  +  +  23/3)/2]  ’  ^  ’ 

Table  3  shows  the  value  of  {agiaoj)  {Dfro)~^^^  for  i,j  =2  —  10.  Using  Eq.  (88)  and  Eq.  (89), 
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4 

5 

6 

■B 

8 

9 

10 
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0.448 

0 

0 

H9I 

Bl 

0 

0.448 

0 

BEH 

II^^H 

mm 

0 

0 

0.023 

0 

0 

B 

0 

0 

0 

0.023 

H9I 

0 

0 

B 

0 

0 

0 

0 

0.023 

0 

0 

0 

0 

B 

0 

-0.014 

0 

0 

8 

0 

0 

mm 

0 

9 

0 

0.006 

0 

10 

HOI 

iBEH 

HEI^I 

0 

0.006 

Table  3.  Covariance  matrix  of  low-order  Zernike  coefficients.  Table  shows  the  value  of 
{aoiaoj)  for  =  2-10,  with  rows  representing  i,  columns  repre¬ 

senting  j. 

for  an  ideal  AO  system  can  be  computed  from  Eq.  (87).  For  an  imaging  system,  Strehl 
performance  is  determined  by  replacing  with  in  Eq.  (66).  Thus,  for  iVth-order 
compensation,  SR  is  determined  by: 


SR  >  exp(-£:^).  (90) 

Since  is  proportional  to  (D/ro)^/^  for  all  AT,  the  value  of  the  relative  aperture  diameter 
D/r^  will  determine  the  Strehl  performance  of  an  imaging  system  for  a  given  order  of  modal 
compensation.  Figure  8  shows  the  Strehl  ratio  improvement  as  a  function  of  compensation 
order  N  for  D/tq  =  5,10,15,20,25,30,  beginning  with  iV  =  2,  and  continuing  through 
N  =  105  (all  modes  with  <  13).  As  is  shown  in  Figure  8,  without  AO  compensation, 
SR  is  near  zero  for  all  values  of  D/tq  considered.  For  D/r^  =  5,  nearly  diffraction- 
limited  performance  can  be  achieved  by  an  ideal  AO  system  with  modal  compensation 
through  N  =  15,  For  D/tq  =  30,  however,  SR  is  appreciably  greater  than  zero  only  for 
N  >  25.  Thus,  the  order  of  modal  compensation  required  in  an  ideal  AO  system  for  nearly 
diffraction-limited  performance  depends  upon  the  size  of  the  system  aperture,  relative  to 
the  atmospheric  coherence  diameter. 
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1 


Figure  8.  Effect  of  modal  compensation  order  N  on  Strehl  ratio  for  the  relative  aperture 
diameter  D/ro- 

2,4  Anisoplanatism  in  adaptive  optics 

As  was  mentioned  in  Section  2.3,  an  ideal  adaptive  optics  system  uses  light  from  the 
object  under  observation  for  wave- front  sensing,  and  compensation  performance  is  given 
by  the  results  shown  in  Figure  8.  When  a  source  other  than  the  object  being  observed 
is  employed  in  wave-front  sensing,  however,  compensation  performance  is  degraded.  This 
effect  is  referred  to  as  anisoplanatism.  As  it  applies  to  adaptive  optics,  the  term  “anisopla¬ 
natism”  is  borrowed  from  the  terminology  for  fixed  optical  systems,  where  the  isoplanatic 
region  is  defined  as  that  region  in  which  the  imaging  system  point  spread  function  can  be 
considered  translation  invariant  [9].  Anisoplanatism  in  adaptive-optical  systems  implies 
that  the  point  spread  function  of  the  fixed  optical  system  plus  the  atmosphere  is  not  trans¬ 
lation  invariant  [45].  As  pointed  out  by  Fried  [21],  the  term  “anisoplanatism”  refers  to  a 
variety  of  effects  brought  on  by  the  stochastic  nature  of  atmospheric  turbulence. 
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In  an  AO  system,  anisoplanatism  occurs  when  the  turbulence-induced  phase  mea¬ 
sured  by  the  WFS  does  not  match  the  turbulence-induced  phase  incident  in  the  compensa¬ 
tion  aperture  from  the  object  of  interest.  This  occurs  whenever  light  from  the  WFS  beacon 
and  the  object  sample  different  volumes  of  turbulence.  Anisoplanatic  effects  are  present  if 
there  is  a  spatial  separation  between  the  object  and  beacon,  a  spatial  separation  between 
the  wave-front-sensing  aperture  and  the  wave-front-compensation  aperture,  or  when  time 
delays  in  the  system  cause  the  beacon  phase  and  the  object  phase  to  be  only  partially 
correlated  due  to  atmospheric  winds  or  motion  of  the  system  components.  Therefore,  in 
adaptive  optics,  anisoplanatism  is  a  spatial  and  temporal  effect  present  in  every  system  to 
some  degree. 

Conventionally,  anisoplanatic  effects  are  classified  into  two  types:  angular  and  focal 
anisoplanatism.  Angular  anisoplanatism  arises  when  the  beacon  has  an  angular  displace¬ 
ment  Of,  in  the  field  of  view  from  the  object  being  imaged,  as  is  depicted  in  Figure  9.  One 
of  the  principle  advantages  of  the  LGS  over  the  NGS  is  that  it  can  be  positioned  anywhere 
in  the  field  of  view.  Thus,  the  angular  separation  between  the  beacon  and  object  may  be 
taken  to  zero.  In  this  case,  an  AO  system  may  still  suffer  from  focal  anisoplanatism,  as  the 
beacon  is  located  at  a  lower  altitude  than  the  object.  Hence,  light  from  the  object  passes 
through  more  turbulent  atmosphere  than  the  light  coming  from  the  beacon.  An  example 
of  focal  anisoplanatism  in  an  AO  system  is  shown  in  Figure  9. 

Fried  quantified  the  general  effects  of  angular  anisoplanatism  on  fixed  imaging  and 
beam-projection  systems  and  characterized  the  magnitude  of  the  angular  anisoplanatic 
effect  by  the  parameter  Oq,  which  is  referred  to  as  the  isoplanatic  patch  [22]  or  isoplanatic 
angle,  and  is  calculated  as: 


0a  = 


/0-\2  -L 

2.905  (  —  )  /  dzCl{z)z^l^ 

V  A  /  Jo 


1  -3/5 


(91) 


For  the  HV-21  profile  in  Table  1,  0o  —  7  /xrad  (1.4  arcsec).  Since  the  HV-21  profile  yields 
ro  =  5  cm  and  0o  =  7  /xrad,  the  HV-21  profile  is  often  referred  to  as  HV5/7.  Fried  defined 
the  isoplanatic  angle  such  that  for  Of,  =  0^,  the  aperture- averaged  residual  phase  variance 
is  1  rad^,  which  implies  that  SR  0.4.  Fried’s  definition  of  0o  has  been  shown  to  be  a 
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Figure  9.  Examples  of  angular  and  focal  anisoplanatism  in  an  adaptive  optics  system. 
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pessimistic  indicator  of  AO  system  performance,  as  it  includes  the  effect  of  the  piston  phase 
component  in  quantifying  the  isoplanatic  angle  [74].  Since  imaging  system  performance 
is  insensitive  to  piston  phase,  actual  angular  anisoplanatic  degradation  is  less  severe  than 
predicted  by  Fried’s  analysis.  Further  analysis  efforts  have  explored  the  degrading  effects 
of  various  forms  of  spatial  and  temporal  anisoplanatism  on  wave-front  compensation  in 
adaptive-optical  systems  [23, 82-84, 87, 93] . 

Anisoplanatic  effects  can  also  be  quantified  using  modal  analysis  techniques.  Over 
the  past  several  years,  many  works  have  been  published  in  which  the  correlative  properties 
of  modal  coefficients  for  various  forms  of  anisoplanatism  have  been  quantified  [17,51,62,74, 
76,94].  When  examining  wave-front  compensation,  the  authors  of  these  studies  have  used 
modal  correlation  calculations  in  order  to  compute  the  aperture-averaged  residual  phase 
variance,  thus  establishing  limits  for  effective  compensation  when  system  performance  is 
degraded  by  turbulence-induced  anisoplanatism. 

To  quantify  anisoplanatism  in  modal  compensation,  consider  the  aperture-averaged 
residual  phase  variance  when  the  phase  from  the  beacon,  designated  (j>t,  is  used  as  the 
conjugate  phase.  In  this  case,  Eq.  (86)  is  written  as: 

f  dpW(?)  (rAm)  -  Y,  (4>  +  Y  ((«.<  -  «n)’) ,  (92) 

i=2  i=2 

where  is  the  ith  Zernike  coefficient  of  the  beacon  phase: 

au  =  J  (93) 

Since  (l>t  is  only  partially  correlated  with  (po,  then  ati  ^  a^i,  and  ((a^j  -  abif)  >  0  for  all  i. 
It  follows  that: 

N 

^  V  {(^oi  ^bi)  )  ^  0,  (94) 

i=2 

and  since  this  term  is  additive  in  Eq.  (92),  is  always  increased  when  anisoplanatism  is 
present  in  the  system,  leading  to  a  decreased  value  for  the  Strehl  ratio. 
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The  expression  for  in  Eq.  92  may  be  simplified  by  expanding  the  third  term  to 


yield: 

N 

(95) 

i=2 

Prom  Eq.  (95),  it  is  apparent  that  the  value  of  {a^iai^i)  (the  correlation  of  the  object  and 
beacon  phase  coefficients)  determines  the  magnitude  of  the  anisoplanatic  effect.  Prom  the 
analysis  methods  of  Stone  et  al  [74]  and  Molodij  and  Rousset  [51],  may  be  quantified 
according  to  Eq.  95.  The  corresponding  Strehl  ratio  is  shown  in  Pigure  10  for  modal  com¬ 
pensation  through  AT  =  6, 10, 15, 21,  with  I)/ro  =  5  and  the  HV-21  profile.  Pigure  10(a) 
shows  the  degradation  of  SR  with  increasing  NGS  beacon  offset  angle,  which  has  been  nor¬ 
malized  to  Pried’s  isoplanatic  angle  given  in  Eq.  (91).  These  results  demonstrate  that  for 
piston-removed  phase,  the  Strehl  ratio  may  be  significantly  larger  than  0.4  at  9^16^  =  1. 
Por  modal  compensation  at  a  particular  NGS  offset,  anisoplanatism  may  be  partially  over¬ 
come  by  increasing  the  order  of  correction,  N,  Pigure  10(b)  shows  the  degradation  of  SR 
with  decreasing  LGS  altitude  due  to  focal  anisoplanatism.  As  the  LGS  altitude  increases, 
system  performance  approaches  that  of  an  NGS  adaptive  optics  system  with  0^  =  0.  Al¬ 
though  the  Strehl  ratio  calculation  has  been  performed  over  the  entire  range  of  LGS  altitude 
from  10  km  to  100  km,  LGS  beacons  may  only  be  formed  over  a  limited  subset  of  these 
ranges.  Rayleigh  beacons  can  be  formed  from  10  km  to  30  km,  whereas  sodium  beacons 
are  formed  near  90  km  altitude.  The  results  shown  in  Pigure  10(b)  illustrate  the  advantage 
of  the  sodium  LGS  over  the  Rayleigh  LGS  due  to  decreased  focal  anisoplanatism. 

Although  anisoplanatism  takes  many  forms,  its  effect  is  universal;  that  is,  turbulence- 
induced  anisoplanatism  degrades  the  performance  of  an  adaptive  optics  system.  While  the 
forms  of  anisoplanatism  are  diverse,  all  forms  are  quantified  by  the  correlative  properties  of 
the  turbulence-induced  phase.  Any  AO  system  in  which  the  phase  incident  upon  the  system 
aperture  is  only  partially  correlated  with  the  conjugate  phase  placed  on  the  deformable 
mirror  will  suffer  anisoplanatic  effects.  Therefore,  instead  of  investigating  a  particular 
form  of  anisoplanatism  (angular,  focal,  temporal,  etc.)  as  it  applies  to  a  given  system,  this 
work  concentrates  instead  on  the  manifestation  of  anisoplanatism  (broadly  defined)  in  that 
system,  and  its  effect  on  system  performance. 


48 


(a) 


(b) 

Figure  10.  Degradation  of  Strehl  ratio  due  to  anisoplanatism.  (a)  Angular  anisoplanatism 
effect  on  Strehl  ratio  for  normalized  NGS  beacon  offset  6^,1  Oq,  (b)  Focal  aniso¬ 
planatism  effect  on  Strehl  ratio  as  LGS  beacon  altitude  increases.  Results 
shown  are  for  modal  compensation  through  iV  =  6, 10, 15,21,  with  D/tq  =  5 
using  the  HV-21  profile. 
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2. 5  Summary 

This  chapter  provides  the  background  knowledge  in  atmospheric  and  adaptive  optics 
necessary  to  understand  the  analysis  and  application  chapters  which  follow.  The  effects 
of  turbulence  in  optical  systems  were  discussed,  beginning  with  fundamental  considera¬ 
tions  regarding  the  stochastic  nature  of  index  of  refraction  fluctuations  in  the  atmosphere, 
which  is  characterized  by  the  spatial  PSD  Wn'(k),  By  considering  the  effect  of  small  per¬ 
turbations  in  the  index  of  refraction  on  wave  propagation,  it  is  was  shown  that  a  random 
phase  is  induced  by  turbulence  as  an  electromagnetic  wave  travels  through  the  atmosphere. 
The  stochastic  nature  of  the  turbulence-induced  phase  is  also  specified  by  its  spatial  PSD 
W^(k)^  which  is  calculated  directly  from  Wn>{k),  Equivalently,  the  statistical  properties  of 
turbulence- induced  phase  may  be  quantified  by  the  phase  structure  function  D^{s)^  which 
indicates  the  mean-square  difference  in  phase  for  points  separated  by  a  distance  s.  It  was 
shown  that  Dfp{s)  scales  with  the  atmospheric  coherence  diameter  Tq,  which  is  used  to 
quantify  the  magnitude, of  turbulence  effects  in  optical  systems  according  to  the  relative 
aperture  diameter  D/tq. 

Prom  the  basic  results  of  frequency  analysis  of  optical  systems,  a  framework  was 
developed  in  which  turbulence  effects  in  imaging,  beam-projection,  and  interferometric 
systems  can  be  understood.  It  was  shown  that  the  performance  of  imaging  and  beam- 
projection  systems  may  be  quantified  by  a  Strehl  ratio  (SR),  and  that  SR  >  exp(— 
where  is  the  aperture-averaged  phase  variance.  For  interferometric  systems,  a  fringe 
Strehl  ratio  (FSR)  was  defined,  which  is  given  by  FSR  =  exp(— 5^/2),  where  6^  is  the 
aperture-averaged  differential  phase  variance  between  the  baseline  apertures. 

Adaptive-optical  systems  were  discussed,  and  the  components  common  to  most  AO 
systems  were  described,  reviewing  issues  related  to  the  use  of  natural  and  laser  guide  stars. 
The  concept  of  modal  compensation  in  AO  systems  was  presented,  illustrating  that  as  D/tq 
increases,  higher-order  modal  compensation  is  required  to  maintain  system  performance. 
The  performance  of  an  AO  system  is  also  degraded  by  turbulence-induced  anisoplanatism. 
For  NGS  adaptive  optics,  angular  anisoplanatism  reduces  the  Strehl  ratio  with  increasing 
angular  separation  between  the  beacon  and  the  object  of  interest.  For  LGS  adaptive 
optics,  focal  anisoplanatism  reduces  the  Strehl  ratio  as  the  beacon  altitude  is  decreased. 
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In  modal  compensation  applications,  anisoplanatism  requires  higher-order  correction  in 
order  to  maintain  performance.  Anisoplanatism  affects  all  non-ideal  AO  systems  to  some 
degree,  since  it  arises  whenever  the  phase  from  the  object  of  interest  and  the  conjugate 
phase  placed  on  the  deformable  mirror  are  only  partially  correlated.  In  the  context  of 
modal  analysis,  anisoplanatic  effects  are  quantified  by  the  correlation  of  beacon  and  object 
phase  coefficients. 
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III.  Optimal  Compensation  for  Anisoplanatic  Systems 

3. 1  Introduction 

This  chapter  builds  upon  the  background  in  atmospheric  and  adaptive  optics  pre¬ 
sented  in  Chapter  II  through  an  in-depth  analysis  of  atmospheric  compensation  in  AO 
systems  degraded  by  turbulence-induced  anisoplanatism.  This  analysis  begins  in  Section 

3.2  by  considering  the  compensation  performance  of  an  adaptive-optical  interferometric 
system.  Results  obtained  by  considering  this  system  are  also  applicable  to  imaging  and 
beam-projection  systems.  When  anisoplanatism  is  present  in  an  interferometric  system, 
compensation  performance  is  shown  to  be  degraded  from  the  limiting  performance  of  an 
ideal  system.  Using  optimal  compensation,  part  of  the  performance  degradation  may  be 
recovered.  In  order  to  construct  the  optimal  estimator  of  conjugate  phase,  the  modal  cor¬ 
relation  of  object  and  beacon  wave-front  phase  must  be  quantified.  An  analysis  geometry 
is  introduced  in  Section  3.3  to  model  the  most  general  form  of  anisoplanatism  required 
for  imaging,  beam-projection,  and  interferometric  system  analysis.  For  this  generalized 
geometry,  the  cross-correlation  of  the  object  and  beacon  wave-front  phase  coefficients  is 
calculated.  The  resulting  cross-correlation  expression  is  then  validated  by  comparison  with 
published  theoretical  results  and  by  reduction  and  analysis  of  laboratory  phase  measure¬ 
ments.  The  chapter  is  summarized  in  Section  3.4. 

3.2  Optimal  modal  compensation  in  adaptive  optics 

The  discussion  of  atmospheric  compensation  in  Chapter  II  focused  on  the  perfor¬ 
mance  of  an  ideal  AO  system,  demonstrating  how  performance  is  degraded  when  aniso¬ 
planatism  is  present.  Anisoplanatism  causes  beacon  phase  coefficients  to  be  only  partially 
correlated  with  the  phase  coefficients  for  the  object  of  interest.  As  a  result,  the  residual 
phase  after  wave-front  correction  is  increased,  resulting  in  reduced  system  performance.  If 
the  correlation  of  object  and  beacon  phase  coefficients  can  be  quantified  for  the  form  of 
anisoplanatism  present  in  the  system,  then  optimal  estimation  theory  [43]  can  be  applied  to 
the  wave-front  compensation  problem.  In  the  context  of  optimal  estimation  theory,  beacon 
phase  measurements  are  not  used  directly,  but  rather  the  wave-front  correction  applied  to 


52 


aperture  1 
radius  R 


aperture  2 
radius  R 


Figure  11.  Schematic  diagram  of  an  AO  system  for  an  optical  interferometer.  If  all  system 
components  associated  with  aperture  2  are  removed,  the  remaining  compo¬ 
nents  represent  a  conventional  AO  system  for  imaging  and  beam-projection 
applications. 

the  DM  is  estimated  from  the  beacon  phase  measurements,  given  the  correlation  of  the 
beacon  and  object  phase.  The  result  is  a  partial  recovery  of  the  performance  degradation 
due  to  anisoplanatism. 

3.2.1  Wave-front  compensation  performance.  The  analysis  of  compensation  per¬ 
formance  for  anisoplanatism  in  AO  systems  in  the  context  of  optimal  estimation  theory 
proceeds  in  parallel  with  the  optical  system  performance  analysis  of  Subsection  2.2.5.  Here, 
we  consider  the  effect  of  wave-front  compensation  applied  to  light  incident  in  each  aperture 
of  an  optical  interferometer,  as  shown  in  Figure  11.  Again,  if  all  system  components  as¬ 
sociated  with  aperture  2  are  removed,  the  remaining  components  represent  a  conventional 
AO  system  for  imaging  and  beam-projection  applications. 
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As  was  shown  in  Chapter  II,  the  performance  of  the  interferometric  system  de¬ 
pends  upon  the  aperture-averaged  differential  phase  variance  For  modal  compensation 
through  order  iV,  the  aperture-averaged  residual  differential  phase  variance,  designated  6]^ 
is  calculated  from  Eq.  (74)  as: 

i=2 

N 

+  ^([Aa„i-Aa„if),  (96) 

where  Aao^  and  are  defined  as  the  difference  in  Zernike  phase  coefficients  between 
aperture  1  and  aperture  2  for  the  object  and  conjugate  phase,  respectively: 

Ao-oi  =  Ctoli  ^o2n 

^CLoi  =  ^oli  ^o2i‘ 

In  the  definition  of  Aa^i  and  Aa^^,  and  are  the  object  phase  coefficients  in  each 
aperture,  whereas  cLou  and  ao2i  are  the  conjugate  phase  coefficients  applied  to  DM  1  and 
DM  2,  respectively.  Note  that  the  form  of  in  Eq.  (96)  is  the  same  as  in  Eq.  (86) 
with  <^o  replaced  by  {(j)o\  “  ^02))  cioi  replaced  by  Aa^i,  and  hoi  replaced  by  Aa<,^.  Thus,  by 
use  of  Eq.  (76),  for  A^th  order  atmospheric  compensation  in  an  interferometer,  the  fringe 
Strehl  ratio  is  given  by: 

FSR  =  exp(-i^^).  (99) 

The  goal  of  an  AO  system  is  to  maximize  the  optical  performance,  which,  according  to 
Eq.  (99),  is  equivalent  to  minimizing  5^.  Considering  the  form  of  Eq.  (96),  6]^  is  minimized 
whenever  ([Attoi  “  Aao^]^)  is  minimized  for  all  z.  For  an  ideal  AO  system,  ASof  Aa^i 
for  all  z,  and  therefore  ([Aa^i  —  Aa^i]^)  0.  Thus,  in  an  ideal  AO  system,  reaches  its 

limiting  value  given  by: 

=  /  dpW{^  {[4>ol{R^  -  cl>02{R^f)  -  g  {{^aoif)  .  (100) 
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For  a  non-ideal  AO  system,  the  differential  compensation  will  not  generally  be  equal 
to  the  differential  phase  coefficient  Aa^j,  and  will  be  larger  than  the  limiting  value  given 
in  Eq.  (100). 

If  optimal  modal  wave-front  compensation  is  defined  as  that  correction  corresponding 
to  a  minimization  of  then  the  optimality  criterion  is  to  minimize  ([Aa^i  —  Aa^i]^}  for 
all  i.  Put  another  way,  we  must  find  Aa^i  such  that  ((Aao^  —  Aa^i)^)  is  minimized.  Thus, 
setting  Acioi  =  Aa^i  provides  optimal  wave-front  compensation.  The  optimality  criterion 
of  minimizing  ((Aa^t  —  AaoiY)  corresponds  to  developing  the  minimum-mean-square-error 
(MMSE)  estimator  [43]  of  Aa^^.  The  MMSE  estimate  Aa^i  is  an  estimate  based  upon  some 
measured  data  that  is  related  statistically  to  Aa^i-  The  form  of  the  estimator  and  the 
minimal  value  of  ((Attof  —  Aao^)^)  are  determined  completely  by  the  statistical  relationship 
between  Aa^i  and  the  data  used  to  estimate  that  differential  phase  coefficient. 

As  discussed  in  Chapter  II,  for  a  typical  modal  compensation  scheme,  measured 
coefficients  from  the  AO  beacon  are  applied  directly  in  the  wave-front  correction.  In  our 
notation,  this  is  to  say  that: 

Adf^i  =  dfjii  ““  Ub2i5  (101) 

where  stnd  ab2i  are  the  modal  coefficients  of  order  i  for  the  beacon  phase  in  aperture  1 
and  aperture  2,  respectively.  Due  to  anisoplanatism,  Aa^i  and  Aa^.^  will  only  be  partially 
correlated,  and  the  MMSE  estimate  Aa^i  will  not,  in  general,  be  equal  to  Aa^^  Thus, 
the  use  of  Aa^i  directly  for  wave-front  compensation  will  not  minimize  6^  for  any  order  of 
modal  correction.  This  implies  that  the  wave-front  compensation  can  be  optimized  only 
by  estimating  coefficients  Aa^f  given  measurements  Aa^i,  not  by  using  the  measurements 
directly  as  the  conjugate  phase.  This  fact  has  been  pointed  out  previously  by  Valley  and 
Wandzura  [88]  and  explored  in-depth  by  Whiteley  et  al.  [96]  in  the  case  of  single  aperture 
systems,  a  special  case  of  the  more-general  analysis  presented  here. 

3.2,2  Minimum-mean- square- error  (MMSE)  estimation.  In  this  subsection,  the 
form  of  the  modal  differential  coefficient  estimates  Aa^i  which  optimize  the  wave-front 
compensation  is  developed.  Here,  we  define  the  column  vector  y  whose  components  are 
the  differential  modal  coefficients  associated  with  the  object  of  interest  of  order  1  through 


55 


N: 


o 

<3 
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”  0'o2N 

(102) 


The  MMSE  estimate  of  y,  designated  y,is  the  column  vector  whose  components  are  the 
MMSE  estimates  of  the  differential  modal  coefficients: 


y  = 


Attoi 

Aao2 

Attos 


Aa 


oN 


(103) 


The  vector  x  is  composed  of  the  corresponding  beacon  differential  phase  coefficient  mea¬ 
surements: 


Aaj>i 

^611 

—  ft  621 

Aa;,2 

«612 

“  ft622 

X  = 

Aa63 

^613 

““  ft623 

> 

9 

1 _ _ 

^blN 

“  ft627V 

(104) 


Given  the  measurement  vector  x,  we  can  form  the  MMSE  estimator  of  y  [43].  The  estimate 
y  is  the  conditional  expectation  given  the  measurement  vector  x: 


y  =  (ylx) ,  (105) 

where  the  expectation  is  with  respect  to  the  posterior  probability  density  function  p(y|x). 
In  order  to  construct  the  MMSE  estimator  of  y,  the  joint  PDF  of  x  and  y  must  be  known. 
Since  the  Zernike  coefficients  of  the  object  and  beacon  phase  are  zero-mean  Gaussian,  and 
the  components  of  x  and  y  are  the  difference  of  these  coefficients,  then  x  and  y  are  also 
zero-mean  and  jointly  Gaussian  [36]. 
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For  X  and  y  zero-mean  and  jointly  Gaussian,  it  has  been  shown  [43]  that  the  MMSE 
estimate  of  y  is  given  by: 

y  =  C„,C-ix.  (106) 

where  Cyx  and  Cxx  are  the  covariance  matrices  of  y  with  x  and  x  with  x,  respectively, 
which  for  zero  mean  vectors  y  and  x  are  given  by: 

Cyx  =  (yx^),  (107) 

a,  =  (xx"’),  (108) 

where  is  the  vector  transpose  of  x.  The  form  of  the  estimator  given  in  Eq.  (106)  will 
be  used  throughout  the  remainder  of  this  paper.  Note  that  y  =  x  only  if  CyxC~^  =  L 

The  value  of  given  in  Eq.  (96)  depends  upon  the  choice  of  conjugate  differential 
phase  coefficients  Ahoi.  In  the  case  of  MMSE  estimation,  Aa^i  =  Aa^i.  These  coefficient 
estimates  are  the  elements  of  y,  and  minimize  the  so-called  Bayesian  mean  square  error 
{{Aaoi  —  AaoiY)  [43],  Accordingly,  the  summation  over  i  is  minimized  as  well,  and  may 
be  computed  directly  as: 

-  Aa,i)2)  =  TR  [Cyy  -  GyxC-,^Cxy]  .  (109) 

i=l 

where  TR[  •  ]  designates  the  matrix  trace  operation  Cyy  is  the  covariance  matrix  of  y  with 
y,  and  C^y  is  the  matrix  transpose  of  Cyx* 

If  the  beacon  measurements  are  used  directly  for  wave-front  compensation,  then 
Adoi  =  Aa^i  in  Eq.  (96).  The  differential  phase  coefficients  Aa^i  are  the  elements  of  x,  and 
will  not  generally  minimize  ((Aaoi  “  Aaj,^)^),  nor  will  the  summation  over  i  be  minimized. 
However,  the  summation  can  be  quantified  using  the  same  notation  as  above: 

N 

Y,  {(Aa<,i  -  =  TR  [Cyy  +  C^x  -  2C,,] .  (110) 

t=i 
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Using  the  matrix  expression  for  the  MMSE  estimator  performance  given  in  Eq.  (109),  <5^ 
is  given  by: 

MMSE:  6%  =  l dp  -  TR  .  (Ill) 

If  the  measured  beacon  differential  phase  coefficients  are  used,  then  from  Eq.  (110),  is 
given  by: 

measured:  S%  =  J  dpW{^  {[ct>oiiR^  -  <f>o2{R^?)  -  TR  [2C„,  -  C,,] .  (112) 

The  performance  of  imaging  and  beam-projection  systems  is  quantified  by  calculating 
the  aperture-averaged  residual  phase  variance  which  is  computed  from  Eq.  (Ill)  and 
Eq.  (112)  by  replacing  ((/>oi  -  ^^2)  with  and  by  setting  ao2i,  Cih2i  =  0  in  the  definition  of 
y  and  x.  In  this  case,  the  subscript  1  is  dropped  from  aou  and  and  the  object  and 
beacon  phase  coefficients  are  designated  a^i  and  respectively.  Thus,  and 

In  addition,  the  piston  component  is  excluded  in  y  and  x,  and  all  summations 
begin  with  i  =  2.  The  form  of  the  MMSE  estimator  is  unchanged,  and  the  optimal  modal 
compensation  theory  may  be  applied  directly  to  single-aperture  systems. 

3,3  Spatial  cross- correlation  of  Zernike  coefficients 

In  the  Subsection  3.2.2,  the  MMSE  estimate  of  the  object  differential  phase  coef¬ 
ficients  y  given  beacon  differential  phase  coefficient  measurements  x  depended  upon  the 
covariance  matrices  Cyx  and  Cxx*  Since  —  (fioim  ^o2m)  Q'Hd  =  (n^in  —  ^b2n)j  the 
mn-th  element  of  Cy^  is  given  by: 

^o2m){^bln  “■  ^62n)} 

i^olm^bln)  {^olm^b2n)  (^o2m®i>ln)  “b  {^o2m^b2n)  •  (113) 
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Likewise,  since  Xm  =  {dbim  ”  (^b2m)  and  —  ai,2n)^  the  mn-th  element  of  C^x  is 

given  by: 

~  {(^61m  ”  ®62m)(^61n  “  ®62n)) 

“  i^blm^bln)  (^Mm^62n)  {^b2m^bln)  (^62m^62n)  • 

Similarly,  the  matrix  elements  of  Cyy  correspond  to  the  cross-correlation  of  the  appropriate 
Zernike  coefficients.  Thus,  in  order  to  compute  the  matrix  elements  of  the  MMSE  estimator 
given  in  Eq.  (106),  and  to  measure  compensation  performance  as  specified  in  Eq.  (Ill)  and 
Eq.  (112),  the  cross-correlation  of  the  appropriate  Zernike  coefficients  must  be  computed. 

3.3.1  Analysis  geometry.  Evaluation  of  the  covariance  matrices  Cxx-,  and 
Cyx  involves  the  calculation  of  modal  cross-correlations  which  have  the  following  general 
form: 

{0'uli^v2j)  ~  -®a«itOv2j  ^^2)5  (H^) 

where  Ba^uav2j  represents  the  spatial  cross-correlation  of  auu  with  as  a  function  of 
the  optical  source  locations  '^sv  and  the  aperture  locations  ^^2-  Spatial  cross- 
correlations  of  this  form  address  all  types  of  anisoplanatism  present  in  imaging,  beam- 
projection,  and  interferometric  systems.  Thus,  to  develop  an  expression  for 
consider  a  geometry  with  two  spatially-separated  apertures  and  two  spatially-separated 
sources,  as  depicted  in  Figure  12,  This  geometry  shows  two  apertures  of  radius  iZ,  with 
locations  given  by  r^i  fa2-  The  separation  of  the  apertures  is  designated  s  =  fa2  —  Tai. 
Two  optical  sources  are  located  by  and  r^y.  The  position  vectors  of  the  two  apertures 
and  the  two  sources  share  a  fixed  coordinate  system  specified  by  the  unit  vectors  x,  y,  and 
z. 

A  single  turbulent  atmospheric  layer,  parallel  to  the  x  —  y  plane,  is  indexed  by  the 
subscript  I,  located  by  the  vector  i),  and  has  thickness  Az;.  Light  rays  propagating  from 
the  source  to  each  point  in  the  aperture  will  intersect  the  turbulent  layer  at  height  z/,  thus 
defining  the  projection  of  each  aperture  in  that  layer,  as  is  indicated  in  Figure  12. 
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source  u 


Figure  12.  Analysis  geometry  used  for  computing  the  spatial  cross-correlation  of  Zernike 
expansion  coefficients. 

Consider  the  vectors  Rpi  and  Rp2^  which  locate  points  in  aperture  1  and  aperture  2, 
respectively,  relative  to  the  aperture  centers.  The  projection  of  Rpi  for  source  u  and  the 
projection  of  Rp^  for  source  v  into  the  Ith  layer  are  designated  by  q^u  and  qy2h  respectively. 
The  projected  aperture  vectors  q^u  and  qy^i  are  related  to  Rpi^  Rp2^ 
zi  as  follows: 


9uli  (1 

(116) 

Qv21  =  (1  —  Av2l)Rp2} 

(117) 
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where  the  layer  scaling  factors  Auu  and  Ay2i  are  given  by: 


A^ii 

Avv 


zi  -  (rgi  •  z) 
{r,u-ra\)-z 
zi  -  {ra2  •  z) 
{r>v  -  ra2)  •  z ' 


(118) 

(119) 


An  important  quantity  in  the  development  of  an  expression  for  Ba^uav2j 
jected  separation  of  the  aperture  centers  in  the  Zth  layer,  designated  S/.  In  terms  of  the 
aperture  and  source  position  vectors,  and  the  layer  scaling  factors,  Si  is  given  by: 


Sj  =  s  +  A^2i{r,v  -  Ki)  -  A^u{r>u  -  ^ai)- 


(120) 


Another  important  quantity  is  the  vector  separation  of  the  projected  aperture  vectors  in 
the  Ith  layer,  designated  As/.  As  can  be  seen  in  figure  12,  As^  is  related  to  s;,  and 
q^2i  as  follows: 

=  q^2i  +  Si  —  quu-  (121) 


If  source  u  is  at  an  infinite  range,  and  source  v  remains  at  a  finite  range,  then: 

=  0,  (122) 

Si  =  s  + A^2i{rsv-ra2) 

-[^1  -  (^ai  •>S)](tan0uCos^ui +  tan0„sin<;6^,y),  (123) 

where  0^  is  the  angular  direction  of  source  u  with  respect  to  the  unit  vector  normal  to 
aperture  1,  and  (f)^  is  the  angular  direction  of  the  source  u  with  respect  to  x.  Similarly,  if 
source  v  is  at  an  infinite  range,  and  source  w  is  at  a  finite  range,  then: 

A,21  =  0,  (124) 

Si  =  s  +  [zi~  (fa2  •  z)]  (tan  0y  cos  <j)yX  +  tan  6y  sin  (f)yy) 

-Auu{rsu-rai),  (125) 
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where  9y  is  the  angular  direction  of  source  v  with  respect  to  the  unit  vector  normal  to 
aperture  2,  and  (j)^  is  the  angular  direction  of  source  v  with  respect  to  x.  If  source  u  and 
source  v  both  lie  at  an  infinite  range,  then: 


=  0,  (126) 

A,21  =  0,  (127) 

Si  =  s+[zi-  (fa2  •  z)]  [(tan  Oy  cos  (f>y)x  +  (tan  6y  sin  (j)y)y] 

-[zi  -  (r^ai  •  i)][(tan0^cos(/>^)x  +  {tdm Oy  sin (j>u)y],  (128) 

3.3,2  Spatial  cross-correlation  calculation.  Using  the  analysis  geometry  intro¬ 
duced  in  Section  3.3.1,  we  now  derive  an  expression  for  the  spatial  cross-correlation  of 
the  Zernike  expansion  coefficients.  The  optical  phase  from  source  u  in  aperture  1  and  the 
optical  phase  from  source  v  in  aperture  2,  designated  (l)yi{Rpi)  and  (l>v2{Rp2)j  respectively, 
may  be  expanded  with  respect  to  the  Zernike  polynomials: 


% 

(129) 

<t>v2{Rp2)  =  0‘v2iZj{p2)i 

(130) 

3 


where  a^u  and  ay2j  are  the  Zernike  coefficient  of  the  phase  from  source  u  in  aperture  1 
of  order  i  and  the  Zernike  coefficient  of  the  phase  from  source  v  in  aperture  2  of  order  jf, 
respectively.  These  coefficients  are  obtained  by  projecting  the  phase  onto  the  corresponding 
Zernike  polynomial  as  follows: 


fluii  =  j dpi  (l>^i{Rpi)Zi{pi)W{pi),  (131) 

o-vij  =  j  dp2  (f>v2iRp2)Zj{p2)W{p2),  (132) 

where  W{^  is  the  aperture  weighting  function  defined  in  Eq.  (65). 
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The  spatial  cross-correlation  of  Zernike  coefficient  auu  with  Zernike  coefficient 
is  given  by: 

^auiiavZj  “ 

Substituting  Eq.  (131)  and  Eq.  (132)  into  Eq.  (133)  yields  the  following: 

Ba.ua.,,  =  (/ d^-l  j  dp2  •  (134) 

Since  (f>ui  and  0^,2  sire  the  only  random  quantities  in  Eq.  (134),  the  expectation  operator 
may  be  taken  inside  the  integral  as  follows: 

Ba.ua.„  =  J  dpi  I  dp2  (<^„l(iipi)^„2(i2p2))^i(Pl)^i(p2)W^(Pl)W"(p2).  (135) 

Thus,  to  evaluate  the  integral  of  Eq.  (135),  {<f>ui{Rpi)<Pv2{Rp2))  must  first  be  evaluated. 
Using  a  layered  atmospheric  model,  and  (f>^2  stre  expressed  as: 


(l>uiiRpi)  =  (y) 

(136) 

(i>v2{Rp2)  =  (y)  «/2(^;2J2)> 

(137) 

where  and  Az/j  cire  the  thickness  of  the  atmospheric  layers  indexed  by  and  ^2?  and 
and  n\^{qv2i2)  are  the  index  of  refraction  fluctuations  in  the  layers  indexed  by  li 
and  I2  at  position  and  qv2h^  respectively.  If  the  index  of  refraction  fluctuations  are 
uncorrelated  layer  to  layer  [61,68,81],  that  is, 

{'^'h{9uih)n'iAQv2h))  oc  (138) 

where  6  is  the  Kronecker  delta,  it  follows  that: 

{<i>uiiRpi)<f>v2{Rp2))  =  i'^'i(9uii)n'i{qv2i))  ■  (139) 
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Recognizing  that  the  quantity  {n\{qy,ii)n\{qy2i))  is  the  spatial  autocorrelation  of  the  index 
of  refraction  fluctuations  in  the  /th  layer,  designated  that: 

(140) 

At  this  point,  we  note  that  the  quantity  in  the  square  braces  of  Eq.  (140)  is  the  spatial  au¬ 
tocorrelation  of  the  optical  phase  in  the  /th  layer,  designated  ^^(guih  ?v2i)-  Thus,  Eq.  (140) 
may  be  written: 

{(j>u\{Rpl)^v2{Rp2))  =  ^^B<t>{^uU'i(lv2l)^  (141) 

For  spatially  wide-sense  stationary  statistics,  will  be  a  function  only  of  the  vector 
separation  of  points  specified  by  the  projected  aperture  vectors  quu  and  q^^i  [32].  Prom 
Eq.  (121),  the  vector  separation  is  As/,  and  it  can  be  said  that: 

B^Asl)  =  B^{q,2i  +  Si  -  quuY  (142) 

Hence,  we  conclude  that: 

{<l>ui{RM^v2{Rp2))  =  B^{q^2i  +  Si  —  quii)^  (143) 

/ 

Substituting  Eq.  (143)  into  Eq.  (135),  yields  the  following  expression  for  Ba^-^iav2j- 

=  j  J  *^^2  ]^B^{qy2l  +  Sl-quu)  Zi{pi)Zj{p2)W{px)W{p2).  (144) 

Exchanging  the  order  of  summation  and  integration  and  transforming  the  integration  vari¬ 
ables  to  the  projected  aperture  vectors  quu  and  q„2i  yields  the  following: 


{<i>ul{Rpi)<t>v2{Rp2))^J2  (x)  ^^"Bn’{quU,qv2l) 


Bauiidvij  ^al  ?  )  ^02) 

=  —  Auii)(l  —  Av2i)]~^  J  ^QuU  j  ^qv2i  B^{q^2i  +  Si  —  q^u) 

‘[iJ(l-A„u)J  ni?(l-A„2,)J  U(1-A„2,)J 


.(145) 
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By  employing  Parseval’s  theorem  and  the  Wiener-Khinchin  theorem  [32],  and  the  scaling 
and  shifting  properties  of  Fourier  transforms  [30],  as  in  detailed  in  Appendix  A,  Eq,  (145) 
can  be  written  as: 

■®a„iiOi,2i  (Tsu^  '^sv^  ^a2) 

=  ^  [  dk  W^{k,  z,)  exp(- j27rfc  •  s,)Qi  [j?(l  -  Qj  [i?(l  -  >  (146) 

where  fc  is  a  two-dimensional  spatial  frequency  vector,  W(i,{k^Zi)  is  the  spatial  PSD  of 
the  optical  phase  in  the  turbulent  layer  at  height  Z/,  Si  is  the  separation  of  the  projected 
aperture  centers  given  in  Eq.  (120),  and  Qi{f)  is  the  Fourier  transform  of  Zi{^W{^. 

Integrating  Eq.  (146)  over  the  azimuthal  component  of  fc,  as  detailed  in  Appendix 
B,  yields  the  following  result: 


'®Oui»av2i  ^al  ?  ^02) 

4 

TT 


^  i 

^*^(ni+l)[(l  “  ~~  Ay2l)^] 

poo 

Jo 


X  < 


(147) 


where  D  =  2R\s  the  diameter  of  the  aperture,  Jk(-)  is  the  Bessel  function  of  the  first  kind 
of  order  k,  x  =  irDk,  si  =  |sj|,  and  the  functions  F,  G,  and  H  are  defined  as: 


F  =  [(ni  +  l)(Tij -|-l)]i(-l)i("‘+"^■)2t^"^(*’”*«+*”•^°)^(-l)’”^  (148) 

G  =  (-l)^(’”‘''‘’”’^cos[(mi  +  mj)05, 

+7r/4  {(1  -  ^„,,o)[(-l)‘  -  1]  +  (1  -  «m,o)[(-l)^‘  -  1]}]  ,  (149) 

H  =  (_i)lN<-»"ilcos[(mi-mj)e„ 

+7r/4  {(1  -  ^„,,o)l(-l)‘  -  1]  -  (1  -  5m,o)[(-l)^'  -  1]}]  ,  (150) 
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where  ni,mi  and  nj,mj  are  the  radial  and  azimuthal  orders  of  the  ith  and  jth  Zernike 
polynomial,  respectively,  and  Og,  is  the  angle  between  sj  and  x. 

The  expression  for  in  Eq.  (147)  is  valid  for  any  wide-sense  stationary  and 

isotropic  statistical  model  for  the  turbulence- induced  phase  given  by  W^.  For  this  work, 
the  von  Karman  spectrum  given  in  Eq.  (31)  is  employed.  According  to  Eq.  (36),  W^{k  = 
xfiTD,zi)  is  given  by: 

W4xl^D,z,)  =  (0''’x/3r(8/3)  (x»  +  (151) 

where  Xq  =  ttB/Lq,  Substituting  Eq.  (151)  into  Eq.  (147)  yields  the  following  expression 


^duiiCLv^j  i'^su^  ^al )  ^02) 

roo 

a  f  (x’ +  x|)-“'V, [^x] 

Jo 

X«^(ni+l)[(l  “  ^ull)j^]J{nj-\-l)[0-  ^  . 

X  <  /»oo 

+H  f(x^  +  xg)-“/«J|„,_„,|[^a;] 

Jo 

^  ^  ^ull  )^]*^(n,+l)[(l  - 

According  to  Eq.  (152),  the  total  cross-correlation  can  be  viewed  as  the  sum  of  layer 
contributions  weighted  by  the  factor  Wi  given  by  Eq.  (34).  The  layer  contribution  is 
determined  by  the  relative  projected  aperture  separation  Si/D^  and  by  the  relative  outer 
scale  Lq/D  using  the  parameter  Xq,  Eq.  (152)  serves  throughout  the  remainder  of  this  work 
as  the  primary  means  of  calculating  the  cross-correlation  of  the  Zernike  coefficients  for  a 
given  set  of  aperture  and  source  locations.  This  expression  and  the  accompanying  geometry 
in  Figure  12  can  be  adapted  to  model  all  types  of  spatial  and  temporal  anisoplanatic  effects 
considered  in  this  research. 

3.3,3  Theoretical  validation.  To  verify  the  validity  of  the  spatial  cross-correlation 
expression  developed  here,  Eq.  (152)  can  be  compared  directly  with  published  results  for  a 
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specific  case  of  the  analysis  geometry.  Takato  and  Yamaguchi  [76]  developed  an  expression 
for  the  spatial  correlation  of  Zernike  coefficients  with  finite  outer  scale.  Their  development 
was  for  plane-wave  propagation  from  a  single  source,  with  two  apertures  separated  by  s. 
In  the  notation  of  the  present  analysis,  this  corresponds  to  Ay2i  =  0  and  Si  =  s.  Thus, 
the  spatial  cross-correlation  is  derived  from  Eq.  (152)  as: 


(a,o,)W  = 


poo 

G  /  J(„.+i)(x) J(„.+i)(x) 

<  ''°foo  (153) 

+H  /  +  x2)-“/6 [^X]  J(„;+i)(x) J(„.+i)(x), 

Jo 


where  the  source  and  aperture  subscripts  have  been  dropped  from  the  Zernike  coefficients, 
as  the  sources  are  coincident  and  the  aperture  is  implied  by  the  the  subscripts  i  (aperture 
1)  and  j  (aperture  2).  For  a  given  aperture  separation  direction,  (a^aj)  depends  upon  the 
relative  aperture  separation  s/D  and  the  relative  outer  scale  Lq/D, 


For  separation  in  the  x  direction,  Takato  and  Yamaguchi  [76]  expressed  results  for 
the  autocorrelation  {i  =  j)  of  the  Zernike  modes  in  terms  of  a  correlation  coefficient  defined 
as: 


Paiaj  (^) 


(aiQj)  (g) 

yjiciai)  (0)  (aja^)(0) 


(154) 


Using  the  expression  for  (a^aj)  given  in  Eq.  (153),  the  results  of  reference  [76]  were  repli¬ 
cated  for  piston,  tilt,  and  defocus  correlation  with  Lq/D  =  1, 10,  20,  50, 100.  These  results 
are  shown  in  Figure  13.  The  autocorrelation  plots  in  Figure  13  demonstrate  the  general 
correlative  properties  characteristic  of  the  Zernike  coefficients  for  finite  outer  scale.  For 
all  modes,  the  correlation  coefficient  is  1  for  coincident  apertures,  and  drops  off  as  the 
separation  increases.  The  correlation  coefficient  falls  more  slowly  for  larger  relative  outer 
scale  values.  Of  the  low-order  Zernike  aberrations,  piston  and  tilt  are  the  only  modes  for 
which  the  correlation  coefficient  is  significantly  larger  than  0  for  s/D  >  1. 


Using  Eq.  (153),  the  cross-correlation  (i  ^  j)  of  the  Zernike  coefficients  can  also 
be  quantified.  The  correlation  coefficient  Pa^aj  was  computed  according  to  Eq.  (154)  for 
i  =  2,  J  =4  and  i  =  3,  j  =  5  for  the  same  relative  outer  scale  values  considered  in  Figure  13. 
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Figure  13.  Spatial  autocorrelation  of  low-order  Zernike  coefficients.  Results  are  for  plane- 
wave  propagation  from  a  single  source  with  relative  outer  scale  Lq/D^  and 
relative  aperture  separation  s/D, 

The  results  of  these  calculations  are  shown  in  Figure  14.  These  results  demonstrate  an 
important  characteristic  of  the  cross-correlation  of  Zernike  modes.  While  the  mode  com¬ 
binations  considered  here  (z  =  2,  j  =  4  and  z  =  3,j  =  5)  are  uncorrelated  for  s/Z)  =  0 
(also  true  for  infinite  outer  scale,  see  Table  3),  the  magnitude  of  the  cross-correlation  in¬ 
creases  through  some  maximum  value  as  s/jD  increases,  before  going  to  0  as  s/D  oo. 
This  cross-correlation  effect  will  play  a  significant  role  when  applying  the  optimal  modal 
compensation  theory  developed  earlier,  as  the  enhanced  cross-correlation  between  modes 
will  affect  compensation  performance  when  anisoplanatism  is  present  in  an  AO  system. 
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D 


(a) 


(b) 

Figure  14.  Spatial  cross-correlation  of  low-order  Zernike  coeflScients.  (a)  Cross¬ 
correlation  of  a2  (x-tilt)  with  a4  (defocus)  and  (b)  cross-correlation  of  (y- 
tilt)  with  as  astigmatism.  While  these  modes  are  uncorrelated  for  s/D  =  0, 
the  magnitude  of  the  cross-correlation  increases  through  some  maximum  value 
Bs  s/D  increases. 
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Integrals  of  the  form  contained  in  Eq.  (153)  can  be  evaluated  analytically  [71,76]. 
The  analytic  expression  is  a  series  solution  for  which  convergence  is  may  be  slow.  Thus, 
in  practice,  it  is  often  easier  to  evaluate  these  integrals  numerically,  as  was  done  in  refer¬ 
ence  [76]  and  for  the  results  shown  in  Figure  13  and  Figure  14.  Throughout  the  remainder 
of  this  work,  numerical  integration  is  used  to  evaluate  the  generalized  cross-correlation 
expression  in  Eq.  (152). 

3.3.^  Experimental  validation.  The  spatial  cross-correlation  expression  given 
in  Eq.  (152)  can  also  be  validated  by  comparison  with  experimental  results.  By  invoking 
Taylor’s  frozen-flow  hypothesis  [61]  to  account  for  the  effects  of  bulk  motion  of  the  turbulent 
medium,  the  spatial  cross-correlation  expression  in  Eq.  (152)  may  be  converted  into  an 
expression  for  the  temporal  cross-correlation  of  Zernike  coefficients,  designated  i2o.o^  (r), 
where  r  is  the  time  delay  between  measurements  [37,94].  In  this  case,  Si  =  ViT  in  Eq.  (152) 
for  a  non-zero  flow  velocity  in  the  Ith.  layer  specified  by  the  vector  vi.  If  is  the  zth  Zernike 
coefficient  in  the  aperture  at  time  t,  and  aj  is  the  jth  Zernike  coefficient  in  the  aperture 
at  time  t  +  r,  then  the  temporal  cross-correlation  of  these  coefficients,  designated  Raiaj  is 
given  by: 

W  =(£)’'"  V3r(8/3)2->/»  (II)  F 

\  G  r^ix^  +  Xo2)-“/« (^x)  J(„,+i)(x) 

xy^tz;(<  (155) 

For  comparison  with  theory,  experimental  measurements  of  turbulence- induced  phase 
were  obtained  in  the  laboratory  using  the  optical  apparatus  in  Fig.  15.  In  this  apparatus,  a 
2  cm  diameter  collimated  laser  beam  is  propagated  through  a  channeled  flow  of  turbulence 
generated  by  the  mixing  of  helium  and  nitrogen  gases.  The  velocity  of  the  flow  along  the 
channel  was  approximately  4  m/s.  The  laser  beam  was  centered  approximately  10  cm 
from  the  exit  of  the  turbulence  generator,  where  the  mixing  layer  of  the  two  gases  has 
just  expanded  to  the  width  of  the  channel.  A  Hartmann  wave-front  sensor  was  used  to 
measure  the  local  gradient  of  the  phase  induced  in  the  beam  upon  propagation  through 
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Figure  15.  Optical  bench  used  to  obtain  phase  measurements  for  laser  beam  propagation 
through  a  turbulent  mixing  layer. 


the  turbulent  flow.  A  thorough  description  of  this  apparatus  and  the  experiments  from 
which  the  data  was  obtained  can  be  found  in  reference  [28]  and  reference  [67].  Phase 
measurements  were  obtained  through  a  least-squares  reconstructor  [67],  resulting  in  an 
ensemble  of  255  phase  realizations,  each  represented  as  a  41  X  41  array  at  1  msec  intervals. 
A  sample  of  these  experimental  phase  measurements  is  shown  in  Figure  16. 


From  the  reconstructed  phase  measurements,  the  Zernike  coefficient  ai  at  time  t  is 
calculated  by  projecting  onto  the  Zernike  polynomial  Zi{xm^yn)  ^ts  follows: 


ai{t) 


E41 

m— 1  Z-/n=l 


Vn]  mj  Vn) 


(156) 


where  the  ordered  pair  {xm,yn)  designates  the  physical  location  specified  by  row  m  and 
column  n  in  the  data  array.  Likewise,  the  Zernike  coefficient  at  time  t  +  r  is  calculated 
by: 

aj{t  +  T 


)  = 


r:n=l  Vn'i  ^ 

E'LlEtxZKXm,yn) 


(157) 
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T  =  0 


T  =  1  msec 


Figure  16.  Sample  of  laboratory  measurements  of  turbulence-induced  phase.  Determin¬ 
istic  optical  system  effects  have  been  removed.  Phase  measurements  were 
obtained  on  a  41  x  41  grid  at  1  msec  intervals. 
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Figure  17.  Turbulent  flow  speed  profile  for  the  5-layer  discrete  turbulence  model  employed 
in  calculating  the  theoretical  values  of  the  Zernike  coefficient  cross-correlation. 
The  flow  velocity  in  each  layer  is  taken  to  be  along  the  major  axis  of  the  flow 
channel.  The  flow  speed  in  each  layer  is  normalized  to  the  maximum  speed  in 
the  channel,  t;max- 

Using  these  calculated  Zernike  coefficients,  the  estimated  temporal  cross-correlation  of  ai{t) 
with  aj{t+pTd),  designated  where  p  indexes  the  number  of  time  steps  between 

the  phase  measurements,  was  calculated  as  follows: 

-  255— p 

RaiaiipTd)  =  255 -p  ^  “  l)^d]aj[(3  -  1  +p)T<i].  (158) 

Thus,  RaiajiP'^d)  was  determined  at  =  1  msec  intervals  for  the  reconstructed  phase 
measurements  obtained  in  the  laboratory. 

Eq.  (155)  requires  that  the  turbulent  flow  be  modeled  in  a  finite  number  of  discrete 
layers,  and  that  the  parameters  tf/,  Wi,  and  Lq  be  specified  for  each  of  these  layers.  Based 
on  laboratory  measurements,  the  velocity  in  each  of  5  discrete  layers  was  taken  to  be  in 
direction  of  the  major  axis  of  the  flow  channel.  In  these  analysis  efforts,  the  normalized 
flow  speed  profile  shown  in  Fig.  17  was  employed,  where  the  speed  of  the  turbulent  flow  in 
each  layer  is  divided  by  the  maximum  speed  t^max  in  the  channel.  To  simplify  the  analysis 
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mode  i 

mode  j 

^max  (W®) 

Lo  (m) 

Lo/D 

2  (x-tilt) 

2  (x-tilt) 

4.8 

0.035 

1.75 

3  (y-tilt) 

3  (y-tilt) 

6.4 

0.020 

1.00 

ItlCTHUl 

4  (defocus) 

3.5 

2.50 

gWEgHTil 

8  (coma) 

3.0 

1.75 

3  (y-tilt) 

7  (coma) 

3.5 

0.035 

1.75 

Table  4.  Maximum  flow  speed  (vmax)  and  outer  scale  (Lq)  parameters  selected  to  fit  the 
theoretical  cross-correlation  results  to  the  experimental  data. 

efforts,  an  equal  weighting  for  each  layer  was  assumed;  that  is,  Wi  =  0.2.  The  outer  scale 
Lo  is  also  assumed  to  be  the  same  for  all  layers. 

The  temporal  cross-correlation  of  the  Zernike  coefficients  was  estimated  from  the 
experimental  data  according  to  Eq.  (158)  and  calculated  using  the  theoretical  expression 
given  in  Eq.  (155).  For  each  of  the  results  presented  in  this  section,  the  maximum  flow 
speed,  t^maX)  Q^nd  the  outer  scale  Lq  were  chosen  to  fit  the  theoretical  curves  with  the 
experimental  measurements.  These  quantities  determine  the  free  parameters  Vi  and  Xq 
in  Eq.  (155).  The  values  of  -Vniax  stnd  Lq  chosen  for  each  combination  of  modes  i  and 
j  presented  in  this  section  are  summarized  in  Table  4.  The  velocities  resulting  from  the 
fit  correspond  well  with  the  mean  velocity  of  4  m/s  measured  in  the  laboratory.  The 
magnitude  of  the  outer  scale  values  resulting  from  the  fit  {Lq/D  ^1  —  3)  indicate  that 
finite  outer  scale  effects  are  present  in  the  data  and  that  the  turbulence  cannot  be  modeled 
adequately  with  an  infinite  outer  scale  (Kolmogorov)  spectrum. 

The  experimental  estimates  for  the  temporal  autocorrelation  of  the  orthogonal  com¬ 
ponents  of  tilt  (Zernike  2  and  3)  are  shown  with  error  bars  in  Figure  18(a)-(b).  For  refer¬ 
ence,  the  theoretical  autocorrelation  is  shown  by  the  solid  curves.  Both  the  experimental 
and  theoretical  autocorrelation  have  been  normalized  to  yield  a  correlation  coefficient  in 
the  range  -1  to  1.  From  these  results,  it  is  apparent  that  the  temporal  autocorrelation 
of  the  orthogonal  components  of  tilt  are  well-modeled  by  Eq.  (155)  for  r  <  5  msec.  For 
T  >  5  msec,  it  is  obvious  that  the  frozen-flow  hypothesis  upon  which  Eq.  (155)  is  con¬ 
structed  does  not  fully  account  for  the  temporal  evolution  of  the  turbulence. 

The  experimental  and  theoretical  cross-correlation  of  tilt  with  defocus  (Zernike  4)  and 
orthogonal  components  of  coma  (Zernike  7  and  8)  is  shown  in  Figure  19.  Only  the  cross- 
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(b) 

Figure  18.  Comparison  of  theoretical  and  experimental  temporal  autocorrelation  results 
for  tilt  components,  (a)  Temporal  autocorrelation  of  the  x-tilt  (Zernike  2)  and 
(b)  temporal  autocorrelation  of  the  y-tilt  (Zernike  3)  for  the  reconstructed 
phase  measurements  studied.  Experimental  results  are  shown  with  error  bars, 
and  theoretical  results  are  indicated  by  the  solid  curves. 
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Figure  19.  Comparison  of  theoretical  and  experimental  temporal  cross-correlation  results. 

Experimental  results  are  shown  with  error  bars,  and  theoretical  results  are 
indicated  by  the  solid  curves. 
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correlation  of  x-tilt  with  defocus  is  shown  as  the  cross-correlation  of  y-tilt  with  defocus 
is  zero  for  all  values  of  r.  Similarly,  x-tilt  is  uncorrelated  with  Zernike  7  and  y-tilt  is 
uncorrelated  with  Zernike  8. 

As  is  shown  by  the  theoretical  cross  correlation  curves,  tilt  is  uncorrelated  with 
defocus  at  r  “  0,  as  is  conventionally  reported  [54]  (also  see  Table  3).  As  r  increases, 
however,  the  magnitude  of  the  cross-correlation  increases  through  some  maximum  value 
before  going  to  zero  for  larger  values  of  r.  The  experimental  cross-correlation  of  tilt  with 
defocus  follows  the  same  trend,  and  is  well-modeled  by  theory  over  a  range  of  r  values.  The 
theoretical  curves  in  Figure  19  also  demonstrate  that  the  tilt  components  are  correlated 
with  orthogonal  coma  components  at  r  =  0  as  the  conventional  analysis  indicates  [54]. 
Tilt  remains  correlated  with  coma  for  r  >  0  before  becoming  uncorrelated  at  larger  values 
of  r. 

While  the  theoretical  cross-correlation  results  shown  in  Figure  19  generally  agree 
with  the  experiment  considering  measurement  error,  there  is  a  noted  difference  between 
theory  and  experiment  at  smaller  values  of  r.  As  this  was  not  the  case  for  the  autocor¬ 
relation  results  in  Figure  18,  validity  of  the  frozen  flow  hypothesis  cannot  be  questioned 
in  explaining  this  disagreement.  Instead,  differences  between  the  experimental  and  theo¬ 
retical  results  are  due  to  statistical  anisotropy,  which  has  been  noted  to  exist  in  mixing 
layers  [28,29]. 

S.Ji.  Summary 

This  chapter  extends  the  analysis  of  optical  system  performance  degraded  by  turbulence- 
induced  anisoplanatism.  By  considering  the  form  of  the  aperture-averaged  residual  differ¬ 
ential  phase  variance  for  A^th-order  modal  compensation,  system  performance  was 
shown  to  be  suboptimal  when  beacon  modal  phase  measurements  are  used  directly  for 
wave-front  correction.  To  achieve  optimal  (but  not  necessarily  ideal)  performance,  the 
modal  correction  must  be  estimated  from  the  beacon  phase  coefficients,  given  knowledge 
of  the  correlation  of  object  and  beacon  phase.  The  optimality  criterion  determined  by 
the  form  of  corresponds  to  MMSE  estimation,  where  the  mean  value  of  the  squared 
difference  between  the  object  and  conjugate  phase  coefficients  is  minimized  for  all  modes. 
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The  formalism  and  notation  associated  with  MMSE  estimation  was  introduced,  and  the 
performance  of  the  estimator  was  quantified  in  terms  of  covariance  matrices  for  the  beacon 
and  object  phase  coefficients. 

In  order  to  calculate  the  MMSE  estimator  matrix,  the  spatial  cross-correlation  of 
the  object  and  beacon  phase  coefficients  must  be  quantified.  To  model  the  full  range  of 
anisoplanatic  effects  in  imaging,  beam-projection,  and  interferometric  systems,  the  cross¬ 
correlation  of  Zernike  coefficients  in  spatially-separated  apertures  from  spatially-separated 
optical  sources  must  be  quantified.  An  analysis  geometry  was  introduced  to  model  these 
effects,  and  the  cross-correlation  calculation  was  performed  for  wide-sense  stationary  and 
isotropic  phase  statistics  with  finite  outer  scale.  For  a  layered  atmospheric  model,  the 
total  cross-correlation  is  a  weighted  sum  of  the  cross-correlation  in  each  layer.  For  each 
layer,  the  cross-correlation  depends  upon  the  ratio  5//J9,  where  s/  is  the  projected  aperture 
separation  in  the  Zth  layer,  and  D  is  the  diameter  of  the  apertures.  The  cross-correlation 
also  depends  upon  the  relative  outer  scale  Lq/D, 

The  expression  resulting  from  the  cross-correlation  calculation  may  be  tailored  for 
analysis  of  a  system-specific  geometry.  To  test  the  validity  of  the  cross-correlation  analysis, 
the  results  of  this  calculation  were  compared  to  published  theoretical  results  [76],  and  were 
found  to  be  in  exact  agreement.  The  expression  for  the  cross-correlation  also  predicts 
that  while  most  low-order  Zernike  modes  are  uncorrelated  at  zero  spatial  separation,  the 
magnitude  of  the  cross-correlation  increases  through  some  maximum  value  before  going  to 
zero  as  the  spatial  separation  becomes  large  relative  to  the  aperture  diameter.  To  test  this 
prediction  and  to  further  confirm  the  validity  of  the  cross-correlation  expression,  laboratory 
measurements  of  turbulence-induced  phase  were  analyzed.  Experimental  measurements 
of  the  temporal  autocorrelation  where  found  to  be  in  close  agreement  with  theoretical 
predictions  when  frozen-flow  assumptions  were  valid.  The  cross-correlation  enhancement 
effect  predicted  by  the  theory  was  also  confirmed  by  the  experimental  results. 
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IV.  Ground-Based  Adaptive- Optical  Imaging 

4 . 1  Introduction 

In  this  chapter,  the  optimal  modal  compensation  theory  developed  in  Chapter  III  is 
applied  to  ground-based  adaptive-optical  imaging  systems.  Adaptive  optics  improve  the 
resolution  of  ground-based  imaging  systems  by  conjugating  the  turbulence-induced  phase 
incident  in  the  system  aperture.  When  anisoplanatism  is  present  in  the  AO  system,  the 
imaging  performance  is  degraded,  since  the  incident  and  conjugate  phase  are  only  partially 
correlated.  By  accounting  for  the  statistics  of  the  turbulence-induced  phase,  part  of  the 
performance  degradation  may  be  recovered  through  optimal  compensation. 

The  application  of  the  optimal  modal  compensation  theory  is  discussed  in  Section  4.2. 
The  minimum-mean-square-error  (MMSE)  estimator  for  the  optimal  modal  phase  coeffi¬ 
cients  is  discussed,  and  the  performance  of  the  estimator  is  given  in  terms  of  the  statistics 
of  the  object  and  beacon  phase.  The  spatial  correlation  of  the  object  and  beacon  phase 
is  given  by  a  special  case  of  the  general  correlation  expression  in  Eq.  (152).  Thus,  both 
angular  and  focal  anisoplanatism  effects  may  be  accounted  for  in  the  optimal  estimator. 

The  performance  of  natural-guide-star  (NGS)  and  laser-guide-star  (LGS)  compensa¬ 
tion  with  the  optimal  estimator  is  discussed  in  Section  4.3  and  Section  4.4.  For  NGS  com¬ 
pensation,  the  performance  degradation  resulting  from  angular  anisoplanatism  is  shown  to 
be  reduced  when  the  MMSE  estimator  of  object  phase  coefficients  is  employed.  Thus,  the 
NGS  beacon  angular  offset  can  be  increased  without  degrading  imaging  performance.  For 
the  LGS  beacon  aided  by  a  NGS  beacon  for  tilt  compensation,  angular  and  focal  anisopla- 
natic  effects  are  reduced  by  the  optimal  estimator.  Thus,  the  LGS  beacon  alititude  can  be 
decreased  without  degrading  imaging  performance. 

Section  4.5  discusses  how  tilt  compensation  performance  may  be  enhanced  by  in¬ 
corporating  higher-order  modes  in  a  MMSE  estimator  for  wave- front  tilt.  For  NGS  tilt 
compensation,  the  maximum-permissible  angular  separation  between  the  object  and  NGS 
is  significantly  increased.  Recovery  of  object  tilt  components  using  only  LGS  beacon  mea¬ 
surements  is  discussed.  In  this  scheme,  while  only  higher-order  modal  measurements  can 
be  obtained  from  the  LGS,  the  cross-correlation  of  the  higher-order  modes  with  tilt  com- 
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ponents  can  be  utilized  in  the  MMSE  estimator  to  provide  an  estimate  of  wave-front  tilt. 
The  performance  of  the  tilt  estimator  is  quantified,  demonstrating  its  ability  to  enhance 
imaging  performance  as  well  as  practical  limits  to  its  application.  Section  4.6  provides  a 
brief  summary  of  the  chapter. 


Application  of  theory  to  ground-based  adaptive  optics 

The  optimal  modal  compensation  theory  discussed  in  Section  3.2  may  be  applied 
directly  to  ground-based  imaging  with  an  AO  system.  Recall  that  a  lower  bound  on  the 
Strehl  ratio  of  an  imaging  system  is  given  in  Eq.  (90)  as  SR  >  exp(—  The  aperture- 
averaged  residual  phase  variance  is  computed  using  Eq.  (86): 

m)  -  E  (4> + E  ((»..  - «..)“)  •  (159) 

i=2  i=2 

The  Zernike  coefficients  of  the  conjugate  phase  are  represented  by  and  are  based  upon 
measurements  of  the  beacon  phase  Suppose  the  modal  coefficients  are  measured 
from  the  beacon  phase  incident  in  the  system  aperture: 


dbi  = 


I  dpWi^  mR^  Zii^- 


(160) 


Recognizing  that  €%  has  the  same  mathematical  form  as  6%  in  Eq.  (96),  and  considering  the 
optimal  modal  compensation  theory  developed  in  Chapter  III,  is  minimized  by  setting 
^oi  =  where  hoi  is  the  MMSE  estimate  of  aoi  given  af,i.  For  imaging  applications,  we  are 
concerned  only  with  piston-removed  phase.  Thus,  the  vector  y  corresponding  to  Eq.  (102) 
is  composed  of  object  phase  coefficients  of  order  2  through  AT,  whereas  x  is  composed  of 
the  corresponding  Zernike  phase  coefficients  measured  from  the  NGS  or  LGS  beacon: 


^02 

<^b3 

y  = 

X  = 

(^oN 

O^bN 

(161) 
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The  components  of  y  are  the  MMSE  estimates  of  the  object  phase: 

^o2 
^o3 

y  tto4 

^oN 

According  to  Eq.  (106),  y  is  related  to  x  by: 

y  =  C'vxC'->,  (163) 

where  Cy^  is  the  covariance  matrix  of  y  with  x  and  Cxx  is  the  covariance  matrix  of  x  with 
X,  as  in  Eq.  (107)  and  Eq.  (108),  respectively.  When  the  optimal  coefficients  are  used  in 
modal  compensation,  then  according  to  Eq.  (109),  is  given  by: 

MMSE:  6%  =  I  dp  W{^  {<f>l{Rp))  -  TR  [CyxC^^Cxy]  .  (164) 

If  the  measured  beacon  phase  coefficients  are  used  directly  to  form  the  conjugate  phase 
(no  optimal  estimation),  then  according  to  Eq.  (110),  is  given  by: 

measured:  4  =  j  dp  W{p)  {<t>l{Rp))  -  TR  [2Cyx  -  C^x]  ■  (165) 

The  matrix  elements  of  Cyy^  C^x  and  Cy^  are  given  by  evaluating  the  spatial  cross¬ 
correlation  expression  of  Eq.  (152)  for  a  geometric  specification  consistent  with  a  ground- 
based  imaging  system.  This  special  case  of  the  generalized  analysis  geometry  in  Figure  12 
is  shown  in  Figure  20,  where  Vai  =  ^a2  =  For  the  single-aperture  case,  the  aperture 
separation  5  =  0,  and  the  projected  aperture  separation  Si  is  given  by  Eq.  (120). 

4-3  Natural- guide- star  compensation 

To  quantify  NGS  compensation  performance,  must  be  calculated  according  to 
Eq.  (164)  for  MMSE-estimates  of  the  object  phase,  and  by  Eq.  (165)  for  employing  beacon 


(162) 
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source  v 


Figure  20.  Special  case  of  analysis  geometry  for  ground-based  imaging  systems. 


phase  measurements  directly.  The  first  term  of  is  the  aperture-averaged  piston-removed 
phase  variance,  which  for  finite  outer  scale  Lq  is  given  by: 


jApW{^{<t>l{R^)  =  (10  V3r(8/3)  2-^/3  (1^) 

^  ^  (a;) 

[10  7o  ^ 


(166) 


where  Xq  =  irD/Lo,  as  in  Eq.  (152).  Since  each  term  of  Eq.  (164)  and  Eq.  (165)  is  propor¬ 
tional  to  {D/ro)^^^  [see  Eq.  (152)],  compensation  performance  may  be  quantified  by  the 
normalized  aperture-averaged  residual  phase  variance,  £^(E)/ro)“®/®.  When  normalized  in 
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this  fashion,  will  depend  only  on  the  parameter  Xq,  which  is  fully  specified  by  the  rela¬ 
tive  outer  scale  Lq/D.  For  the  results  which  follow,  aim  aperture  diameter  is  assumed, 
and  angular  offsets  are  normalized  to  the  isoplanatic  angle  Oq  defined  in  Eq.  (91). 

Figure  21(a)  shows  the  value  of  £%{D/ro)~^^^  for  compensation  using  an  NGS  with 
angular  offset  6b/Oo  =  0  to  10  for  N  =  6,10,15,21,  and  Lq/D  =  100.  Similarly,  Fig¬ 
ure  22(a)  shows  the  value  of  e]^{D for  compensation  using  an  NGS  with  angular 
offset  Ob/O^  =  0  to  2  for  iV  =  28,36,45,55,  and  =  50.  The  dashed  horizontal  lines 

in  Figure  21(a)  and  Figure  22(a)  correspond  to  =  1  rad^  for  each  value  of  the  relative 
aperture  diameter  D/tq  considered.  For  each  relative  aperture  diameter,  the  normalized 
error  must  be  less  than  the  value  indicated  for  Strehl  ratio  performance  greater  than  0.37. 
Performance  results  for  the  MMSE  estimator  are  shown  by  the  solid  curves  for  each  correc¬ 
tion  order.  Performance  results  for  coefficients  measured  directly  from  the  NGS  are  shown 
by  the  dashed  curve.  As  is  illustrated  in  Figure  21(a)  and  Figure  22(a),  the  normalized 
error  from  the  MMSE  estimator  is  always  less  than  or  equal  to  the  normalized  error  using 
phase  coefficients  directly  as  the  conjugate  phase. 

Figure  21(b)  shows  the  Strehl  ratio  of  an  AO  system  with  D/ro  =  5  corresponding  to 
the  results  in  Figure  21(a).  The  MMSE  estimator  is  shown  to  yield  imaging  performance 
greater  than  or  equal  to  the  imaging  performance  when  beacon  phase  coefficients  are  used 
directly  in  the  compensation.  For  the  orders  of  correction  considered  in  Figure  21(b),  the 
difference  between  the  imaging  performance  of  the  MMSE  compensation  and  the  imaging 
performance  using  beacon  measurements  directly  is  appreciable  for  9b /Oq  >  2. 

Figure  22(b)  shows  the  Strehl  ratio  performance  of  an  AO  system  with  =  10 
corresponding  to  the  results  in  Figure  22(a).  Again,  the  MMSE  estimator  yields  imaging 
performance  greater  than  or  equal  to  the  imaging  performance  when  beacon  phase  coeffi¬ 
cients  are  used  directly  in  the  AO  compensation.  For  the  orders  of  correction  considered 
in  Figure  22(b),  the  difference  between  the  imaging  performance  of  the  MMSE  compensa¬ 
tion  and  the  imaging  performance  using  beacon  measurements  directly  is  appreciable  for 
9b  19^  >  0.5.  For  N  =  55,  the  MMSE  estimator  can  provide  equivalent  imaging  perfor¬ 
mance  at  beacon  offset  angles  approximately  10%  larger  than  when  beacon  measurements 
are  used  directly  in  the  compensation. 
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(b) 

Figure  21.  Comparison  of  lower-order  NGS  optimal  compensation  performance  with  con¬ 
ventional  compensation  performance,  (a)  Value  of  for  the  NGS 

beacon  offset  0j,/0o  using  the  MMSE  estimator  (solid  hne)  and  using  bea¬ 
con  measurements  directly  (dashed  line).  Horizontal  lines  indicate  levels  at 
which  =  1  rad^  for  each  value  of  D/rQ.  (b)  Corresponding  Strehl  ratio 
performance  with  relative  aperture  diameter  £)/ro  =  5. 
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Figure  22.  Comparison  of  higher-order  NGS  optimal  compensation  performance  with 
conventional  compensation  performance,  (a)  Value  of  for  the 

NGS  beacon  offset  Ob/Oo  using  the  MMSE  estimator  (solid  line)  and  using 
beacon  measurements  directly  (dashed  line).  Horizontal  lines  indicate  levels 
at  which  =  1  rad^  for  each  value  oi  D/tq.  (b)  Corresponding  Strehl  ratio 
performance  with  relative  aperture  diameter  D/tq  =  10. 
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Figure  23.  Schematic  diagram  of  an  AO  system  with  LGS  compensation  aided  by  an 
off-axis  NGS  for  tilt  compensation. 

Natural- guide- star  aided  laser- guide- star  compensation 

When  investigating  LGS  compensation,  we  must  consider  the  effect  of  using  an  off- 
axis  NGS  for  tilt  compensation  along  with  an  on-axis  LGS  for  higher-order  modes.  An  AO 
system  in  which  the  LGS  is  aided  by  an  NGS  for  tilt  compensation  is  shown  in  Figure  23. 
The  LGS  beacon  is  formed  at  some  finite  altitude  and  the  NGS  has  an  angular  offset 
0b  with  respect  to  the  optical  axis.  In  this  configuration,  the  beacon  phase  coefficient 
measurements  are  decorrelated  from  the  object  phase  coefficients  a^i  due  to  both  angular 
anisoplanatism  (NGS  effect)  and  focal  anisoplanatism  (LGS  effect).  The  form  of  the  beacon 
measurement  vector  does  not  change,  but  the  tilt  components  and  are  associated 
with  the  NGS. 
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One  conceptual  difficulty  that  arises  in  modeling  LGS  effects  with  the  geometry 
shown  in  Figure  20  is  the  notion  of  specifying  a  fixed  position  vector  for  the  LGS.  When 
an  on-axis  laser  beam  propagates  upward  through  the  atmosphere  to  form  the  beacon,  a 
random  deflection  is  induced  by  the  turbulence.  Thus,  the  angle  associated  with  the  vector 
location  of  the  beacon  is  a  random  quantity.  Tyler  [85]  has  estimated  the  variance  of  this 
random  deflection,  cr^,  to  be  given  by: 

=  0,182(A/Dy(D/ro)^/\  (167) 

For  D  =  1  m,  A  =  500nm,  and  D/vq  <  10,  ae  is  less  than  2  /xrad.  Since  the  angular 
offset  between  the  LGS  and  tilt  NGS  considered  in  these  results  is  more  than  an  order  of 
magnitude  larger  than  cr^,  the  random  deflection  of  the  beam  is  negligible,  and  the  LGS  is 
considered  to  be  fixed  on-axis  at  the  specified  altitude. 

Figure  24(a)  shows  the  value  of  ej^(D/ro)~^^^  for  compensation  using  an  LGS  at 
altitude  10  to  30  km  for  N  =  6,10,15,21,  with  NGS  offset  0i,/0o  =  5  and  Lq/D  =  100. 
Similarly,  Figure  25(a)  shows  the  value  of  for  compensation  using  an  LGS  at 

altitude  10  to  30  km  for  N  =  28, 36, 45, 55,  with  NGS  offset  9b/0o  =  5  and  Lq/D  =  50.  As 
was  the  case  in  the  NGS  compensation  results,  the  dashed  horizontal  lines  in  Figure  24(a) 
and  Figure  25(a)  correspond  to  =  1  rad^  for  each  value  of  the  relative  aperture  diameter 
D/tq  considered.  Performance  results  obtained  with  the  MMSE  estimator  are  shown  by 
the  solid  curves.  The  dashed  curves  indicate  the  compensation  performance  when  beacon 
phase  measurements  are  used  directly  for  the  conjugate  phase.  These  results  show  that  the 
normalized  error  associated  with  the  MMSE  estimator  is  always  less  than  the  normalized 
error  associated  with  using  phase  coefficients  directly  at  all  LGS  altitudes,  and  for  all  orders 
of  correction.  The  results  in  Figure  24(a)  show  that  the  relative  error  reduction  with  the 
MMSE  estimator  is  small  for  the  lower  orders  of  correction  {N  =  6,10,15,21).  For  the 
higher  orders  of  correction  {N  =  28,36,45,55)  shown  in  Figure  25(a),  the  relative  error 
reduction  with  the  MMSE  estimator  is  greater,  with  enhanced  compensation  performance 
evident  at  all  LGS  altitudes. 
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(b) 

Figure  24.  Comparison  of  lower-order  LGS  optimal  compensation  performance  with  con¬ 
ventional  compensation  performance,  (a)  Value  of  for  an  on- 

axis  LGS  at  altitude  10  to  30  km  and  a  tilt-compensation  NGS  offset  by 
0^1  Oq  =  5  using  the  MMSE  estimator  (solid  line)  and  using  beacon  mea¬ 
surements  directly  (dashed  line).  Horizontal  lines  indicate  levels  at  which 
=  1  rad^  for  each  value  of  (b)  Corresponding  Strehl  ratio  perfor¬ 

mance  with  relative  aperture  diameter  D/ro  =  5. 
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Figure  25.  Comparison  of  higher-order  LGS  optimal  compensation  performance  with  con¬ 
ventional  compensation  performance,  (a)  Value  of  for  an  on- 

axis  LGS  at  altitude  10  to  30  km  and  a  tilt-compensation  NGS  offset  by 
=  5  using  the  MMSE  estimator  (solid  line)  and  using  beacon  mea¬ 
surements  directly  (dashed  line).  Horizontal  lines  indicate  levels  at  which 
=  1  rad^  for  each  value  of  D/ro-  (b)  Corresponding  Strehl  ratio  perfor¬ 
mance  with  relative  aperture  diameter  D/ro  =  10. 
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Figure  24(b)  shows  the  Strehl  ratio  of  the  LGS/NGS  AO  system  with  D/r^  =  5 
corresponding  to  the  results  in  Figure  24(a).  Figure  25(b)  shows  the  Strehl  ratio  of  the 
LGS/NGS  AO  system  with  D/ro  =  10  corresponding  to  the  results  in  Figure  25(a).  For 
all  orders  of  correction  and  at  all  LGS  altitudes,  the  MMSE  compensation  yields  better 
imaging  performance  than  when  beacon  measurements  are  used  directly  for  compensation. 
For  the  higher  correction  orders  [N  —  28,36,45,55),  the  MMSE  estimator  reduces  the 
effect  of  focal  anisoplanatism,  recovering  some  of  the  imaging  performance  lost  by  using 
a  lower-altitude  LGS.  For  N  =  55,  the  MMSE  estimator  can  provide  equivalent  imaging 
performance  at  LGS  altitudes  approximately  20%  lower  than  when  beacon  measurements 
are  used  directly  in  the  compensation. 

4^5  Tilt  compensation  enhancement 

For  ground-based  AO  imaging  applications,  accurate  compensation  for  wave- front  tilt 
is  critical,  since  the  variance  of  the  tilt  phase  is  large  compared  to  higher-order  modes.  As 
stated  in  Eq.  (88),  the  aperture- averaged  piston-removed  phase  variance  for  Kolmogorov 
turbulence  (infinite  outer  scale)  is  1.0299(Z)/ro)^/^.  From  the  covariance  matrix  shown  in 
Table  3,  the  sum  of  the  variance  of  the  tilt  components  (a^2)  +  (^03)  =  0,896(i)/ro)^/^. 
Thus,  the  wave- front  tilt  components  account  for  nearly  87%  of  the  total  aperture-averaged 
piston-removed  phase  variance.  For  finite  outer  scale  (von  Karman)  turbulence,  the  per¬ 
centage  of  the  piston-removed  phase  variance  represented  by  tilt  will  be  smaller  than  for 
Kolmogorov  turbulence,  but  will  still  be  large  compared  to  all  higher-order  modes. 

When  anisoplanatism  is  present  in  an  AO  system,  tilt  compensation  performance 
will  be  degraded.  Using  the  optimal  modal  compensation  theory,  however,  some  of  the 
performance  degradation  may  be  recovered  when  the  partial  correlation  of  the  tilt  com¬ 
ponents  is  properly  accounted  for  in  the  MMSE  estimator  for  tilt.  The  MMSE  estimator 
also  provides  a  framework  for  further  enhancing  tilt  compensation  performance.  If  the 
object  phase  coefficient  vector  y  contains  only  tilt  components  ao2  and  a^a,  and  the  beacon 
measurement  vector  contains  tilt  plus  higher-order  modes  a62,a63?G^64  •  •  •  then  any 

cross-correlation  between  tilt  and  the  higher-order  modes  will  further  reduce  the  value  of 
given  in  Eq.  (164)  over  the  value  of  when  only  tilt  measurements  are  employed. 
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As  discussed  in  Subsection  3.3.3,  while  tilt  is  uncorrelated  with  most  higher-order 
modes  for  concurrent  apertures,  the  magnitude  of  cross-correlation  of  tilt  and  many  of  the 
higher-order  modes  is  enhanced  for  moderate  values  of  the  relative  aperture  separation 
s/D.  A  similar  effect  occurs  when  there  is  an  angular  separation  between  the  object  and 
the  beacon  in  an  AO  system.  This  angular  offset  produces  an  effective  aperture  separation 
in  the  turbulence  layers  above  the  system  aperture,  as  depicted  in  Figure  20.  The  resulting 
cross-correlation  between  the  object  wave-front  tilt  and  the  beacon  wave-front  higher-order 
modes  can  be  used  to  enhance  tilt  compensation  performance  from  an  off-axis  beacon,  as 
is  illustrated  in  the  remainder  of  this  section. 

4^5 A  Single-layer  turbulence  model  The  analysis  associated  with  the  cross¬ 
correlation  between  tilt  and  higher-order  modes  for  an  angular  separation  between  the 
object  and  beacon  can  be  simplified  and  the  results  can  be  generalized  if  a  single-layer  model 
is  invoked  for  the  turbulence  [93].  For  a  single- layer  model,  the  turbulence  is  regarded  as 
being  concentrated  at  a  particular  altitude  above  the  aperture.  Since  the  turbulence  is 
concentrated  within  a  single  layer,  the  weight  for  that  layer  is  wi  =  1.  Note  that  many 
of  the  turbulence  profiles  presented  in  Table  1  are  highly  concentrated  in  a  single  layer. 
For  instance,  for  the  HV-21  profile,  the  layer  at  altitude  200  m  has  an  associated  weight 
of  0.8902. 

The  aperture  and  source  geometry  for  a  single-layer  turbulence  model  is  shown  in 
Figure  26.  The  turbulence  layer  is  located  at  an  altitude  zi  above  the  aperture  of  diameter 
D.  The  beacon  is  located  at  an  altitude  hi,  above  the  aperture,  and  is  displaced  by  the  angle 
01,  from  the  direction  of  the  object  (the  optical  axis  of  the  system).  Consider  identifying 
the  distant  object  in  Figure  26  with  source  u  in  Figure  20,  and  identifying  the  beacon 
with  source  v.  Accordingly,  since  the  object  is  at  an  infinite  range,  the  layer  scaling  factor 
associated  with  the  object  is  A^u  =  0.  From  Eq.  (119),  the  layer  scaling  factor  associated 
with  the  beacon  is  Ay2i  =  zi/ hi,.  For  Oi,  in  the  x  direction,  the  projected  aperture  separation 
in  the  layer  is  given  by  Si  =  Zi  tan^e,.  Thus,  from  Eq.  (152),  the  cross-correlation  of  on-axis 
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Figure  26.  Geometry  for  a  single-layer  turbulence  model.  Projected  aperture  separation 
is  given  by  $i  =  2;tan05  for  turbulence  layer  altitude  zi  and  beacon  angle  6^. 


object  phase  coefficients  with  the  off-axis  beacon  coefficients  is  given  by: 

G  r^ix^  + 

Jo 

X-^(ni+l)[a;]J(n,+l)[(l  -  Zl/hb)x] 

rOO 

+H  / 

Jo 

X-/(ni+l)MJ'(n,+l)[(l  -  Zl/hb)x], 

where  F,  G,  and  H  are  evaluated  according  to  Eq.  (148)-(150)  with  0^^  =  0.  As  can  be  seen 
from  Eq.  (168),  for  a  given  mode  combination  the  cross-correlation  is  parameterized  by 
ratio  of  the  layer  altitude  to  the  beacon  altitude,  zi/hi,,  the  projected  separation  relative  to 
the  aperture  diameter,  (z/tan^^)/!),  and  Xq  =  ivD/Lo,  The  parameter  Zi/h^^  characterizes 
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focal  anisoplanatism,  (^/tan^^)/!)  characterizes  angular  anisoplanatism,  and  xq  accounts 
for  the  finite  outer  scale  of  the  turbulence. 

4-5,2  Higher-order  modal  cross-correlation.  The  cross-correlation  between  the 
object  wave-front  tilt  and  modal  beacon  measurements  is  shown  in  Figure  27.  In  these 
figures,  Ooi  {i  —  2,3)  represent  the  tilt  components  of  the  object  wave-front,  whereas 
{j  =  2,3,...,  15)  represent  the  modal  measurements  from  the  beacon  for  tilt  through 
Zernike  15.  The  cross-correlation  is  expressed  by  the  correlation  coefficient  Paoiabj  computed 
as  in  Eq.  154.  The  autocorrelation  results  (i  =  j)  shown  in  Figure  27  illustrate  the 
expected  behavior  of  maximum  correlation  at  =  0  and  decreasing  correlation  with 
increasing  beacon  angular  offset.  As  indicated  in  these  figures,  at  =  0  all  higher-order 
modes  considered  here  are  uncorrelated  with  tilt  except  the  orthogonal  components  of 
coma  {j  =  7, 8).  The  magnitude  of  the  correlation  coefficient  for  all  other  modes  increases 
as  (z/tan^b)/!?  increases,  passing  through  some  maximum  value  before  going  to  0  at 
large  values  of  (z/  tan^^)/!).  The  enhanced  cross-correlation  between  tilt  and  higher-order 
modes  for  (z/  tan^^)/!?  >  0  indicates  that  the  higher-order  measurements  from  an  off-axis 
beacon  may  be  used  in  the  MMSE  estimator  for  the  object  tilt  components  to  improve 
tilt-compensation  performance. 

4-5,3  Natural- guide- star  tilt  compensation.  To  quantify  the  tilt-compensation 
performance  enhancement  obtained  by  including  higher-order  modes  in  the  MMSE  tilt 
estimator,  the  value  of  the  aperture-averaged  residual  phase  error  after  tilt  correction  is 
computed.  This  quantity  is  designated  (subscript  T  indicates  tilt  compensation  only), 
and  is  equal  to  for  =  3  calculated  by  use  of  Eq.  (164).  Thus,  is  proportional 
to  {D/roY^^  [see  Eq.  (166)  and  Eq.  (168)],  and  compensation  performance  may  be  quan¬ 
tified  by  the  normalized  aperture-averaged  residual  phase  variance  after  tilt  correction, 
eUD/ro)-^/\ 

Figure  28(a)  shows  the  value  of  e\{D  for  an  off-axis  NGS  employed  as  the 

AO  beacon  with  Zi/hi  =  0.  These  results  are  obtained  for  the  single-layer  turbulence 
model  with  Lq/I?  =  100  and  beacon  angular  offsets  (z/tan06)/-D  in  the  range  from  0 
to  0.5.  The  dashed  curve  shows  the  compensation  performance  when  the  NGS  beacon 
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Figure  27.  Cross-correlation  of  tilt  with  higher-order  modes  for  angular  beacon  offset 
with  single-layer  turbulence  model.  Results  are  shown  for  the  object  wave- 
front  tilt  components  with  beacon  measurements  through  Zernike  15  with 
=  and  Lq/D  =  100  (uncorrelated  modes  omitted). 
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(a) 


(b) 

Figure  28.  NGS  tilt  compensation  performance  with  MMSE  estimator  including  higher- 
order  modes.  (a)  Value  of  e\,{D for  the  NGS  beacon  offset 
(z;  tan^i,)//)  using  the  MMSE  estimator  for  tilt  (solid  lines)  and  using  tilt 
measurements  directly  (dashed  line),  (b)  Corresponding  Strehl  ratio  perfor¬ 
mance  with  relative  aperture  diameter  D/tq  =  4.  Results  shown  are  for 
L^jD  =  100,  Zifhi  =  0. 
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tilt  measurements  are  used  directly  for  wave-front  tilt  correction  (no  MMSE  estimation 
employed).  The  solid  curves  show  the  effect  of  using  the  MMSE  tilt  estimator  with  the 
modal  beacon  measurements  indicated.  These  results  illustrate  that  if  only  tilt  measure¬ 
ments  are  used  in  the  MMSE  estimator,  optimal  compensation  will  reduce  the  value  of 
at  larger  values  of  the  NGS  offset.  When  higher-order  modes  are  included  in 
the  tilt  estimator,  e\{D may  be  reduced  significantly,  especially  at  larger  values  of 
(z;  tan0e,)/i5. 

Figure  28(b)  shows  the  Strehl  ratio  performance  corresponding  to  each  curve  in  Fig¬ 
ure  28  for  D/ro  =  4.  For  including  modal  measurements  through  Zernike  15  in  the  MMSE 
tilt  estimator,  the  Strehl  ratio  remains  within  1%  of  the  limiting  performance  (at  zero  an¬ 
gular  offset)  out  to  {zi  tdin0i,)/D  =  0.15.  For  using  only  tilt  measurements  in  the  estimator, 
the  corresponding  limiting  performance  angle  is  {zi  ta,n  Of,) /D  =  0.05  (only  1/3  as  great). 
For  the  largest  NGS  offset  angle  considered,  using  the  beacon  measurements  directly  for 
tilt  compensation  results  in  SR  <  0.02.  By  incorporating  beacon  modal  measurements 
through  Zernike  15  in  the  tilt  estimator,  the  Strehl  ratio  is  increased  by  nearly  5  times  to 
SR  0.1.  As  the  results  in  Figure  28  illustrate,  for  a  fixed  NGS  offset,  the  tilt  compensa¬ 
tion  performance  may  be  enhanced  by  including  higher-order  modes  in  the  tilt  estimator. 
Correspondingly,  an  equivalent  level  of  tilt  compensation  performance  may  be  obtained  at 
larger  beacon  offset  angles  when  an  optimal  compensation  is  applied  by  the  AO  system. 

4-5.4  Laser- guide- star  tilt  estimation.  As  was  discussed  in  Section  2.3,  the  LGS 
signal  does  not  contain  tilt  phase  information  since  the  tilt  components  of  the  LGS  wave- 
front  are  naturally  compensated  by  the  reciprocal  propagation  from  the  beacon  to  the 
receiving  aperture.  This  imposes  a  limitation  on  the  use  of  the  LGS  in  adaptive  optics.  In 
practice,  therefore,  the  LGS  is  used  for  higher-order  compensation,  with  tilt  compensation 
provided  by  means  of  a  nearby  NGS. 

Within  the  past  few  years,  research  has  been  conducted  to  develop  methods  of  mea¬ 
suring  tilt  from  the  LGS.  Belen’kii  has  proposed  several  techniques  [5-7]  in  which  LGS  tilt 
is  sensed  by  observing  the  upward  propagation  of  the  laser  beam  with  auxiliary  telescopes 
separated  from  the  transmitter  in  transverse  directions.  Belen’kii  has  argued  that  the  tilt 
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contribution  from  the  downward  path  to  the  auxiliary  telescope  may  be  reduced  by  av¬ 
eraging  over  the  angular  extent  of  the  beacon,  thus  singling  out  only  the  tilt  induced  in 
the  upward  path.  This  tilt  measurement  corresponds  to  the  tilt  induced  in  light  incident 
from  the  astronomical  object  of  interest.  Belen’kii  has  also  proposed  a  modification  to  this 
LGS  tilt  sensing  method  [8],  which  involves  averaging  the  LGS  image  over  the  position  of 
the  auxiliary  telescope  in  an  array,  and  averaging  over  subaperture  positions  within  the 
auxiliary  telescopes.  Ragazzoni  et  al.  [57-59]  have  also  proposed  techniques  for  LGS  tilt 
sensing  that  involve  observing  or  transmitting  the  beacon  from  auxiliary  telescopes  lo¬ 
cated  in  transverse  directions.  Recently,  Lukin  [50]  has  reported  on  a  method  for  optimal 
tilt  correction  using  tilt  measurements  from  monostatic  and  bistatic  LGS  adaptive  optics 
configurations. 

The  method  for  LGS  tilt  estimation  proposed  in  this  research  is  substantially  dif¬ 
ferent  from  the  methods  of  Belen’kii,  Ragazzoni  et  al.,  and  Lukin.  Since  the  methods 
proposed  previously  attempt  to  measure  tilt  directly  from  the  LGS  signal,  efforts  must  be 
made  to  overcome  the  beam  reciprocity  effect,  thus  leading  to  the  bistatic  configurations 
proposed.  These  configurations  in  turn  induce  tilt  measurement  inaccuracies  which  must 
be  overcome  through  novel  methods  (angular,  multiple  aperture,  and  subaperture  aver¬ 
aging  in  Belen’kii’s  method;  concurrent  NGS  tilt  sensing  in  Ragazzoni’s  method;  optimal 
estimation  in  Lukin’s  method).  In  the  method  proposed  here,  we  ignore  the  tilt  compo¬ 
nents  associated  with  the  LGS  completely,  measuring  only  higher-order  modal  aberrations, 
which  are  assumed  to  be  unaffected  by  beam  reciprocity.  When  a  LGS  is  positioned  off-axis, 
the  correlation  between  the  on-axis  tilt  components  from  an  astronomical  object  and  the 
higher-order  LGS  modes  is  sufficiently  strong  to  provide  for  tilt  sensing  using  the  MMSE 
tilt  estimator.  Thus,  wave-front  tilt  components  are  not  measured  directly  from  the  LGS 
signal,  but  instead,  the  higher-order  measurements  are  used  exclusively  to  estimate  the 
wave- front  tilt. 

According  to  Eq.  164,  the  performance  of  the  MMSE  tilt  estimator  depends  upon 
the  cross-correlation  of  the  object  tilt  components  with  higher-order  beacon  measurements 
(Cyx)  as  well  as  the  cross-correlation  of  the  beacon  modes  (Cg-x).  Considering  the  cross¬ 
correlation  results  in  Figure  27,  the  performance  of  a  tilt  estimator  based  only  on  beacon 
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Figure  29.  OfF-axis  deployment  geometry  for  LGS  tilt  estimation,  (a)  Single  beacon  at 
angle  0^  in  x  direction,  (b)  Two  beacon  configuration  along  the  x  direction, 
both  at  angle  9^,,  (c)  Four  LGS  configuration  with  beacons  offset  along  x  and 
y  axes,  each  at  angle  0^, 


measurements  beginning  with  Zernike  4  will  be  highly  dependent  upon  the  geometry  of 
the  LGS  deployment.  Consider  the  three  off-axis  LGS  deployment  configurations  shown  in 
Figure  29.  These  diagrams  show  (a)  the  LGS  geometry  for  the  case  of  a  single  beacon  offset 
angularly  in  the  x  direction,  (b)  two  beacons  offset  symmetrically  along  the  x  direction, 
and  (c)  four  beacons  offset  along  the  x  and  y  directions,  each  at  angle  0^. 

In  order  to  quantify  performance  for  each  configuration  shown  in  Figure  29,  we  will 
consider  all  wave-front  distortion  higher  than  tilt  to  be  perfectly  compensated  (presumably 
by  an  on-axis  LGS  not  shown).  Therefore,  the  residual  phase  variance  after  tilt  compen¬ 
sation  is  due  only  to  errors  in  estimating  the  actual  tilt  components  associated  with  the 
on-axis  object  of  interest.  In  this  case,  is  given  by: 


e^  =  TR  [Cyy 


Cy.C-^C,y]  , 


(169) 
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Figure  30.  Tilt  compensation  performance  from  a  single  LGS.  Value  of  £^{D/ro)~^^^  has 
been  quantified  using  the  beacon  modal  measurements  indicated  in  the  MMSE 
tilt  estimator.  Results  shown  are  for  Lq/D  =  100  and  =  0.01. 

where  y  represents  the  on-axis  tilt  components  associated  with  the  object  and  x  repre¬ 
sents  higher-order  measurements  associated  with  each  LGS.  Again,  since  all  matrix  ele¬ 
ments  in  Eq.  (169)  are  proportional  to  (D/ro)®/®,  the  normalized  residual  phase  variance 
e‘^{D/ro)~^^^  is  considered  in  the  results  which  follow. 

The  tilt  estimation  performance  for  the  single  LGS  configuration  is  shown  in  Fig¬ 
ure  30.  The  value  of  was  computed  for  (zjtan0(,)/D  in  the  range  0  to  1. 

Each  curve  corresponds  to  the  residual  phase  variance  after  applying  the  tilt  correction  es¬ 
timated  from  the  higher-order  beacon  measurements  indicated.  For  all  mode  combinations 
considered,  £^{D/ro)~^^^  attains  its  smallest  value  at  an  LGS  angle  (z;  tan^i,)/!)  ~  0.5, 
the  exact  value  depending  upon  the  particular  modal  measurement  combination  employed. 
The  general  effect  of  including  higher-order  modes  through  Zernike  15  is  to  reduce  the  value 
of  the  normalized  residual  phase  variance.  Notice  that  inclusion  of  the  coma  components 


99 


(Zernike  7  and  8)  reduces  the  error  at  (z/tanfl^)/!)  =  0,  since  these  modes  are  correlated 
with  the  tilt  components  at  zero  angular  offset  (see  Figure  27). 

Tilt  estimation  performance  results  are  shown  in  Figure  31(a)-(b)  for  the  multiple- 
LGS  beacon  configurations  depicted  in  Figure  29(b)-(c).  Each  curve  represents  the  per¬ 
formance  obtained  by  incorporating  the  modes  indicated  into  the  measurement  vector, 
beginning  with  Zernike  4  (defocus)  and  continuing  through  Zernike  15.  These  results 
illustrate  that  increasing  the  number  of  LGS  beacons  will  enhance  the  performance  of 
the  tilt  estimator,  as  the  higher-order  modes  from  each  beacon  are  correlated  with  the 
on- axis  tilt  components.  The  smallest  value  of  e\{D Jr occurs  for  beacon  offsets  of 
(;2;/ tan0fe)/D  =  0.4  —  0.6,  depending  upon  the  modes  included  in  the  measurement  vector 

X. 

The  magnitude  of  the  turbulence  outer  scale  also  has  a  significant  effect  on  tilt 
compensation  performance  using  the  LGS  measurements.  Recall  that  the  cross-correlation 
of  Zernike  coefficients  in  Eq.  (168)  depends  upon  the  parameter  Xq  =  itD/Lq,  Thus,  the 
value  of  will  vary  depending  upon  the  value  of  the  relative  outer  scale  LofD,  For  the 
four-LGS  beacon  configuration,  estimating  tilt  by  measuring  Zernike  modes  4  to  15  from 
each  beacon,  was  quantified  for  Lq/D  =  10,20,50,100.  These  results  are  shown  in 
Figure  32(a).  As  indicated  in  this  figure,  decreases  at  all  values  of  {zita,n6i^)/D  with 
decreasing  values  of  Lq/D.  Figure  32(b)  shows  the  Strehl  ratio  performance  corresponding 
to  each  curve  in  Figure  32(a)  for  the  relative  aperture  diameter  D/tq  =  4.  For  Lq/D  =  100, 
the  Strehl  ratio  is  no  greater  than  0.02  for  any  beacon  offset  angle.  For  smaller  values 
of  Lq/D^  increased  Strehl  ratio  performance  is  noted  at  all  beacon  offset  angles.  For 
Lq/D  =  10,  the  Strehl  ratio  is  over  0.3  for  {zitsin9t)/D  ^  0.6.  Considering  that  the 
majority  of  experimental  studies  [2]  have  measured  outer  scale  values  in  the  range  from 
5  m  to  100  m,  smaller  values  of  Lq/D  may  often  be  the  case  in  practice.  Thus,  the  LGS 
tilt  estimation  method  considered  here  may  offer  a  significant  gain  in  imaging  performance 
when  tilt  information  is  otherwise  unavailable. 

As  was  mentioned  in  Subsection  4.5.1,  the  parameter  zi/ht  characterizes  the  effect 
of  focal  anisoplanatism  on  the  Zernike  coefficient  cross-correlation  in  Eq.  (168).  This 
parameter  specifies  the  altitude  of  the  effective  turbulence  layer  relative  to  the  altitude 
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(a) 


(b) 

Figure  31.  Tilt  compensation  performance  for  multiple  LGS  configurations  for  beacon 
measurements  indicated  on  each  curve,  (a)  for  two-LGS  con¬ 
figuration,  and  (b)  for  four-LGS  beacon  configuration.  Results 

shown  are  for  Lq/D  =  100  and  zi/h^  =  0.01. 
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D 


(a) 


(b) 

Figure  32.  Effect  of  finite  outer  scale  on  LGS  tilt  compensation  performance,  (a)  Nor¬ 
malized  residual  phase  variance  and  (b)  Strehl  ratio  for  beacon 

angular  offset  (z/  tBJi0t)/D,  Results  shown  for  four-LGS  beacon  configuration 
using  Zernike  modes  4  to  15  in  the  measurement  vector  for  tilt  estimation 
with  Zi/hi  =  0.01. 
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Figure  33.  Effect  of  focal  anisoplanatism  on  LGS  tilt  compensation  performance.  Results 
are  shown  for  the  four-LGS  beacon  configuration  using  Zernike  modes  4  to  15 
in  the  measurement  vector  for  tilt  estimation  with  Lq/jD  =  100. 

of  the  beacon.  For  the  results  presented  thus  far,  zi/hf,  =  0.01  has  been  employed  in  all 
calculations.  This  value  of  2/ //if,  was  chosen  since  effective  layer  altitudes  are  typically  on 
the  order  of  0.1  km  to  1  km,  whereas  LGS  beacon  altitudes  are  on  the  order  of  10  km 
to  100  km.  If  the  effective  layer  altitude  is  higher,  or  the  beacon  altitude  is  lower,  then 
focal  anisoplanatism  will  effect  the  performance  of  LGS  tilt  estimation.  This  effect  was 
quantified  by  calculating  the  value  of  for  zifh^)  —  0.01,0.1,0.2,0.5  using  the  four-LGS 
beacon  configuration,  estimating  tilt  by  measuring  Zernike  modes  4  to  15  from  each  beacon. 
These  results  are  shown  in  Figure  33,  illustrating  that  focal  anisoplanatism  degrades  the 
performance  of  LGS  tilt  estimation,  since  increases  for  increasing  values  of  zi/hi^.  This 
performance  degradation  comes  about  due  to  the  reduced  correlation  of  the  on-axis  object 
tilt  with  higher-order  beacon  modal  measurements  in  the  MMSE  tilt  estimator. 
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Figure  34.  Beacon  angular  offset  for  peak  LGS  tilt  estimation  performance.  Each  curve 
shows  the  turbulence  layer  altitude  zi  and  beacon  offset  6t  consistent  with 
{ziia,n0i,)/D  =  0.6  for  each  value  of  aperture  diameter  D  indicated. 

For  tilt  estimation  using  the  four-LGS  beacon  configuration  with  Zernike  modes  4 
to  15  in  the  measurement  vector,  tilt  estimation  performance  is  maximized  at  a  non-zero 
beacon  offset  angle,  as  is  illustrated  in  Figure  31(b)  and  Figure  32.  This  angle  is  determined 
by  the  combined  effect  of  the  enhanced  cross-correlation  between  the  on-axis  object  tilt 
and  the  off-axis  higher-order  beacon  measurements.  As  is  illustrated  in  Figure  32,  for  the 
range  of  relative  outer  scale  values  considered,  the  tilt  compensation  performance  reaches 
its  peak  near  (^/tan^^)//?  =  0.6.  Thus,  to  achieve  the  greatest  possible  tilt  compensation 
performance  with  the  method  discussed  here,  the  LGS  beacon  offset  angle  should  be  such 
that: 

9t  =  tan"^  ("~)  • 

For  a  given  aperture  diameter  D,  Eq.  (170)  establishes  a  relationship  between  Zi  and  9t,. 
This  relationship  is  shown  graphically  in  Figure  34  for  several  values  of  D.  Figure  34 
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indicates  that  for  =  1  km,  should  be  on  the  order  of  10  arcsec  when  D  =  0.1  m. 
When  £)  =  10  m,  however,  0^  should  be  approximately  2,000  arcsec  (over  0.5  deg)  for 
zi  =  \  km.  Intermediate  aperture  diameters  will  require  to  be  approximately  100  arcsec 
to  1,000  arcsec.  Since  1  arcsec  =  4.85  /^rad,  and  the  isoplanatic  angle  0o  is  nominally  7//rad 
for  the  HV-21  profile,  then  for  peak  LGS  tilt  compensation  performance  may  be  as 
much  as  several  hundred  times 

4,6  Summary 

In  this  chapter,  the  optimal  modal  compensation  theory  developed  in  Chapter  III 
was  applied  to  AO  compensation  for  ground-based  imaging  systems.  The  turbulence- 
induced  phase  from  an  NGS  or  LGS  beacon  is  only  partially  correlated  with  the  phase 
from  the  object  of  interest  when  the  beacon  is  offset  from  the  optical  axis  or  at  finite 
altitude.  Thus,  the  wave-front  compensation  performance  obtained  by  applying  beacon 
phase  measurements  directly  in  an  AO  system  is  sub-optimal.  Instead,  to  achieve  optimal 
compensation  performance,  the  object  phase  must  be  estimated  using  the  MMSE  estimator 
given  modal  beacon  measurements. 

For  the  off-axis  NGS,  the  MMSE  estimator  resulted  in  compensation  performance 
that  was  equal  to  or  better  than  the  compensation  performance  using  beacon  measurements 
directly  for  all  orders  of  correction  and  at  all  beacon  offset  angles.  For  the  correction  orders 
=  6, 10, 15, 21  with  aim  aperture,  the  MMSE  estimator  provides  enhanced  performance 
for  beacon  offsets  greater  than  20o«  For  N  =  28, 36, 45, 55  with  aim  aperture,  the  MMSE 
estimator  provides  enhanced  performance  for  beacon  offsets  greater  than  0.5  0o*  For  the 
highest-order  correction  studied,  the  MMSE  estimator  can  provide  equivalent  imaging 
performance  at  NGS  offset  angles  10%  larger  than  when  NGS  measurements  are  used 
directly  for  compensation. 

For  the  case  of  an  on-axis  LGS  with  an  off-axis  NGS  for  tilt  compensation,  the  MMSE 
estimator  resulted  in  enhanced  compensation  performance  over  using  beacon  measurements 
directly  throughout  the  entire  range  of  LGS  altitudes  considered.  The  relative  performance 
enhancement  was  greater  for  the  higher  correction  orders.  The  MMSE  estimator  partially 
reduces  the  effect  of  focal  anisoplanatism.  For  the  highest-order  correction  studied,  the 
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MMSE  estimator  can  provide  equivalent  imaging  performance  at  altitudes  20%  lower  than 
when  beacon  measurements  are  using  directly  for  wave-front  compensation. 

The  MMSE  estimator  was  also  employed  in  order  to  enhance  tilt  compensation  from 
an  off-axis  NGS.  Using  the  MMSE  estimator  for  object  tilt  components,  higher-order  modes 
can  be  included  in  the  beacon  measurement  vector.  For  an  off-axis  NGS,  many  of  these 
higher-order  modes  are  correlated  with  the  on-axis  tilt  components.  Thus,  including  these 
modes  in  the  tilt  estimator  results  in  increased  Strehl  ratio  performance  over  using  tilt 
measurements  only.  For  tilt  compensation  from  an  off-axis  NGS  using  beacon  Zernike 
modes  2  to  15,  it  was  shown  that  the  maximum  beacon  offset  angle  for  nearly-limiting 
performance  is  3  times  larger  than  when  only  tilt  measurements  are  used.  Equivalently, 
for  a  given  beacon  offset,  the  Strehl  ratio  is  increased  by  incorporating  higher-order  modes. 

While  tilt  components  cannot  be  measured  conventionally  from  the  LGS  wave-front 
due  to  beam  reciprocity  effects,  the  MMSE  tilt  estimator  permits  the  tilt  components 
to  be  estimated  using  only  higher-order  measurements  from  the  LGS  beacon.  Since  the 
higher-order  modes  from  an  off-axis  LGS  beacon  are  correlated  with  on-axis  tilt,  the  off- 
axis  higher-order  beacon  measurements  enable  estimation  of  the  on-axis  tilt  components. 
The  performance  of  such  an  estimator  depends  upon  the  geometric  configuration  of  the 
LGS  beacons,  which  beacon  modes  are  used,  the  relative  turbulence  outer  scale  Lq/D,  and 
the  extent  of  focal  anisoplanatism  degradation.  It  was  shown  that  for  a  constellation  of 
four  LGS  beacons  offset  orthogonally  from  the  optical  axis,  using  measurements  of  Zernike 
modes  4  to  15  from  each  beacon,  the  Strehl  ratio  can  be  over  0.3  for  D/tq  =  4  and  Lq/D  = 
10.  Performance  is  degraded  for  larger  values  of  Lq/D.  The  tilt  estimation  performance  is 
degraded  further  when  focal  anisoplanatism  reduces  the  cross-correlation  between  tilt  and 
the  higher-order  modal  measurements.  Peak  LGS  tilt  estimation  performance  occurs  at 
an  angle  which  depends  upon  the  diameter  of  the  aperture  and  the  height  of  the  effective 
turbulence  layer.  This  peak  angle  was  shown  to  be  significantly  larger  than  the  isoplanatic 
angle  Oq  for  typical  AO  systems. 
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V,  Airborne  Laser  Beam-Steering  Adaptive  Optics 

5.1  Introduction 

For  a  system  such  as  Airborne  Laser  (ABL)  [11,26],  which  attempts  to  project  a 
laser  beam  through  the  atmosphere  onto  a  distant  object  over  an  extended  period  of 
time,  turbulence-induced  phase  aberrations  degrade  system  performance.  Random  angular 
beam  deflections  result  in  reduced  energy  deposition  on  the  surface  of  the  target  object, 
whereas  spatial  dispersion  reduces  the  incident  beam  energy  density.  Therefore  in  ABL 
target  engagements,  beam  steering  should  account  for  the  effects  of  turbulence-induced 
wave-front  tilt.  This  chapter  describes  the  application  of  the  optimal  modal  compensation 
theory  developed  in  Chapter  III  to  beam-steering  adaptive  optics  for  the  ABL. 

Due  to  the  finite  temporal  response  of  adaptive-optical  systems,  steering  commands 
must  be  based  upon  time-delayed  wave-front  measurements.  The  dynamics  of  the  ABL 
target  engagement  scenario  and  turbulence-induced  anisoplanatism  cause  the  time-delayed 
wave-front  measurements  to  be  decorrelated  from  the  required  wave-front  tilt  correction. 
This  decorrelation  leads  to  suboptimal  beam-steering  performance  when  time-delayed 
wave- front  measurements  are  used  directly  for  compensation.  In  Section  5.2,  optimal 
modal  compensation  theory  is  applied  to  the  ABL  target  engagement  scenario  for  opti¬ 
mal  compensation  by  an  AO  system,  illustrating  the  necessity  to  quantify  the  temporal 
cross-correlation  of  Zernike  coefficients  for  a  dynamic  system.  Section  5.3  discusses  how 
aperture  and  source  motion  effects  are  modeled  in  the  analysis  geometry  associated  with 
the  optimal  compensation  theory.  ABL  target  engagement  results  are  presented  in  Section 
5.4,  demonstrating  that  beam-steering  performance  is  enhanced  by  incorporating  tilt  plus 
higher-order  time-delayed  measurements  in  the  MMSE  estimator  for  wave-front  tilt.  The 
chapter  is  summarized  in  Section  5.5. 

5.2  Application  of  theory  to  Airborne  Laser  beam- steering  adaptive  optics 

As  was  discussed  in  Chapter  II,  the  performance  of  a  beam-projection  system  such 
as  the  Airborne  Laser  is  also  quantified  by  a  Strehl  ratio  which  depends  upon  the  aperture- 
averaged  phase  variance.  When  a  conjugate  phase  is  applied  to  the  wave-front  prior  to 
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beam  propagation,  the  residual  phase  variance  determines  the  performance.  In  the  case  of 
adaptive-optical  beam  steering  control,  a  conjugate  tilt  phase  is  applied  to  the  outgoing 
beam,  and  the  performance  of  the  system  depends  upon  the  aperture-averaged  residual 
phase  variance  after  tilt  correction,  designated  and  calculated  as: 

4  =  1  dpw{^  {[URp.t)  -  URp,t)f)  •  (171) 

The  corresponding  lower  bound  on  Strehl  ratio  is  given  by: 

SR  >  exp(— g:^).  (172) 

The  tilt  correction  phase  t)  may  be  expanded  in  terms  of  the  orthogonal  components 

of  tilt,  ao2(0  aosit)  (Noll’s  ordering  [54]).  By  employing  the  orthonormality  property 
of  the  Zernike  polynomials  on  the  unit  circle,  may  be  expressed  as: 

3  3 

dpW{^  {<t>l{Rp,  t))  -  ^  ^  (Ki(i)  -  ,  (173) 

t=2  i=2 

where  is  the  coefficient  of  the  expansion  of  (f>o{Rp^t)  with  respect  to  the  Zernike 

polynomials: 

aoiit)  =  I  dpW{^  URp,t)Zi{^-  (174) 

The  expression  for  has  the  same  form  as  6%  in  Chapter  III,  and  therefore  the  opti¬ 
mal  modal  compensation  theory  may  also  be  applied  to  AO  compensation  in  ABL  beam 
steering. 

In  ABL  target  engagements,  the  conjugate  tilt  coefficients  ao2{t)  and  hosit)  are  de¬ 
termined  by  measurements  from  a  beacon  source.  Due  to  the  finite  temporal  response  of 
the  AO  system,  a  time  delay  exists  between  the  beacon  measurements  at  time  t  —  r  and 
the  tilt  correction  applied  to  the  outgoing  laser  beam  at  time  t.  The  measured  beacon  tilt 
coefficients,  designated  at,2(t  —  r)  and  -  r),  are  determined  by  the  projection  of  the 
beacon  phase  ^^(JRp,  t  —  r)  onto  the  appropriate  Zernike  polynomials: 

G'biit  -t)  =  J  dp  W {^(j)f,{Rp^  t  -  r)Zi{^,  for  i  =  2, 3.  (175) 


108 


Due  to  aperture  and  source  motion  within  the  delay  period  r,  — r)  will  be  decorrelated 
from  aoi{t)  [94].  Thus,  if  the  beacon  measurements  are  used  directly  for  the  conjugate  tilt, 
that  is  aoi{t)  =  then  the  compensation  will  be  suboptimal.  Optimal  performance 

may  be  obtained  by  estimating  the  coefficients  aoi{t)  given  measurements  abi{t  —  r)  using 
an  MMSE  estimator  for  the  object  wave-front  tilt  components. 


The  notation  associated  with  MMSE  tilt  estimation  follows  that  of  the  MMSE  es¬ 
timator  introduced  in  Subsection  3.2.2.  For  the  beam-steering  application,  however,  only 
the  wave-front  tilt  components  need  to  be  estimated.  The  tilt  acquired  by  the  laser  beam 
as  it  propagates  to  the  target  object  is  described  by  the  tilt  vector  y  whose  elements  are 
the  orthogonal  components  of  the  wave-front  tilt  at  time  t: 


y  = 


ao2{t) 

O-osit) 


(176) 


The  beacon  measurement  vector  x  may  contain  tilt  plus  higher-order  modal  components 
measured  at  time  t  —  t: 


062  (<  -  t) 

063  (<  -  r) 


X  = 


064  (<  -  t) 


(177) 


OftAr(t  -  t) 


The  tilt  estimate  vector  y  applied  at  time  t  based  on  measurements  at  time  t  —  r  will 
contain  two  elements  corresponding  to  the  MMSE  estimates  of  the  tilt  acquired  by  the 
propagating  laser  beam: 


y  = 


002  (*,  r) 


Oo3(i,T) 


(178) 


According  to  Eq.  (106),  y  is  related  to  y  and  x  by: 


yx^xx 


(179) 
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When  the  MMSE  estimator  is  used  in  the  tilt  compensation,  according  to  Eq.  (109),  is 
given  by: 

4  =  1  dp  W{^  {4ARP,  t))  -  TR  [Cy,C;^C,y] .  (180) 

The  MMSE  estimate  of  the  phase  coefficients  y  given  time-delayed  phase  coefficient 
measurements  x  depends  upon  the  covariance  matrices  Cyx  and  Cxx-  Since  =  Oo(m+i)(t) 
and  Xn  =  ai,(„+i)(t  -  r),  the  mn-th  element  of  Cyx  is  given  by: 

[^!/®]mn  ~  (®'o(m+l)(^)®6(n+l)(i  “  '^))  =  m  +  l)a'(n-|-l)  (r),  (181) 

where  i?o(^+i)a(n.hi)('^)  the  temporal  cross-correlation  of  ao(m+i)(0 

Likewise,  since  Xm  =  <i6(m+i)(t  “  t)  and  Xn  =  at(n+i)(t  —  t),  the  mn-th  element  of  is 

given  by: 

“  (^6(m+l)(t  “  ““  ■^<i(Tn-|-i>0(n+i) (^)*  (182) 

Similarly,  the  matrix  elements  of  Cyy  and  C^y  correspond  to  the  temporal  cross-correlation 
of  the  appropriate  Zernike  coefficients. 

5.5  Temporal  cross-correlation  for  aperture  and  source  motion 

The  temporal  cross-correlation  of  Zernike  coefficients  when  aperture  and  source  mo¬ 
tion  are  present  has  been  developed  previously  by  Whiteley  et  al.  [94].  The  analysis 
geometry  and  cross-correlation  results  are  summarized  in  this  section.  Figure  35  shows  an 
aperture  of  radius  R  whose  position  vector  at  times  separated  by  a  period  r  is  designated 
by  fa(t  -  r)  and  fa(t)-  An  optical  source  is  located  at  times  separated  by  a  period  r  by  the 
vectors  fg{t  —  r)  and  fs(t),  respectively.  Atmospheric  layers  are  indexed  by  the  subscript  Z, 
located  at  height  Zi^  have  thickness  and  a  wind  velocity  Vi.  At  a  given  time,  light  rays 
propagating  from  the  source  to  each  point  in  the  aperture  will  intersect  each  layer,  thus 
defining  the  temporal  projections  of  the  aperture  in  that  layer,  as  is  indicated  in  Figure  35. 
Due  to  the  aperture  and  source  motion  during  the  time-delay  period  r,  the  analysis  geom¬ 
etry  in  Figure  35  corresponds  exactly  with  the  analysis  geometry  for  the  interferometric 
system  in  Figure  12.  This  correspondence  indicates  that  the  anisoplanatism  introduced  by 


no 


source  at 
time  t-z 


Figure  35.  Analysis  geometry  used  for  computing  the  temporal  cross-correlation  of 
Zernike  coefficients  for  aperture  and  source  motion. 


time-delay  in  a  dynamic  system  is  a  generalized  form  of  spatial  anisoplanatism.  Therefore, 
the  spatial  cross-correlation  analysis  associated  with  Figure  12  may  be  applied  directly  to 
Figure  35  with  proper  accounting  for  the  aperture  and  source  motion. 

Consider  the  case  where  the  aperture  and  source  velocities,  designated  Va  and  are 
constant  with  respect  to  the  period  r.  Accordingly,  the  position  vectors  of  the  aperture 
and  source  at  time  t  —  r  and  time  t  are  related  as  follows: 


ra{t)  =  rait-T)  +  VaT, 

(183) 

rs{t)  =  f,{t-T)+V,T, 

(184) 
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Comparing  the  analysis  geometry  in  Figure  35  to  Figure  12,  we  recognize  that  fa{t  —  r) 
corresponds  to  aperture  location  fa{t)  corresponds  to  aperture  location  fa2^  rs{t  —  r) 
corresponds  to  source  location  f\,i,  and  f^(t)  corresponds  to  source  location  7^32.  The 
effective  projected  aperture  separation  in  the  /th  layer  as  a  function  of  r,  designated 
is  obtained  using  Eq.  (120),  with  ra2j  ^ai  and  ^51  related  by: 

ra2  =  Kl  +  VaT,  (185) 

rs2  =  r,i  +  v,T.  (186) 


By  employing  Taylor’s  frozen-flow  hypothesis  to  account  for  the  temporal  effects  of  the 
wind  velocity  v/,  s/  '(r)  is  given  by  [94]: 

Si  '(r)  =  VaT  +  [Ai{t)  -  A((0)](f;  -  fa)  +  Mt){v,  -  Va)T  -  ViT,  (187) 


where  the  layer  aperture  scaling  factor  A/(r)  is  given  by: 

,  .  {Zl  -fg-Z)-  {Vg  •  Z)t 

‘  {f,  -fa)-Z  +  [(■?,  -  Vg)  •  z]t  ’ 


(188) 


and  fa,  fg  are  understood  to  be  the  position  vectors  of  the  aperture  and  source  at  time 
t  —  T, 

An  expression  for  the  temporal  cross-correlation  Raiaj  ('^)  is  obtained  using  the  spatial 
cross-correlation  expression  in  Eq.  (152)  by  substituting  Si  '(r)  for  S;,  substituting  A;(r) 
for  Ay2h  and  substituting  Ai{0)  for  Auu* 


Raiaj{r;fg,f,;vg,v,;vi) 


x*^(ni+i){[l  -  ^j(0)]aj}*^(n,+i){[l  -  ^/(T)]a;} 

,  X  J(„.+1){[1  -  A,(0)]x}J(„.+i){[l  -  x4,(r)]x}. 


(189) 
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Figure  36.  ABL  target  engagement  scenario.  This  figure  illustrates  the  particular  aper¬ 
ture  and  target  positions  and  velocities  used  to  quantify  temporal  cross¬ 
correlation  of  the  Zernike  coefficients  used  in  the  MMSE  tilt  estimator.  The 
position  of  the  target  object  is  specified  relative  to  the  position  of  the  aper¬ 
ture  by  the  parameter  d,  which  represents  the  distance  between  the  target 
and  aperture. 


The  functions  F,  G,  and  H  are  computed  as  in  Eq.  (148)-(150),  with  replaced  by 
which  is  the  angle  between  5/  '(r)  and  x. 


5.4  Beam-steering  performance  enhancement 

The  expression  for  the  temporal  cross-correlation  given  in  Eq.  (189)  may  be  used 
to  quantify  the  elements  of  the  MMSE  tilt  estimator  in  Eq.  (179).  In  order  to  perform 
these  calculations,  a  certain  geometric  specification  must  be  made.  The  geometry  used  to 
obtain  the  results  presented  in  this  section  is  shown  in  figure  36.  This  geometry  models 
the  expected  target  engagement  for  an  ABL  system  [10]. 

In  this  geometry,  the  aperture  is  located  at  a  height  of  13,000  m  above  a  fixed  origin 
specified  by  the  unit  vectors  x,  y,  and  z.  The  x  vector  is  in  the  direction  of  the  aperture’s 
velocity  vector,  Va,  whose  magnitude  is  220  m/s.  The  target  object  is  located  at  a  height 
of  15,300  m.  The  component  of  the  target  position  vector  in  the  y  direction  is  0.  The 
X  component  of  the  target’s  position  is  specified  by  the  parameter  d,  representing  the 


113 


distance  between  the  transmitting  aperture  and  the  target  object.  Since  the  difference  in 
the  z  components  of  Tq  and  is  fixed  to  be  2,300  m,  and  the  y  component  of  is  always 
0,  then  the  x  component  of  Tq  is  given  by: 

fo-x=  -^^277(2300"^.  (190) 

For  all  calculations,  the  aperture  diameter  is  1.75  m,  and  the  wind  velocity,  vi  was  set  to 
0,  and  d  =  100  km.  The  velocity  of  the  target  object,  is  1,000  m/s  in  the  z  direction. 

In  ABL  target  engagements,  a  beacon  signal  (measured  at  time  t  —  r)  is  generated 
by  laser  reflection  from  the  leading  surface  of  the  target  object.  The  system  attempts  to 
steer  the  transmitted  beam  at  time  t  to  the  same  physical  location  from  which  the  beacon 
was  received,  with  the  forward  motion  of  the  target  determining  the  impact  point  on  the 
object’s  body  (the  target’s  displacement  is  v^r).  Thus,  it  may  be  said  that  fo{t)  =  ^^(t  —  r). 
This  scenario  corresponds  to  Vg  =  0,  which  implies  that  rg{t)  =  rg{t  —  r)  in  Figure  35. 
Thus,  the  temporal  dependence  of  the  cross-correlation  is  determined  primarily  by  the 
aperture  velocity  Va^ 

5,4^ i  Turbulence  model  The  expression  for  Ra^aji'^)  given  in  Eq.  (189)  requires 
the  optical  turbulence  to  be  modeled  in  a  finite  number  of  discrete  layers  with  the  pa¬ 
rameters  Wi^  and  Lq.  For  the  ABL  target  engagement  considered  here,  we  have  chosen 
to  use  Hufnagel- Valley  (HV-21)  Cl{z)  profile.  The  weights  shown  in  Table  1  are  not  ade¬ 
quate  for  modeling  the  turbulence  since  is  distributed  differently  between  the  aperture 
and  target  than  was  the  case  for  the  HV-21  weights.  Over  the  altitudes  considered  here, 
decreases  in  a  nearly  linear  fashion  with  z.  Hence,  the  atmosphere  was  divided  into 
10  evenly-spaced  layers  between  the  aperture  altitude  (13,000  m)  and  the  target  altitude 
(15,300  m).  The  altitude  and  weight  for  each  layer  is  shown  in  Table  5.  For  the  results 
shown  in  this  chapter,  a  constant  value  of  Lq  =  100  m  was  chosen  for  all  layers. 

5,4^^  Higher-order  modal  cross-correlation.  Since  the  projected  aperture  sepa¬ 
ration  Si  '(r)  is  non-zero  for  r  >  0,  then  the  cross-correlation  effect  discussed  in  Chap¬ 
ter  III  will  also  be  manifested  in  the  cross-correlation  of  tilt  with  higher-order  modes  in 
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layer  1 

altitude  Zi  (m) 

Wi 

layer  1 

altitude  Zi  (m) 

Wi 

1 

13,209 

0.1281 

6 

14,255 

2 

13,418 

0.1216 

7 

14,464 

0.0904 

3 

8 

14,673 

0.0846 

4 

13,836 

9 

14,882 

0.0791 

5 

14,045 

0.1024 

10 

15,091 

0.0738 

Table  5.  Turbulence  modeling  for  ABL  target  engagement.  Altitudes  and  weighting  fac¬ 
tors  associated  with  10  discrete  layers  used  to  model  the  atmosphere  between 
the  aperture  and  the  target. 

an  ABL  target  engagement.  Using  Eq.  (189)  and  the  target  engagement  geometry  shown 
in  Figure  36,  it  can  be  demonstrated  that  for  many  higher-order  modes,  the  magnitude 
of  cross-correlation  with  tilt  increases  from  0  at  zero  time  delay  through  some  maximum 
value  before  decreasing  to  0  at  long  time  delays.  Thus,  including  these  higher-order  modes 
in  the  MMSE  tilt  estimator  will  enhance  its  performance  at  moderate  time  delays. 

The  cross-correlation  between  the  tilt  acquired  by  the  transmitted  beam  upon  prop¬ 
agation  to  the  target  object  and  time-delayed  modal  measurements  from  a  beacon  source 
are  shown  in  Figure  37.  In  these  figures,  (z  =  2,3)  represent  the  tilt  components  for 
propagation  to  the  target  object  at  time  t,  whereas  (j  =  2,3,  ...,15)  represent  the 
modal  measurements  from  the  beacon  for  tilt  through  Zernike  15  at  time  (t  —  r).  The 
cross-correlation  has  been  normalized  as  indicated  on  the  figures  to  express  a  correlation 
coefficient  in  the  range  -1  to  1.  The  autocorrelation  of  the  tilt  components  [i  =  j)  shown 
in  Figure  37  illustrate  the  expected  behavior  of  maximum  correlation  at  r  =  0  and  de¬ 
creasing  correlation  with  increasing  time  delay.  As  indicated  in  these  figures,  at  r  =  0  all 
higher-order  modes  considered  here  are  uncorrelated  with  tilt  except  the  orthogonal  com¬ 
ponents  of  coma  {j  =  7, 8).  The  magnitude  of  the  correlation  coefficient  for  all  other  modes 
increases  as  r  increases,  passing  through  some  maximum  value  before  going  to  0  at  large 
r.  The  enhanced  cross-correlation  between  tilt  and  higher-order  modes  for  r  >  0  indicates 
that  time-delayed  higher-order  measurements  may  be  used  in  the  MMSE  tilt  estimator  to 
improve  its  performance. 

5,4*3  Tilt  compensation  performance  results.  To  study  the  effect  of  higher- 
order  modes  in  the  MMSE  tilt  estimator,  let  us  consider  the  value  of  as  the  modal 
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Figure  37.  Cross-correlation  of  tilt  with  higher-order  modes  for  ABL  target  engagement. 

Correlation  results  are  shown  for  the  tilt  acquired  by  the  transmitted  beam 
upon  propagation  to  the  target  object  at  time  t  with  beacon  measurements 
through  Zernike  15  at  time  t  —  r  (uncorrelated  modes  omitted). 
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measurements  are  included.  For  the  target  engagement  scenario  shown  in  Figure  36,  the 
aperture-averaged  piston-removed  phase  variance  is  given  by: 

I  dpWi^  =  0.3828  .  (191) 

Since  all  elements  of  the  covariance  matrices  C^y  are  proportional  to 

[see  Eq.  (189)],  e\  may  be  normalized  for  the  relative  aperture  diameter  and  calculated  as 

follows: 

/n\-5/3 

4(^-j  =  0.3828  -(^-j  TR[(7,,C-iC,,].  (192) 

The  normalized  aperture-averaged  residual  phase  variance  was  quantified  accord¬ 
ing  to  Eq.  (192)  for  a  range  of  time  delay  r  using  various  mode  combinations  in  the 
measurement  vector  x.  The  result  of  these  calculations  is  shown  in  Figure  38.  For  all 
correction  schemes,  as  r  ^  0,  approaches  its  limiting  value  of  0.1342,  indi¬ 

cating  ideal  tilt  correction.  For  r  >  0,  the  normalized  error  increases  with  increasing  r. 
When  higher-order  modal  measurements  are  incorporated  into  the  MMSE  tilt  estimator, 
however,  increases  more  slowly  with  r  than  when  tilt  measurements  are  used 

alone. 

To  compare  performance  for  the  correction  schemes  studied,  consider  the  time  delay 
at  which  each  scheme  yields  e\{P =  0.1356;  that  is,  the  time  delay  at  which  the 
normalized  error  is  1%  larger  than  the  limiting  value  of  0.1342.  Let  us  call  this  time 
delay  the  tilt  correction  interval^  and  designate  it  r^.  If  tilt  correction  corresponding  to 
E‘^{D/ro)~^/^  <  0.1356  is  required  for  ABL  beam  steering,  then  represents  the  maximum 
interval  between  AO  correction  updates.  Thus,  we  may  think  of  fc  =  I/tc  as  the  operating 
frequency  or  bandwidth  of  the  AO  system.  The  values  of  and  fc  corresponding  to 
the  results  in  Figure  38  are  presented  in  Table  6.  As  higher  order  modes  are  included 
in  the  MMSE  tilt  estimator,  Tc  is  increased,  corresponding  to  a  reduced  value  for  /c. 
The  bandwidth  reduction,  as  a  percentage  of  the  bandwidth  using  tilt  measurements  only 
(modes  2-3),  is  shown  in  the  last  column  of  Table  6.  When  tilt  plus  higher  order  modal 
measurements  through  Zernike  15  are  included  in  the  tilt  estimator,  is  over  3  times  longer 
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Figure  38.  Normalized  aperture- averaged  residual  phase  variance  after  tilt  correction  for 
the  ABL  target  engagement.  The  value  of  is  normalized  to  (D/vo)^^^  as 
a  function  of  beacon  measurement  time  delay  r  for  the  MMSE  tilt  estimator 
using  the  modal  measurements  indicated.  The  dashed  line  shows  a  normalized 
error  that  is  1%  larger  than  the  limiting  value  of  0.1342. 

than  when  tilt  measurements  alone  are  included.  Correspondingly,  the  AO  beam-steering 
bandwidth  is  reduced  by  over  70%. 

For  a  given  value  of  the  relative  aperture  diameter  D/ro,  the  Strehl  ratio  performance 
of  an  ABL  beam-steering  AO  system  may  be  quantified  according  to  Eq.  (172).  The  results 
shown  in  Figure  38  for  were  applied  to  a  propagation  scenario  for  which 

D/tq  =  5,  The  corresponding  Strehl  ratio  performance  is  shown  in  Figure  39  as  a  function 
of  the  beacon  measurement  time  delay  r. 

As  is  illustrated  in  Figure  39,  greater  beam-projection  performance  is  maintained  at 
longer  beacon  measurement  time  delays  when  higher-order  modes  are  incorporated  into 
the  MMSE  tilt  estimator.  At  a  fixed  value  of  r,  the  Strehl  ratio  is  increased  by  using 
the  correlated  higher-order  modes.  In  order  to  maintain  a  particular  value  of  the  Strehl 
ratio,  longer  measurement  time  delays  may  be  used  in  the  system  if  higher-order  modes 
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MMSE  modes 

Tc  (msec) 

fa  (kHz) 

bandwidth  reduction 

2-3 

0.3839 

2.6048 

... 

2-4 

0.4742 

2.1088 

19.04% 

2-6 

1.0161 

0.9842 

62.22% 

2-10 

1.1065 

0.9038 

65.30% 

2-15 

1.3323 

0.7506 

71.18% 

Table  6.  ABL  AO  bandwidth  reduction  using  the  MMSE  estimator.  The  tilt  correction 
interval  Tc  corresponds  to  value  of  e\{D that  is  1%  larger  than  the 
limiting  value  of  0.1342  using  modal  measurements  in  the  MMSE  tilt  estimator 
as  indicated. 


are  utilized  in  estimating  the  tilt  components.  The  limiting  value  of  the  Strehl  ratio 
(approximately  0.14)  after  tilt  compensation  for  the  relative  aperture  diameter  D/ro  —  5 
indicates  that  in  many  practical  target  engagements,  higher-order  compensation  may  be 
required  in  the  ABL  AO  system  to  maintain  beam-projection  performance.  If  higher-order 
compensation  is  required,  the  MMSE  estimator  will  provide  for  enhanced  performance 
when  the  effect  of  measurement  time  delay  is  significant. 


5.5  Summary 

Beam  steering  for  ABL  target  engagements  should  include  adaptive-optical  tilt  cor¬ 
rections  to  compensate  for  the  degrading  effect  of  phase  aberrations  induced  by  turbu¬ 
lence  upon  propagation  from  the  system  aperture  to  the  target  object.  Due  to  the  finite 
bandwidth  of  any  such  AO  system,  the  tilt  correction  will  be  based  upon  time-delayed 
measurements  of  the  beacon  phase.  The  analysis  presented  in  this  chapter  indicates  that 
due  to  anisoplanatism,  tilt  correction  is  suboptimal  when  time-delayed  tilt  measurements 
are  used  directly  in  the  beam-steering  AO  system.  Optimal  correction  is  realized  only 
when  tilt  components  are  obtained  using  an  MMSE  estimator  given  time-delayed  beacon 
measurements.  The  MMSE  estimator  permits  the  use  of  tilt  plus  higher-order  modal  mea¬ 
surements  in  forming  the  tilt  estimate.  The  elements  of  the  MMSE  estimator  may  be 
quantified  by  calculating  the  cross-correlation  of  target  object  and  beacon  Zernike  modes 
as  a  function  of  the  time  delay  r. 
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Figure  39.  Strehl  ratio  performance  of  ABL  beam-steering  AO  system  for  D/r^  =  5. 

Performance  is  quantified  as  a  function  of  the  measurement  time  delay  r  for 
MMSE  tilt  estimation  using  the  modal  measurements  indicated. 

The  analysis  presented  in  this  chapter  for  an  ABL  target  engagement  demonstrates 
that  while  tilt  is  uncorrelated  with  most  higher-order  modes  at  r  =  0,  the  correlation  passes 
through  some  extremum  at  r  >  0  before  going  to  0  at  large  values  of  r.  Thus,  incorporating 
higher-order  modes  in  the  tilt  estimator  will  enhance  tilt  correction  performance.  For  the 
ABL  target  engagement  scenario,  tilt  correction  performance  was  quantified  when  various 
mode  combinations  were  included  in  the  MMSE  tilt  estimator.  These  results  illustrate 
that  enhanced  performance  is  obtained  by  including  delayed  modal  measurements  through 
Zernike  15.  Hence,  longer  measurement  time  delays  may  be  used  while  maintaining  a 
normalized  residual  phase  variance  that  is  1%  larger  than  its  limiting  value.  Correspond¬ 
ingly,  the  ABL  beam-steering  AO  bandwidth  may  be  reduced  when  higher-order  modes  are 
included  the  MMSE  tilt  estimator.  By  including  beacon  measurements  through  Zernike 
15,  the  beam-steering  AO  bandwidth  is  reduced  by  over  70%  from  the  bandwidth  when 
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using  tilt  measurements  alone.  At  a  fixed  bandwidth,  the  Strehl  ratio  is  increased  by 
incorporating  higher-order  modes  into  the  MMSE  estimator  for  wave-front  tilt. 
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VL  Adaptive- Optical  Interferometry 


6. 1  Introduction 


Michelson  stellar  interferometry,  also  referred  to  as  amplitude  interferometry,  pro¬ 
vides  for  imaging  of  astronomical  objects  with  higher-resolution  than  can  be  achieved  using 
single  aperture  systems  [32,60,77].  This  imaging  technique  is  based  upon  the  Van  Cittert- 
Zernike  theorem  [9],  which  relates  the  visibility  and  phase  of  a  fringe  in  the  interferometer 
(referred  to  as  the  complex  visibility)  to  the  normalized  Fourier  transform  of  the  irradiance 
distribution  associated  with  the  object  of  interest.  Mathematically,  the  Van  Cittert-Zernike 
theorem  implies  that  [32] : 


IAo(^| 


Jd7fIo{rf)exp{jj^s-Tf) 

Jdffloirf) 


(193) 


where  Po  is  the  normalized  mutual  intensity  of  the  object  of  interest,  s  is  the  baseline  of 
the  interferometer,  Iq  is  the  object  irradiance,  A  is  the  observation  wavelength,  and  L  is 
the  distance  from  the  object  to  the  interferometer  plane.  Recall  from  Eq.  (71)  that  when 
an  interferometer  is  unaffected  by  turbulence,  then  the  fringe  visibility  V  =  Prom 

Eq.  (193),  it  can  be  concluded  that  V  for  a  given  baseline  s  is  equal  to  the  modulus  of  the 
normalized  object  spectrum  evaluated  at  the  spatial  frequency  s/XL,  Thus,  measurement 
of  the  fringe  visibility  at  many  baseline  values  provides  for  reconstruction  of  the  object 
spectrum. 

Interferometers  have  been  constructed  and  operated  at  many  locations  throughout 
the  world  for  the  purpose  of  interferometric  imaging  [13,52,73].  The  use  of  large  apertures 
in  interferometric  applications  at  optical  wavelengths  is  limited  by  the  degrading  effects 
of  atmospheric  turbulence.  When  atmospheric  turbulence  is  present,  the  fringe  pattern 
in  the  interferometer  moves  randomly,  corrupting  the  relationship  between  the  phase  of 
the  fringe  and  the  phase  of  the  Fourier  transform  of  the  object  irradiance,  and  reducing 
fringe  visibility  measurements  for  long-exposure  fringe  measurement  [66].  To  enhance  in¬ 
terferometer  performance  in  the  presence  of  atmospheric  turbulence,  adaptive-optics  phase 
compensation  has  been  proposed  [48]  and  implemented.  Many  of  these  applications  involve 
fringe-tracking  methods  for  control  of  the  AO  systems  [12,46]. 
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In  this  chapter,  practical  limitations  to  the  use  of  AO  compensation  in  interferometry 
are  explored.  The  effects  of  angular  anisoplanatism  are  investigated  in  Section  6.2.  The 
performance  degradation  related  to  angular  anisoplanatism  is  quantified  for  various  levels 
of  modal  compensation.  In  Section  6.3,  drawing  upon  the  optimal  compensation  theory 
developed  in  Chapter  III,  estimation  theory  is  applied  to  AO  compensation  for  interfer¬ 
ometry  to  partially  recover  the  performance  degradation  due  to  anisoplanatism.  Results 
are  presented  in  Section  6.4  to  quantify  the  visibility  enhancement  obtained  when  optimal 
AO  compensation  is  employed.  The  chapter  is  summarized  in  Section  6.5 

6,2  Angular  anisoplanatism  in  adaptive  optics  for  interferometry 

Interferometric  techniques  are  often  employed  for  high-spatial-resolution  imaging  of 
astronomical  objects.  These  objects  are  typically  dim,  and  therefore  may  not  be  bright 
enough  to  provide  for  wave-front  sensing.  In  this  case,  a  nearby  NGS  may  be  used  as 
the  AO  beacon.  Since  there  is  an  angular  separation  between  the  beacon  and  object  of 
interest,  turbulence-induced  anisoplanatism  will  affect  AO  compensation  performance  in 
each  aperture  of  the  interferometer.  To  quantify  the  effect  of  angular  anisoplanatism  in 
an  interferometric  system,  consider  the  system  diagram  in  Figure  40,  which  shows  an 
adaptive-optical  interferometer  affected  by  angular  anisoplanatism.  For  simplicity,  the 
object  is  located  on  the  zenith,  with  the  beacon  located  at  an  angle  0^  with  respect  to  the 
object  in  the  direction  of  the  aperture  separation  vector  s. 

The  performance  of  the  interferometer  is  quantified  by  the  fringe  Strehl  ratio,  which 
according  to  Eq.  (99)  is  determined  by  the  aperture-averaged  residual  differential  phase 
variance  6%  for  compensation  for  piston  through  the  iVth-order  Zernike  mode.  Recall  that 
the  FSR  decreases  with  increasing  values  of  Suppose  that  the  beacon  modal  phase 
measurements  in  each  aperture  are  used  as  the  conjugate  phase  for  each  DM.  In  this  case, 
the  difference  of  the  conjugate  phase  coefficients  Ad^i  is  given  by  Add  =  ^bu  —  Thus, 
from  Eq.  (96),  is  expanded  as: 
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Figure  40.  Angular  anisoplanatism  in  adaptive  optics  for  interferometry.  When  an  off- 
axis  source  is  used  for  wave  front  sensing,  angular  anisoplanatism  is  introduced 
into  each  adaptive-optics  system,  resulting  in  reduced  fringe  visibility. 
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r  ^ 

=  /  dpW{^  ([<f>oim  -  <f>02{R^?)  -  E  ((^oli  -  ao2if) 

•'  i=2 

N 

”1“  ^  ^  ([(^oli  ^o2i)  ^62i)]  )  j 

i=2 

If  the  terms  inside  the  summation  of  Eq.  (194)  are  expanded,  the  expression  for  is 
simplified  as  follows: 

si  =  jdpW{^l^[\<j>ol{Rp)-<j>02{Rp)]^) 

N 

~~  ^  ^  {2  “  (®oli^62i)  ““  {®o2«^61i)  "I”  {^o2i^b2i)] 

i=l 

“  i^bli)  +  2  (a6Haft2i)  -  («i62i)}  •  (19^) 


Using  the  expression  in  Eq.  (152),  for  the  spatial  cross-correlation  of  Zernike  coefficients 
Rauua^ij^  each  of  the  terms  of  Eq.  (195)  can  be  quantified.  If  the  position  of  the  object 
of  interest  is  designated  f^,  and  the  position  of  the  beacon  is  designated  then  the 
correlation  calculations  required  in  Eq.  (195)  are  as  follows: 


{cioli^bli)  —  -Sauitav2i  (^05  ^ol  ?  ^6?  ^al)j  ifloH^b2^  —  ■®o„iiav2t  (^OJ  )  ^^>5  ^^2)5 


(Uo2i®61i)  —  ■®Oui»av2t  (^05  ^a2)  ^6?  ^al)?  {^o2i^b2i)  -®aui*au2t  (^o?  *^02)  ^02)? 


(^61i)  =  Bauua.2i{n,raun,rai),  {atiiai,2i)  =  Ba^ua.^i{rb,raun,fa2), 


(196) 


and  {al^i)  =  5a,,,a.2i(n,7^a2,n,fa2). 

Eq.  (195)  also  requires  the  aperture-averaged  differential  phase  variance  associated  with 
the  object  of  interest  to  be  quantified.  Including  the  effects  of  finite  outer  scale  Lq^  this 
term  is  calculated  as  follows: 
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(197) 


I  dpW{^  {[<f>oi{R^  -  MR^f) 

-  iff 

X  (§x) 

where  s/D  is  the  relative  aperture  separation  and  Xq  =  tvD/Lq,  as  in  Eq.  (152). 

As  indicated  in  Eq.  (197),  the  aperture-averaged  differential  phase  variance  is  pro¬ 
portional  to  Prom  Eq.  (152),  it  can  be  seen  that  each  term  of  Eq.  (195)  is 

proportional  to  Thus,  6%  is  also  proportional  to  {D/voY^^,  and  compensa¬ 

tion  performance  can  be  quantified  in  terms  of  the  normalized  aperture-averaged  residual 
differential  phase  variance,  6^(i)/ro)“^/^. 

For  a  range  of  beacon  angular  offset  0^,  6%{D/roY^^  was  calculated  according  to 
Eq.  (195)  for  modal  compensation  through  order  AT  =  15,  =  36,  iV  =  66,  N  =  105,  and 

AT  =  153,  corresponding  to  Zernike  polynomial  radial  order  n  =  4,  n  =  7,  n  =  10,  n  =  13, 
and  n  =  16,  respectively.  The  results  of  these  calculations  are  shown  in  Figure  41(a)-(b) 
for  interferometer  baselines  5  =  |^  corresponding  to  s/D  =  4:  and  s/D  —  6.  The  results  for 
interferometer  baselines  s/D  =  8  and  s/D  —  10  are  shown  in  Figure  42(a)-(b).  All  results 
have  been  calculated  with  an  aperture  diameter  D  =  1  m,  and  Lq/D  =  100.  The  beacon 
angular  offset  6f,  has  been  normalized  to  the  isoplanatic  angle  ^05  which  for  the  HV-21 
turbulence  profile  employed  in  these  calculations  is  approximately  7  /zrad,  as  indicated  in 
Table  1.  The  maximum  angle  considered  in  these  results  is  9^,  =  AOq, 

As  can  be  seen  in  each  graph  of  Figure  41  and  Figure  42,  at  each  baseline  and  for 
all  orders  of  compensation,  6%{D/ro)^^^^  increases  with  increasing  beacon  offset.  This  in¬ 
crease  represents  the  magnitude  of  the  angular  anisoplanatic  effect  for  the  adaptive  optical 
interferometer.  For  a  given  beacon  offset,  the  angular  anisoplanatic  effect  may  be  reduced 
by  increasing  the  order  N  of  the  atmospheric  compensation.  For  larger  values  of  9^/00^ 
however,  the  marginal  decrease  in  angular  anisoplanatism  with  increasing  N  is  small. 
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Figure  41.  Effect  of  modal  compensation  on  angular  anisoplanatism  in  adaptive-optical 
interferometry  for  baselines  s/D  =  4  and  s/D  =  6.  Normalized  residual 
differential  phase  variance  6jf{D/roy^^  is  shown  for  relative  beacon  offset  6b/9o 
with  correction  order  N  =  15, 36, 66, 105, 153.  Dashed  lines  show  =  1  rad^ 
for  D/ro  =  5, 10, 15, 20, 25.  Results  are  for  Lq/D  =  100. 
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(b) 

Figure  42.  Effect  of  modal  compensation  on  angular  anisoplanatism  in  adaptive-optical 
interferometry  for  baselines  (a)  s/D  =  8  and  (b)  s/D  =  10.  Normalized 
residual  differential  phase  variance  6lf{D/ro)^^^  is  shown  for  relative  beacon 
offset  Ob/Oo  with  correction  order  N  =  15,36,66,105,153.  Dashed  lines  show 
=  1  rad^  for  D/tq  =  5, 10, 15, 20, 25.  Results  are  for  L^jD  =  100. 
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Consider  the  maximum  angle  for  which  <5^  <  1  rad^,  which  implies  that  FSR  > 
0.6.  This  maximum  angle  specifies  the  range  of  angular  offsets  at  which  the  FSR  is  not 
significantly  reduced  by  turbulence-induced  anisoplanatism  in  the  AO  system.  The  value 
of  corresponding  to  =  1  rad^  is  shown  in  each  graph  of  Figure  41  and 

Figure  42  by  dashed  lines  for  the  relative  aperture  diameters  D/tq  =  5,10,15,20,25,30. 
For  the  20  to  35  cm  Tq  values  experienced  at  the  best  observatory  locations  (for  A  =  500 
nm),  these  relative  aperture  diameters  are  applicable  to  1  to  10  m  class  telescopes  used  as 
baseline  apertures. 

For  the  least  stressing  turbulence  condition  indicated  in  each  graph  of  Figure  41  and 
Figure  42  {D/tq  =  5),  the  lowest-order  correction  considered  {N  =  15)  is  an  adequate 
level  of  compensation  for  Of,  <  1.5  Oq  at  all  baseline  values  studied.  When  the  highest-order 
correction  considered  {N  =  153)  is  employed  for  D/tq  =  5,  a  beacon  offset  of  up  to  3 
results  in  a  FSR  value  no  less  than  0.6. 

For  the  most  stressing  turbulence  condition  indicated  {D/tq  —  30),  even  the  highest 
order  correction  considered  {N  =  153)  is  insufficient  to  provide  <  1  rad^  at  any  beacon 
offset  for  the  baseline  values  s/D  =  4  and  s/D  =  10.  At  s/D  =  6,  compensation  through 
N  =  153  reduces  the  residual  variance  to  a  sufficient  level  only  when  9^  0.  At  s/D  =  8^ 

the  maximum  allowable  beacon  offset  for  N  =  153  is  approximately  0.2  9o. 

It  is  curious  and  perhaps  counter-intuitive  that  6%{D/rQ)~^^^  should  be  smaller  at  the 
intermediate  baseline  values  studied  {s/D  =  6,  s/D  =  8)  than  at  s/D  =  4  for  0b/Oo  <  1,  as 
can  be  seen  by  comparing  the  results  in  Figure  41(b)  and  Figure  42(a)  with  Figure  41(a). 
This  effect  can  be  understood  by  considering  the  limiting  form  of  <5^  as  ^  0?  which 
is  given  in  Eq.  (100),  and  rewritten  here  in  an  expanded  form: 

Aaii-^Aa,.  "  E 

N  N 

~  E  ((®o2i)^)  +  2  ^  {O'oliO-oii)  ■  (198) 

i=2  i=2 

As  can  be  seen  from  Eq.  (198),  the  interaperture  correlation  of  the  Zernike  coefficients  (the 
final  summation)  adds  to  the  limiting  value  of  6^.  For  relative  outer  scale  L^/D  =  100,  the 
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interaperture  correlation  of  the  x-tilt  (Zernike  2)  is  negative  and  increasing  in  magnitude 
through  some  extremum  value  for  intermediate  values  of  s/D  [76]  (see  also  Figure  13). 
Thus,  the  interaperture  term  for  z  =  2  in  Eq.  198  is  negative  and  large  compared  to  all 
other  terms  at  the  intermediate  values  of  sfO,  Therefore,  8]^  is  reduced  from  its  value  at 
smaller  relative  baselines. 

6,3  Application  of  theory  to  optical  interferometric  imaging 

The  results  of  Section  6.2  illustrate  that  angular  anisoplanatism  results  in  reduced 
FSR  values  by  increasing  the  aperture-averaged  residual  differential  phase  between  the 
baseline  apertures.  To  recover  some  of  the  performance  lost  due  to  anisoplanatism,  the 
optimal  compensation  theory  developed  in  Chapter  III  can  be  applied  directly.  Recall  that 
the  optimal  compensation  theory  was  developed  for  the  interferometric  system  shown  in 
Figure  11.  In  Chapter  IV  and  Chapter  V,  the  general  results  of  Chapter  III  were  specialized 
to  imaging  and  beam-projection  systems.  For  an  interferometric  system,  however,  the 
results  of  the  general  analysis  can  be  used  immediately. 

When  modal  beacon  measurements  ai,ii  and  ab2i  are  obtained  by  wave-front  sensing 
in  aperture  1  and  aperture  2,  respectively,  the  minimum  value  of  is  given  by  replacing 
Auoi  with  in  Eq.  (96)  as  follows: 

min(^^)  =  f  d^W{^  -  <f>o2{R^?)  -  E 

^  i=2 

N 

+  E  “  AOoi]^)  ,  (199) 

where  Ad^i  is  the  MMSE  estimate  of  Aa^i  =  {oou  —  tto2i))  the  difference  of  object  phase 
coefficients  between  aperture  1  and  aperture  2.  The  value  of  Ad^i  is  estimated  from  the 
measured  value  of  Aa^i  ~  (a^H  “  <^b2i)  using  the  MMSE  estimator  of  Eq.  (106). 

Since  6%  depends  only  upon  the  difference  of  the  conjugate  phase  applied  to  the  DM 
in  each  aperture,  then  without  loss  of  generality,  the  deformable  mirror  2  may  be  removed 
and  replaced  with  a  plane  mirror,  as  shown  in  Figure  43.  Replacing  the  deformable  mirror 
in  aperture  2  with  a  plane  mirror  corresponds  to  setting  do2i  =  0,  and  hence,  Ad^i  =  cLou, 
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object  beacon  object  beacon 


Figure  43.  Optimal  compensation  for  anisoplanatism  in  adaptive-optical  interferometry. 

For  optimal  compensation,  only  one  DM  is  employed.  The  conjugate  phase 
placed  on  the  DM  is  estimated  from  the  difference  in  phase  coefficients  between 
wave-front  sensor  1  and  wave-front  sensor  2. 
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Therefore,  to  obtain  the  minimal  value  of  6%  given  in  Eq.  (199),  the  MMSE  estimates  of 
Adoi  represent  the  conjugate  phase  that  should  be  placed  on  the  single  deformable  mirror 
for  optimal  compensation  in  an  adaptive-optical  interferometer. 

6.4  Fringe  visibility  enhancement  through  optimal  atmospheric  compensation 

Figure  41  and  Figure  42  illustrate  the  effect  of  angular  anisoplanatism  on  the  aperture- 
averaged  residual  differential  phase  variance  6%,  The  results  were  normalized  for  the  rela¬ 
tive  aperture  D/r^,  For  a  given  value  of  J3/ro,  the  fringe  Strehl  ratio  is  calculated  according 
to  Eq.  (99),  repeated  here  for  clarity: 

FSR  =  exp(-^^^)  .  (200) 

Prom  the  results  in  Figure  41  and  Figure  42,  it  can  be  determined  that  the  relative  aperture 
diameter  D/r^  =  5  requires  modal  compensation  through  approximately  N  =  15.  Consider 
the  FSR  results  shown  in  Figure  44.  This  plot  gives  the  FSR  value  at  s/D  =  8  for  D/ro  =  S 
and  Lq/D  =  100  with  N  —  6^  10, 15,21  when  beacon  measurements  are  employed  directly 
to  form  the  conjugate  phase.  The  beacon  offset  Of,  is  normalized  to  6q,  and  the  results 
were  calculated  using  D  =  1  m.  As  can  be  seen  in  Figure  44,  the  compensation  orders 
indicated  yield  FSR  values  in  the  range  from  0.4  to  approximately  0.8  at  Oi,/Oo  =  0.  The 
FSR  for  all  compensation  orders  falls  off  as  0f,/0Q  increases,  resulting  in  a  fringe  Strehl 
ratio  no  greater  than  0.5  for  beacon  offsets  greater  than  AOq,  At  the  largest  angular  offset 
considered  {0f,/6Q  =  10),  the  FSR  is  less  than  0.05  for  all  compensation  orders  considered. 

It  is  instructive  to  compare  the  effect  of  angular  anisoplanatism  on  the  fringe  Strehl 
ratio  in  an  interferometer  with  its  effect  on  the  Strehl  ratio  in  an  imaging  system.  This 
can  be  done  by  comparing  the  results  shown  in  Figure  44  with  the  angular  anisoplanatism 
degradation  results  of  Figure  10(a)  in  Chapter  II  for  each  order  of  compensation  con¬ 
sidered.  For  convenience.  Figure  45  shows  the  FSR  performance  curve  for  A?”  =  21  in 
Figure  44  plotted  with  the  SR  performance  curve  for  AT  =  21  in  Figure  10(a).  This  graph 
illustrates  that  while  the  FSR  is  greater  than  the  SR  for  Ot/Oo  =  0  (due  to  the  interaper¬ 
ture  tilt  correlation  effect  discussed  previously),  the  FSR  degradation  is  more  severe  than 
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Figure  44.  Effect  of  angular  anisoplanatism  on  fringe  Strehl  ratio  in  adaptive  optics  for 
interferometry.  Value  of  FSR  is  shown  for  relative  beacon  offset  ^it  inter¬ 
ferometer  baseline  s/D  =  8  with  the  compensation  order  indicated.  Results 
are  for  D/r^  =  5  and  L^jD  =  100. 

the  SR  degradation  with  increasing  beacon  offset  angle.  Thus,  angular  anisoplanatism 
imposes  more  severe  beacon  offset  limits  in  interferometric  applications  than  for  imaging 
applications  with  adaptive  optics. 

As  was  illustrated  by  the  results  presented  in  Figure  44,  angular  anisoplanatism  can 
significantly  reduce  the  fringe  visibility  for  beacon  offsets  that  ar^  several  times  the  isopla- 
natic  angle.  The  decreased  fringe  visibility  limits  the  ability  to  recover  object  spectrum 
information  for  interferometric  imaging  applications.  To  mitigate  the  degrading  effects  of 
angular  anisoplanatism,  the  optimal  compensation  theory  may  be  applied,  as  was  discussed 
in  Section  6.3. 

When  the  conjugate  phase  coefficients  are  estimated  from  beacon  modal  measure¬ 
ments,  is  given  by  Eq.  (111).  If  the  measured  beacon  modal  coefficients  are  used  directly 
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Figure  45.  Comparison  of  angular  anisoplanatism  effect  on  performance  of  interferometric 
and  imaging  systems.  The  solid  line  shows  the  fringe  Strehl  ratio  for  an 
interferometer  with  s/D  —  8  and  compensation  order  =  21.  The  dashed 
line  shows  the  Strehl  ratio  of  an  imaging  system  with  compensation  order 
AT  =  21.  Results  are  for  I?/ro  =  5  and  Lq/D  =  100. 


to  form  the  conjugate  phase,  then  6%  is  given  by  Eq.  (112).  Since  6%  is  minimized  by  using 
the  MMSE  estimator,  then  the  FSR  will  be  larger  when  the  estimator  is  employed  than 
when  the  beacon  measurements  are  employed  directly.  To  express  the  increase  in  the  fringe 
Strehl  ratio  that  results  from  optimal  compensation,  the  FSR  gain  GpgR  defined  as: 

exp  [— I  (measured:  o^)J 

From  the  expressions  given  in  Eq.  (Ill)  and  Eq.  (112),  GpgR  expressed  using  the 

vector  notation  defined  in  Chapter  III: 


GpsR  =  exp  j-^TR  [2Cy,  -  C,,  -  | . 


(202) 
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Since  each  element  of  Cy^  and  C^x  is  proportional  to  {D/roY^^^  then  Gpsj^  will  also 
depend  upon  the  relative  aperture  diameter  D/tq, 

The  value  of  Gpgj^  was  quantified  for  beacon  offset  angles  from  Ot/Oo  =  0  to 
^b/^o  =  10  for  D/tq  =  5  and  Lq/D  =  100.  The  results  of  these  calculations  are  shown 
in  Figure  46(a)  and  Figure  46(b)  for  iV  =  6  and  N  =  10,  respectively.  Figure  47(a)  and 
Figure  47(b)  show  the  FSR  gain  results  for  iV  =  15  and  N  =  21,  respectively.  These 
calculations  were  performed  using  D  =  1  m  with  the  4-layer  HV-21  turbulence  profile  in 
Table  1. 

For  modal  compensation  through  iV  =  6  and  N  =  10,  the  plots  in  Figure  46  show 
the  value  of  Gpgj^  for  the  relative  interferometer  baselines  s/D  =  4,6,8,10.  As  can  be 
seen  from  these  plots,  GpsR  is  appreciably  greater  than  1  for  beacon  offset  angles  greater 
than  6q,  The  FSR  gain  increases  as  the  beacon  angular  offset  increases  for  all  relative 
baseline  values.  The  gain  is  also  higher  for  the  lower  values  of  s/J9  at  all  values  of 
indicating  that  the  enhanced  cross-correlation  of  the  beacon  measurements  and  object 
phase  coefficients  at  lower  values  of  s/D  reduces  6]^  using  the  MMSE  estimator.  At  the 
largest  beacon  angular  offsets  considered  {Ob/Oo  =  10),  the  FSR  is  increased  by  15%  to 
approximately  25%,  depending  upon  the  relative  baseline  s/D, 

For  modal  compensation  through  iV  =  15  and  N  =  21,  the  plots  in  Figure  47  show 
the  value  of  Gpgj^  for  the  relative  interferometer  baselines  s/D  =  4,6,8,10.  As  was  the 
case  for  =  6  and  N  =  10,  Gpgj^  is  greater  than  1  for  beacon  offset  angles  greater  than  9q 
at  all  values  of  s/D,  The  FSR  gain  increases  as  the  beacon  angular  offset  increases  for  all 
relative  baseline  values.  Again,  the  gain  is  higher  for  the  lower  values  of  s/D  at  all  values 
of  Ob/9o-  At  9b /9o  =  10,  the  FSR  is  increased  by  over  25%,  for  s/D  =  4.  By  comparing  the 
plots  of  Figure  46  with  the  plots  of  Figure  47,  it  can  be  seen  that  the  FSR  gain  is  increased 
for  the  higher  orders  of  compensation.  This  indicates  that  the  MMSE  estimator  is  able  to 
partially  recover  some  of  the  performance  degradation  due  to  anisoplanatism  even  for  the 
highest  aberration  orders  considered. 

While  the  FSR  gain  is  greatest  for  larger  values  of  9b /9o,  the  actual  increase  in  the 
FSR  may  be  small,  since  the  FSR  is  small  for  these  beacon  angular  offset  values  (see 
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Figure  46.  Fringe  Strehl  ratio  gain  for  compensation  order  (a)  TV  =  6  and  (b)  N  =  10. 

FSR  gain  Gpgj^  is  shown  for  relative  beacon  offset  6b /Oq  at  interferometer 
baselines  s/D  =  4, 6, 8, 10  with  D/tq  =  5  and  Lq/D  =  100. 


Figure  44).  In  fact,  the  FSR  increase  may  be  larger  at  moderate  values  of  06 /^O)  since  the 
FSR  value  is  larger.  Regardless  of  the  FSR  increase,  >  1  indicates  an  increase  in 

the  fringe  visibility,  and  therefore  a  partial  recovery  of  the  performance  degradation  due 
to  angular  anisoplanatism.  It  should  be  pointed  out  here  that  since  GpgR  depends  upon 
D/t-q,  the  increase  in  FSR  is  expected  to  be  greater  for  the  higher  orders  of  atmospheric 
compensation  which  are  required  for  larger  values  of  relative  aperture  diameter 

6,5  Summary 

High-resolution  imagery  is  obtained  from  an  optical  interferometer  by  measuring 
the  visibility  of  the  fringe  obtained  when  light  from  two  apertures  is  combined,  and  the 
resulting  irradiance  pattern  is  measured  on  the  focal  plane  of  the  system.  According  to 
the  Van  Cittert-Zernike  theorem,  the  fringe  visibility  gives  the  normalized  spectrum  of 
the  object  irradiance  distribution  evaluated  at  a  spatial  frequency  which  is  proportional 
to  the  baseline  separation  of  the  apertures.  Thus,  long  baselines  provide  for  high  spatial- 
frequency  sampling  of  the  object  spectrum.  Large  apertures  may  be  employed  to  increase 
the  irradiance  on  the  focal  plane  of  the  system  in  order  to  enhance  visibility  measurements. 
Due  to  the  increased  magnitude  of  turbulence  effects  in  large  apertures  (increased  D/ro), 
however,  the  fringe  visibility  is  decreased,  and  object  spectrum  information  cannot  be 
accurately  measured.  In  this  case,  adaptive  optics  may  be  employed  to  compensate  for 
turbulence  effects,  restoring  the  benefits  of  using  large  apertures  in  the  interferometer. 

If  the  object  being  observed  by  an  interferometer  is  too  dim,  then  it  may  not  be 
bright  enough  to  provide  for  wave-front  sensing  in  the  AO  loop.  In  this  case,  a  nearby 
NGS  beacon  must  be  employed.  The  angular  separation  between  the  object  of  interest 
and  the  AO  beacon  results  in  a  degradation  of  compensation  performance  due  to  aniso¬ 
planatism.  In  this  chapter,  the  angular  anisoplanatic  effect  was  characterized  using  modal 
analysis  techniques,  where  the  performance  degradation  is  quantified  by  the  correlation  of 
the  appropriate  object  and  beacon  Zernike  coefficients.  For  several  values  of  the  relative 
baseline  s/D^  the  aperture-averaged  residual  differential  phase  variance  was  calculated 
over  a  range  of  beacon  angular  offsets  06-  The  maximum-permissible  beacon  angular  offset 
depends  upon  the  relative  aperture  diameter  D/tq  and  the  order  of  compensation  N,  For 
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D/ro  =  5,  the  beacon  angular  offset  limit  is  between  1.5  6q  and  3  Sqj  depending  upon  the 
compensation  order.  For  D/tq  =  30,  the  angular  limit  was  approximately  0.2  Oq  for  the 
highest-order  compensation  considered  {N  —  153).  In  general,  performance  degradation 
due  to  angular  anisoplanatism  is  more  severe  for  interferometric  applications  than  for  imag¬ 
ing  applications,  as  the  fringe  Strehl  ratio  drops  off  more  rapidly  with  increasing  beacon 
offset  angle  than  the  Strehl  ratio  of  an  imaging  AO  system. 

The  optimal  compensation  theory  developed  in  Chapter  III  can  be  applied  to  pro¬ 
vide  enhanced  fringe  visibility  in  an  adaptive-optical  interferometer  suffering  performance 
degradation  due  to  angular  anisoplanatism.  For  optimal  compensation,  the  beacon  mea¬ 
surements  are  not  used  directly  for  compensation  in  each  aperture,  but  instead,  an  optimal 
conjugate  phase  is  formed  in  one  aperture  based  upon  the  difference  of  beacon  measure¬ 
ments  in  each  aperture.  Application  of  the  optimal  conjugate  phase  results  in  an  increased 
FSR,  quantified  by  the  fringe  Strehl  ratio  gain  GpsR'  The  value  of  G^SK  depends  upon 
the  compensation  order  the  relative  baseline  s/D^  and  the  relative  aperture  diameter 
D/tq.  It  was  shown  that  for  A”  =  21  at  the  beacon  offset  Ot/^o  =  10,  optimal  compensation 
results  in  an  FSR  increase  of  15%  to  over  25%,  depending  upon  the  relative  baseline  s/D. 
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VIL  Conclusions  and  Recommendations 


7 A  Introduction 

The  principle  objective  of  this  research  was  to  establish  a  common  framework  for 
the  analysis  of  anisoplanatism  in  imaging,  beam-projection,  and  interferometric  adaptive- 
optical  systems  and  to  apply  optimal  AO  compensation  techniques  in  each  application. 
This  objective  was  motivated  by  the  inadequacy  of  current  analysis  methods  to  address 
the  degrading  effects  of  anisoplanatism  in  emerging  adaptive  optics  applications  such  as 
beam  steering  and  control  for  the  Airborne  Laser  and  atmospheric  compensation  for  high- 
resolution  interferometric  imaging.  Specifically,  this  research  was  intended  to  assess  the 
performance  degradation  caused  by  anisoplanatism,  and  to  quantify  the  performance  im¬ 
provement  that  could  be  obtained  using  optimal  atmospheric  compensation. 

In  Chapter  I  of  this  dissertation,  an  overview  of  anisoplanatism  in  adaptive  optics  was 
presented,  along  with  the  historical  background,  research  motivation,  problem  statement, 
and  key  research  results.  Chapter  II  provided  background  knowledge  in  atmospheric  and 
adaptive  optics,  highlighting  the  effects  of  turbulence  in  optical  systems.  Optimal  com¬ 
pensation  theory  was  laid  out  in  Chapter  III,  including  the  calculation  of  the  spatial  cross¬ 
correlation  of  Zernike  phase  coefficients,  which  characterizes  anisoplanatism  in  a  general 
sense.  The  application  of  the  optimal  compensation  theory  to  ground-based  adaptive  optics 
for  imaging  was  presented  in  Chapter  IV.  Chapter  V  was  devoted  to  exploring  the  optimiza¬ 
tion  of  Airborne  Laser  beam-steering  adaptive  optics.  Angular  anisoplanatic  effects  and 
optimal  compensation  in  adaptive-optical  interferometry  were  discussed  in  Chapter  VI.  In 
the  present  chapter,  the  significant  advances  and  results  of  this  research  are  summarized, 
conclusions  are  drawn,  and  recommendations  for  future  research  are  offered. 

7.2  Significant  advances 

This  research  has  resulted  in  the  following  significant  advances  regarding  the  under¬ 
standing  of  anisoplanatism  in  adaptive  optics: 

•  A  unified  approach  to  the  analysis  of  anisoplanatism  in  imaging,  dynamic  beam- 

projection,  and  interferometric  systems. 
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•  An  optimal  compensation  theory  for  adaptive-optical  systems  suffering  generalized 
anisoplanatic  degradation,  based  on  modal  decomposition  of  the  turbulence-induced 
phase  and  established  optimal  estimation  techniques. 

•  The  first  known  theoretical  and  experimental  validation  of  a  modal  cross-correlation 
enhancement  effect  that  occurs  in  systems  with  moderate  levels  of  anisoplanatic 
degradation. 

•  The  first  published  analysis  of  temporal  anisoplanatism  for  a  dynamic  system  [94], 
and  the  first  analysis  of  angular  anisoplanatism  in  adaptive  optics  for  interferomet¬ 
ric  imaging  [95],  establishing  the  fringe  Strehl  ratio  as  a  new  metric  to  quantify 
turbulence  degradation. 

•  A  new  class  of  tilt-sensing  techniques  for  laser  guide  star  adaptive  optics  to  obtain 
tilt  information  lost  due  to  reciprocal  compensation  in  the  LGS  signal. 

•  The  first  known  use  of  tilt  plus  higher-order  modal  measurements  to  improve  tilt 
compensation  performance  in  anisoplanatic  systems. 

•  A  simplified  system  design  for  adaptive  optics  in  interferometric  imaging,  where  an 
optimal  conjugate  phase  is  placed  on  a  single  deformable  mirror,  based  upon  wave- 
front  sensor  measurements  in  each  aperture. 

7.3  Summary  of  results 

The  application  of  optimal  modal  compensation  in  adaptive-optical  systems  whose 
compensation  performance  is  degraded  by  anisoplanatism  yielded  a  variety  of  results.  Key 
results  associated  with  each  of  the  system  types  considered  are  summarized  in  this  section. 

7.3.1  Ground-based  adaptive- optical  imaging  results. 

•  For  compensation  from  an  off-axis  NGS,  MMSE  estimation  resulted  in  compensation 
performance  improvement  for  all  orders  of  correction  and  at  all  beacon  offset  angles. 
For  the  correction  orders  AT  =  6, 10, 15, 21  with  aim  aperture,  enhanced  performance 
was  noted  for  beacon  offsets  greater  than  2  Oq.  For  N  =  28,36,45,55  with  aim 
aperture,  enhanced  performance  was  noted  for  beacon  offsets  greater  than  0.5  0o*  For 
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compensation  through  AT  =  55,  the  MMSE  estimator  can  provide  equivalent  imaging 
performance  at  NGS  offset  angles  10%  larger  than  when  NGS  measurements  are  used 
directly  for  compensation. 

•  For  compensation  from  an  on-axis  LGS  with  an  off-axis  NGS  for  wave-front  tilt 
measurement,  the  optimal  estimator  resulted  in  enhanced  compensation  performance 
for  all  LGS  altitudes  considered.  The  MMSE  estimator  partially  reduces  the  effect  of 
focal  anisoplanatism.  For  compensation  through  N  =  55,  the  MMSE  estimator  can 
provide  equivalent  imaging  performance  at  altitudes  20%  lower  than  when  beacon 
measurements  are  using  directly  for  wave-front  compensation. 

•  For  an  off-axis  NGS,  including  higher-order  modal  measurements  in  the  tilt  estimator 
resulted  in  increased  Strehl  ratio  performance  over  using  tilt  measurements  only. 
For  tilt  compensation  using  beacon  Zernike  modes  2  to  15,  it  was  shown  that  the 
maximum  beacon  offset  angle  for  nearly-limiting  performance  is  3  times  larger  than 
when  only  tilt  measurements  are  used.  For  a  given  beacon  offset,  the  Strehl  ratio  is 
increased  by  incorporating  higher-order  modes. 

•  Wave-front  tilt  cannot  be  measured  directly  from  the  LGS  signal  due  to  natural  com¬ 
pensation  caused  by  beam  reciprocity.  Nonetheless,  optimal  compensation  methods 
permit  on-axis  wave-front  tilt  components  to  be  estimated  using  only  higher-order 
measurements  from  off-axis  LGS  beacons.  It  was  shown  that  for  a  constellation  of 
four  LGS  beacons  offset  orthogonally  from  the  optical  axis,  using  measurements  of 
Zernike  modes  4  to  15  from  each  beacon,  the  Strehl  ratio  can  be  greater  than  0.3  for 
D/tq  =  4  and  Lq/D  =  10.  Tilt  estimation  performance  is  degraded  by  focal  aniso¬ 
planatism.  Peak  LGS  tilt  estimation  performance  occurs  at  an  angle  which  depends 
upon  the  diameter  of  the  aperture  and  the  height  of  the  effective  turbulence  layer. 
The  peak  angle  can  be  significantly  larger  than  the  isoplanatic  angle  for  typical  AO 
systems. 
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1.3,2  Airborne  Laser  beam- steering  adaptive  optics  results. 


•  Due  to  anisoplanatism,  tilt  correction  is  suboptimal  when  time-delayed  tilt  mea¬ 
surements  are  used  directly  in  the  beam-steering  adaptive  optics  for  the  Airborne 
Laser  system.  Optimal  correction  is  obtained  using  the  MMSE  estimator  given  time- 
delayed  beacon  measurements,  which  may  include  tilt  plus  higher-order  modes. 

•  While  tilt  is  uncorrelated  with  most  higher-order  modes  at  r  =  0,  the  correlation 
passes  through  some  extremum  at  r  >  0  before  going  to  0  at  large  values  of  r.  Thus, 
incorporating  higher-order  modes  in  the  tilt  estimator  will  enhance  tilt  correction 
performance. 

•  For  the  ABL  target  engagement  scenario,  enhanced  tilt  compensation  performance 
is  obtained  by  including  delayed  modal  measurements  through  Zernike  15.  Thus, 
compensation  performance  may  be  maintained  to  within  1%  of  its  limiting  value  at 
longer  measurement  time  delays,  resulting  in  a  reduced  AO  bandwidth  for  the  ABL. 
It  was  shown  that  by  including  beacon  measurements  through  Zernike  15,  the  beam¬ 
steering  AO  bandwidth  is  reduced  by  over  70%.  At  a  fixed  bandwidth,  the  Strehl 
ratio  is  increased  by  incorporating  higher-order  modes  into  the  MMSE  estimator. 

7.3.3  Adaptive- optical  interferometry  results. 

•  The  angular  anisoplanatic  effect  was  characterized  for  adaptive  optics  in  interferomet¬ 
ric  imaging  using  modal  analysis  techniques.  For  D/ro  =  5,  the  beacon  angular  limit 
is  between  1.5  Oq  and  3  9q^  depending  upon  the  compensation  order.  For  D/vq  =  30, 
the  beacon  angular  limit  is  approximately  0.2^o  with  compensation  through  N  =  153. 
This  analysis  demonstrates  that  performance  degradation  due  to  angular  anisopla- 
natism  is  more  severe  for  interferometric  applications  than  for  imaging  applications. 

•  For  optimal  compensation  in  an  interferometer,  beacon  measurements  should  not 
used  directly  to  conjugate  the  phase  in  each  aperture,  but  instead,  an  optimal  phase 
should  be  formed  on  a  deformable  mirror  in  one  aperture  based  upon  the  difference 
of  beacon  measurements  in  each  aperture. 
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•  Application  of  an  optimal  conjugate  phase  results  in  an  increased  fringe  Strehl  ratio, 
quantified  by  the  FSR  gain,  GpgR-  The  value  of  Gpgj^  depends  upon  the  compen¬ 
sation  order  iV,  the  relative  baseline  s/D^  and  the  relative  aperture  diameter  D/tq. 
For  =  21  at  the  beacon  offset  =  10?  optimal  compensation  resulted  in  an 

FSR  increase  of  15%  to  over  25%,  depending  upon  the  relative  baseline  s/D. 

Conclusions  drawn  from  research 

Although  anisoplanatism  is  manifested  in  diverse  ways  for  different  types  of  adaptive- 
optical  systems,  its  physical  origin  is  the  same.  Measured  wave-front  aberrations  and 
the  turbulence-induced  phase  from  the  object  of  interest  are  not  conjugate,  but  instead 
only  partially  correlated  in  anisoplanatic  systems.  This  effect  results  from  differences  in 
the  turbulence  volume  sampled  by  beacon  and  object  light.  This  generalized  view  of 
anisoplanatism  serves  as  a  unifying  concept  in  the  analysis  of  a  broad  class  of  adaptive 
optics  for  imaging,  beam-projection,  and  interferometric  systems. 

In  the  context  of  modal  analysis  techniques,  wave-front  phase  aberrations  are  de¬ 
scribed  by  projecting  the  turbulence-induced  phase  onto  a  basis  set  of  Zernike  polynomials. 
Thus,  anisoplanatism  is  quantified  by  the  spatial  cross-correlation  of  Zernike  coefficients. 
By  considering  a  mathematical  description  of  modal  compensation  performance,  optimal 
estimation  techniques  can  be  applied  to  adaptive  optics.  The  resulting  optimal  modal 
compensation  theory  indicates  that  when  the  performance  of  an  AO  system  is  degraded 
by  anisoplanatism,  the  conjugate  phase  should  be  estimated  from  the  wave-front  measure¬ 
ments.  Otherwise,  direct  application  of  the  wave-front  measurements  in  phase  conjugation 
will  result  in  suboptimal  compensation  performance. 

The  application  of  optimal  modal  compensation  results  in  enhanced  performance  for 
ground-based  AO  imaging  systems,  dynamic  beam-projection  systems  like  the  Airborne 
Laser,  and  adaptive-optical  interferometers  for  high-resolution  image  reconstruction.  For 
ground-based  imaging  systems,  natural  guide  star  beacons  at  larger  angular  offsets  and 
lower-altitude  laser  guide  star  beacons  may  be  employed  without  significantly  reducing 
system  performance.  In  addition,  wave-front  tilt  estimates  can  be  obtained  from  a  laser 
guide  star  even  though  tilt  cannot  be  measured  directly  from  the  LGS  signal.  For  beam- 
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steering  adaptive  optics  in  the  Airborne  Laser,  optimal  tilt  compensation  permits  longer 
measurement  time  delays,  thus  reducing  the  bandwidth  of  such  a  system.  When  adaptive 
optics  are  used  with  a  spatial  interferometer,  optimal  compensation  results  in  increased 
fringe  visibility  at  larger  beacon  offsets,  reducing  the  degrading  effects  of  angular  aniso- 
planatism,  and  leading  to  better  image  reconstruction. 

7.5  Recommendations  for  future  research 

The  research  activity  presented  in  this  dissertation  presents  a  description  of  optimal 
compensation  for  anisoplanatism  that  is  complete  and  consistent.  While  the  research 
answers  many  important  questions,  it  asks  just  as  many,  if  not  more.  There  were  several 
considerations  that  were  outside  the  scope  of  the  current  research  effort  that  could  be 
addressed  in  future  research  activity.  These  issues  are  outlined  in  the  following  paragraphs. 

Turbulence  modeling.  The  generalized  concept  of  anisoplanatism  espoused  in  this 
research  focused  on  aperture  and  source  spatial  separations  as  the  primary  source  of  aniso¬ 
planatism.  Indeed,  this  view  is  consistent  with  all  previous  conventional  analysis  efforts 
that  assume  statistical  isotropy  in  turbulence  modeling.  It  is  recognized,  however,  that  sta¬ 
tistical  anisotropy  will  change  the  correlative  properties  of  all  Zernike  modes,  especially  the 
cross-correlation  between  different  aberration  orders.  In  this  sense,  statistical  anisotropy 
would  give  rise  to  a  performance  degradation  effect  in  AO  systems  similar  to  conven¬ 
tional  or  the  generalized  form  of  spatial  anisoplanatism.  By  invoking  reasonable  models 
for  anisotropic  turbulence,  the  cross-correlation  calculations  for  Zernike  modes  could  be 
extended  to  account  for  this  effect.  As  a  result,  the  optimal  compensation  theory  would 
help  to  mitigate  the  performance  degradation  resulting  from  anisotropic  turbulence. 

Measurement  noise  effects.  In  the  development  of  the  optimal  compensation  theory, 
it  was  assumed  that  apart  from  anisoplanatism,  the  wave-front  sensing  device  was  ideal. 
Therefore,  the  effects  of  sensor  noise  were  not  considered.  If  an  appropriate  noise  model 
for  the  Zernike  coefficients  resulting  from  wave-front  sensor  measurements  were  developed, 
that  model  could  be  incorporated  into  the  MMSE  estimator,  resulting  in  a  reduced  sensi¬ 
tivity  to  noise  effects  caused  by  low  light  levels  or  detector  read  noise. 
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Laser  guide  star  tilt  sensing.  The  method  for  wave-front  tilt  sensing  developed  in 
this  research  is  innovative  in  that  it  takes  advantage  of  the  natural  properties  of  aniso- 
planatism  to  obtain  a  tilt  measurement  indirectly  from  higher-order  modal  measurements. 
The  particular  off-axis  LGS  geometry  investigated  here  is  just  one  of  many  ways  to  utilize 
the  cross-correlation  between  tilt  and  higher-order  modes.  Even  if  a  single  LGS  is  de¬ 
ployed  on-axis,  higher-order  wave-front  measurements  from  spatially-offset  subregions  of 
the  aperture  could  be  employed  to  estimate  the  full-aperture  tilt.  Although  the  LGS  signal 
does  not  contain  full-aperture  tilt  components,  it  does  contain  local  wave-front  gradients 
that  correspond  to  sub-aperture  tilts.  Indeed  these  subaperture  tilts  are  the  fundamen¬ 
tal  signal  used  to  reconstruct  higher-order  modes.  It  follows  therefore  that  wave-front 
sensor  subaperture  tilts  should  be  correlated  with  the  full- aperture  tilt.  In  this  case,  an 
optimal  estimate  of  the  full-aperture  tilt  could  be  formed  on  the  basis  of  subaperture  tilt 
measurements. 

Higher-order  compensation  for  ABL,  The  application  of  optimal  compensation  to 
adaptive  optics  for  the  Airborne  Laser  dealt  only  with  beam-steering  or  tilt  correction. 
While  tilt  compensation  is  the  most  important  aberration  correction  in  this  system,  higher- 
order  compensation  may  be  required,  depending  upon  the  relative  aperture  diameter  D/r^. 
In  this  case,  optimal  higher-order  modal  compensation  would  also  result  in  increased  Strehl 
ratio  for  the  irradiance-on-target  with  extended  propagation  paths.  The  required  modal 
cross-correlations  can  be  obtained  from  the  current  expressions,  with  proper  accounting 
for  the  effects  of  aperture  and  source  motion. 

Laser  guide  star  in  interferometry.  When  considering  the  effect  of  optimal  compen¬ 
sation  in  adaptive  optics  for  interferometric  applications,  the  only  effects  investigated  in 
this  research  were  those  associated  with  angular  offsets  of  a  natural  guide  star.  Since  the 
limits  of  angular  anisoplanatism  in  interferometry  are  more  severe  than  in  conventional 
imaging,  it  is  reasonable  to  assume  that  laser-guide-star  adaptive  optics  may  be  required 
for  reasonable  sky  coverage.  In  this  case,  focal  anisoplanatism  also  will  affect  compensation 
performance.  As  was  the  case  in  conventional  single-aperture  imaging  with  LGS  compen¬ 
sation,  there  will  be  limits  placed  on  the  beacon  altitude  by  focal  anisoplanatism.  These 
limits,  and  the  ability  for  optimal  compensation  to  expand  these  limits  could  be  investi- 
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gated  using  the  current  theory  and  cross-correlation  expressions,  with  proper  modeling  to 
account  for  the  finite  altitude  of  the  LGS  beacon. 
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Appendix  A.  Conversion  of  spatial  integral  to  a  spectral  integral  in  the  Zernike 
coefficient  cross- correlation  calculation 

In  this  appendix,  the  conversion  of  the  spatial  integral  for  in  Eq.  (145)  to 

the  spectral  integral  in  Eq.  (146)  is  detailed.  This  development  begins  with  the  expression 
for  the  spatial  cross-correlation  given  in  Eq.  (145): 


■®Ouiiav2j  (^»tt5  ^02) 
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For  spatial  functions  /(f)  and  5(f),  Parseval’s  theorem  [32]  says  that: 


J  df /(r)/(f)  =  j  df  F{f)G*{f),  (204) 

where  F{f)  and  G{f)  are  the  Fourier  transforms  of  /(f)  and  g{r),  respectively.  The  scaling 
properties  of  Fourier  transforms  applied  to  the  spatial  function  /(r)  may  be  written  as  [30]: 

^t{/(^)}  =  ]6]F(6/),  (205) 

where  FT{-}  denotes  the  Fourier  transform  operation.  Applying  Parseval’s  theorem  and 
the  scaling  properties  of  the  Fourier  transform  to  the  integral  over  q2i  in  Eq.  (203)  yields: 


Rdunavij  (^««)  ^ol)  ^02) 

=  ^[il(l  —  A„ij)]“‘  J  dg„ij  J  die  FT  {5^(g„2j  +  Si  -  g„ij)} 

where  Qj{f)  is  the  Fourier  transform  of  Zj{f)W{r). 

The  Wiener-Khinchin  theorem  says  that  for  a  spatially  wide-sense  stationary  quantity 
x(f),  the  spatial  power  spectral  density  of  x,  designated  Wx{k)  is  the  Fourier  transform  of 
the  spatial  autocorrelation  of  x  [32],  designated  Bx{^r).  Since  the  phase  in  the  Ith  layer 
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is  spatially  wide-sense  stationary,  then,  by  the  Wiener-Khinchin  theorem: 


TT  {B^iqv2t  +  si-  g„i,)}  =  W^{k,  zi)  exp  j2Trk  •  {quu  -  «/)]  , 


(207) 


where  the  exponential  term  arises  from  the  spatial  shift  property  of  the  Fourier  trans¬ 
form  [30].  Thus,  substituting  Eq.  (207)  into  Eq.  (206)  yields: 

^auUO.v2j  ^ai;  f*SV  5  ra2) 

=  j  j  dk  w^{k,  Zi)  exp  [j27rfc  •  (g„i,  -  s,)] 


xZi 


Quil 


w 


Quit 


Q*  [i2(l  -  A„2,)fc]  . 


(208) 


Exchanging  the  order  of  integration  over  q^u  and  fc,  and  rearranging  terms  in  the  integral 
yields  the  following  expression: 


^0>uHO'v2j  (^5W  5  TaU  ra2) 


=  [Ril  -  Aui,)]-^  I  dk  W^ik,  z,)Q*  [/?(1  -  A,2,)k]  exp(-j27rfc  •  5,) 


J  dquii^i 


Quit 


W 


i2(l  — Aui()J  [-^(1  — -Aui/) 


quit 


exp(j27rfc  •  g„i,). 


(209) 


At  this  point,  it  is  recognized  that  the  integral  over  q^u  in  Eq.  (209)  is  the  Fourier  transform 
of  Zi[  ■  ]W[  •  ].  Since,  by  the  scaling  property  of  Fourier  transforms. 


TT 


Wi 


quit 


w 


9uiJ 


LiJ(l-A„w)J  Li2(l-A„H) 

it  is  therefore  concluded  from  Eq.  (209)  that 


I  =  ii(l  -  A^u)Qi  [i?(l  -  A„i,)^]  ,  (210) 


^Q‘uliO'v7j  (XsU^  Tal,  ra2) 


=  dkW^{k,z,)exp{-j2'Kk  ■  si)Qi  R{1  -  Auu)k  Qj  ^R{1  -  Ay2i)k  ,  (211) 


which  is  the  spectral  integral  of  Eq.  (146). 
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Appendix  B.  Integration  of  the  azimuthal  component  of  the  Zernike  coefficient 

cross- correlation  expression 

In  this  appendix,  the  integration  of  Eq.  (146)  with  respect  to  the  azimuthal  compo¬ 
nent  of  k  is  performed  to  yield  the  expression  in  Eq.  (147),  which  involves  one-dimensional 
integration  over  the  radial  component  of  k.  The  development  begins  by  introducing  the  fol¬ 
lowing  expression  for  the  Fourier  transform  of  the  product  of  the  ith  (Noll’s  ordering  [54]) 
Zernike  polynomial  and  the  aperture  weighting  function,  as  a  function  of  the 

radial  and  azimuthal  components  of  k,  designated  k  and  0,  respectively: 

{™i»  +  f(l  -  -  1|}  .  (212) 

where  Ui  and  are  the  radial  and  azimuthal  orders  of  the  ith  Zernike  polynomial,  Jk(-)  is 
the  Bessel  function  of  the  first  kind  of  order  k,  and  6  is  the  Kronecker  delta.  The  expression 
for  9)  in  Eq.  (212)  differs  from  its  conventional  form,  but  is  mathematically  equivalent 
to  the  multi-case  expression  given  by  Noll  [54]. 

Since  Ogj  is  defined  as  the  angle  between  the  vector  si  and  the  unit  vector  x,  and  if 
k  is  taken  to  be  a  vector  whose  angle  with  respect  to  x  is  0,  then  it  may  be  said  that: 

k  ^  Si  =  ksi  cos{9  —  05,).  (213) 

Hence  Eq.  (146)  may  be  written  as: 


^ol)  ^02) 

rOO  r03j+2lT 

=  dkkW^{k,Zi)  J  d0  exp[-j27rA:sjcos(^  -  0,,)] 

xQi  [i?(l  -  Auu)k,  0]  Q*  [il(l  -  A,2,)k,  0] . 


(214) 


Substituting  the  expression  of  Eq.  (212)  for  Qi  and  Qj  into  Eq.  (214)  yields  the  following 
expression: 
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•®aui»Ou2j  ^al )  ^a2) 

=  ^[7r^i?^(l  —  >l.uij)(l  “  ■Ai>2()]~^  [(^i  +  +  1)]  ^ 

I 

rOO 

X  /  dfc  fc-i  Ty^(A:,  ^,)  J„,+i  [27rJ?(l  -  A^u)k]  J„,+i  [27ri?(l  - 
Jo 

/^jj“}"27r  f  TT  *  1 

d0  exp[— j27rA:sj  cos{9  -  ^j,)]  cos  yriiO  +  —(1  -  ^mjo)[(-l)*  ~ 

X  cos  i^rujO  +  ^(1  -  ^TO,o)[(“l)^  ~ 


The  product  of  cosine  factors  in  Eq.  (215)  can  be  expanded  as  follows: 

cos  jmifl  +  ^(1  -  ^mio)[(-l)*  -  1]|  cos  ^nijO  +  ^(1  -  ^m^.o)[(-l)^  -  1]| 

=  i  cos  [(rrij  +  mj)$  +  5]  +  ^  cos  [{mi  -  mj)0  +  h] ,  (216) 

z  z 

where  g  and  h  are  defined  as  follows: 

g  =  I{(l-««.)l(-l)‘-l|+(l-«m,o)l(-l)'-ll}.  (217) 

h  =  j{(l-«„,.)|(-l)‘-l)-(l-«:,,.)l(-l)l-l|}.  (218) 

If  we  now  make  the  substitution  (f)  =  0  -  6,,,  then  the  right  hand  side  of  Eq.  (216)  may 
then  be  expanded  as  follows: 


^  cos  [(m<  +  mj)0  +  s]  +  ^  cos  [(m^  —  mj)0  +  h] 
z  z 

f 

I  cos  [(mi  +  mj)(f>]  cos  [(mj  +  mj)6g,  +  g] 
-i  sin  [(mi  +  mj)(f>]  sin  [(mi  +  mj)0„  +  g] 

=  < 

+1  cos  [(mi  —  mj)(f)]  cos  [(mi  —  mj)0j,  +  /i] 
sin  [(mi  -  mj)(t>]  sin  [(mi  -  mj)9s,  +  h] . 


(219) 
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When  this  expansion  of  the  cosine  terms  is  substituted  into  Eq.  (215),  the  integration  with 
respect  to  6,  designated  here  as  I0  becomes  (using  <f>  =  6  —  6g,): 


r27r 

l0=  d<^  exp[— j27rA:sj  cos 

Jo 

I  cos  [{rrii  +  mj)(l)]  cos  [(mj  +  mj)6s,  +  g] 
sin  [{nii  +  mj)<f>]  sin  [(mj  +  +  g] 

+\  cos  [{rrii  -  mj)4>]  cos  [(m^  -  mj)0,^  +  h] 
[  -5  sin  [(mi  -  mj)4>]  sin  [{rrii  -  mj)0f,  +  h] . 


(220) 


Now,  using  the  following  definition  for  the  Bessel  function  of  the  first  kind  of  order  k  [1]: 


1  r 

Jk{z)  =  — K  /  d<^  exp(j2cos(;6)cos(/c^), 
V— 1  TT  Jo 


(221) 


it  can  be  shown  that: 


/  d(/>  exp 3 z  cos  (j))  cos («;</>) 

Jo 

p27r 

/  d(^  exp(— j2cos(^)sin(/c^) 

Jo 


=  27r(-l)2''J„(2:), 


=  0. 


(222) 

(223) 


Hence,  I0  in  Eq.  (220)  is  evaluated  as: 


I0  =  -IT 


(-l)2('"*'^’”^)*^(mi+m,)(27rfcs,)cos[(mi  +  mj)0„  +3] 

_^^l(27r/cs,)cos[(mi  -mj)0„  +  h]. 


(224) 


Substituting  Eq.  (224)  in  for  the  integral  over  6  in  Eq.  (215),  as  well  as  the  definitions  of  the 
quantities  g  and  h  defined  in  Eq.  (217)  and  Eq.  (218),  we  obtain  the  following  expression 


for  B, 
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^CLuli  ttv2j  8U  3  ^Ol  3  3  fa2) 

TT  Y 


-1 


/  X  \ 

X«^(ni+1)[-R(l  “  ^ull  )27rfc]J(„j,+i)[jR(l 

/°°dk  (  X  \ 

^•^(ni+l)[-^(l  ~  ■^ul/)27rfc]«/(nj+l)[-^(l 


A^2i)'2'!rk] 

A^2i)2nk], 


(225) 


where  the  functions  F,  G,  and  H  are  defined  as: 


F  =  [(nj  +  l)(n_,' +  (226) 

a  =  (-l)t('™i+™i)cos[(mi  +  mj)0„ 

+7r/4  {(1  -  <5„.,o)[(-l)‘  -  1]  +  (1  -  «5„.,o)[(-l)^'  -  1]}]  ,  (227) 

H  =  (-l)tl”*^-'"ilcos[(mi-m^)0„ 

+7r/4{(l  -  5„,,o)[(-l)‘  -  1]  -  (1  -  ^m,o)[(-l)^’  -  1]}]  .  (228) 


By  defining  x  =  irDk^  where  D  =  2i?,  Eq.  (225)  may  be  rewritten  as: 


^(luUO.v2j  (Xsu^  ^al  5  ^3t;3  *^02) 

4 

7r 


=  -■£>-"-F’5][(i-a„h)(i-a„2,)]“' 


X  < 


X-^(n.+l)[(l-  ■^ull)^]*7(nj+l)[(l  ■^v2l^^]i 


(229) 


which  is  the  cross-correlation  expression  given  in  Eq.  (147). 
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Appendix  C.  Spatial  cross-correlation  of  Zernike  coefficients  for  apertures  of 

different  diameter 

The  development  presented  in  Chapter  III  for  the  spatial  cross-correlation  of  Zernike 
coefficients  assumed  that  the  two  apertures  in  the  analysis  geometry  of  Figure  12  were 
of  equal  radius  R,  This  analysis  may  be  modified  to  yield  a  more-general  expression  for 
the  spatial  cross-correlation  in  different-sized  apertures.  If  aperture  1  has  radius  and 

aperture  2  has  radius  i?2,  then  the  Zernike  coefficients  and  ay2j  are  given  by: 


<^uii  —  j  dpi<i)y,i{Ripi)Zi{pi)W{pi) 

(230) 

Q'v^j  =  [  dp2  (t>v2{R2p2)Zj{p2)W{p2)- 

(231) 

The  projected  aperture  vectors  qy^n  and  qy2i  are  modified  accordingly: 

quii  =  (1  —  Auu)RiP\  (232) 

?v2l  =  (1  —  A>^J2l)R2p2^  (233) 

where  the  layer  scaling  factors  A^u  and  Ay2i  remain  unchanged  from  the  definitions  given  in 
Eq.  (118)  and  Eq.  (119),  respectively.  If  Eq.  (232)  and  Eq.  (233)  are  employed  throughout 
the  development  in  Chapter  III  and  Appendix  A-B,  the  steps  in  the  development  are  the 
same,  with  the  quantity  R{1  -  Auu)  replaced  by  i?i(l  —  Auu)  and  R{1  —  Ay2i)  replaced  by 
ii2(l  A^)2i)*  When  the  spectral  variable  is  transformed  from  A:  to  x  in  Appendix  B,  the 
appropriate  transformation  is  x  =  TrDik,  where  the  diameter  of  aperture  1  is  JDi  =  2i?i. 
Thus,  grouping  Di  with  ro,  and  accounting  for  the  variable  transformation,  the  expression 
for  different-sized  apertures  corresponding  to  Eq.  (152)  is  given  by: 
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! 


•®Oul»Ov2i  (X 8u  j  ^ol  j  >  ^a2  ) 

G  /  [^a;] 

«/o 

X  \  ^oo  -  ^ 

Jo 

^  ><*^(ni+l)[(l  -  Aull  )j^]J(nj+i)  [;^(^  ^t;20^]  ) 


(234) 


where  D2  =  2R2.  Thus,  the  cross-correlation  is  normalized  for  the  size  of  aperture  1  rela¬ 
tive  to  the  atmospheric  coherence  diameter,  Di/ro,  and  depends  upon  diametric-disparity 
parameter,  D2/D1. 
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